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ABSTRACT 


The thermal conductivity of the solidified inert gases A, Ne and Kr, 
contained in an Inconel tube, has been measured at temperatures from 2°x 
to near their melting points. Above about 15°K, we find that the thermal 
conductivity, 4, may be represented by the expression AY B/T', where B has 
values of 250 mw em! for A, 75 for Ne, and 270 for Kr. Values calcu- 
lated from approximate theoretical treatments of lattice anharmonicity 
yield values in excess of these by about 30% to 80%. At liquid helium 
temperatures, i.e. below the conductivity maximum, we find that X varies 
approximately as 7? and is small in magnitude, suggesting mean free paths 
of 10-* to 10-*cem for the lattice waves. This may result from the large 
contraction of the solids during cooling which produces severe internal strains 
and probably results in many lattice defects within the crystals. 


$1. INTRODUCTION 


In the absence of free electrons heat is transported through a solid by lattice 
waves. We may write for the conductivity, 2: 


A~4Cvl 


where C is the specific heat per unit volume, v is the velocity of lattice 
waves (i.e. essentially the velocity of sound) and / is their mean free path. 
Static imperfections of the lattice such as impurity atoms or physical 
defects will cause scattering, but even in their absence the mean free 
path will remain finite because of the coupling between the normal modes 
which arises from the anharmonic component of the interatomic forces. 
Many years ago Debye (1914) suggested that at moderately high tempera- 
tures (7 $6), where the specific heat is sensibly constant, the anharmonic 
coupling should cause scattering proportional to the absolute temperature, 
i.e. Acc1/T7. Later Peierls (1929) examined the problem more rigorously 
and showed that only a particular class of anharmonic interactions would 
cause resistance to heat flow. These Umklapp processes, as Peierls 
called them, are quite frequent at ‘high’ temperatures and give rise to a 
thermal resistance proportional to T (i.e. Acc1/7'). At low temperatures, 
however, it becomes increasingly difficult to satisfy the conditions neces- 
sary for Umklapp processes and the probability of their occurrence falls 
off rapidly as exp (—9/bT) whereb~1. Consequently at low temperatures 
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the resistance due to anharmonicity diminishes rapidly, the mean free 
path increasing roughly as exp (6/bT). Ultimately lattice imperfections 
and crystal boundaries will set a limit to the mean free path. These 
different scattering mechanisms have been discussed in detail in the recent 
reviews by Berman (1953) and Klemens (1956), both in relation to the 
theoretical treatment and the results of experimental studies. 

Although Peierls deduced no absolute magnitudes for the resistance 
due to anharmonic coupling, attempts have been made in the past few 
years to predict this theoretically for simplified models of solids. These 
attempts have related \ (in the higher temperature region where A 1/7’) 
to the anharmonicity through such other physical properties as the 
specific heat, C (or Debye characteristic temperature, 6,), the thermal 
expansion coefficient, «, the compressibility, y, and the Gruneisen para- 
meter, y=a/yC. For the alkali halides fair agreement seems to be obtained 
between theory and experiment; for more anisotropic materials like 
Ge, Si, Bi, Sb and Te agreement is less satisfactory ; the theoretical treat- 
ments appear to over-estimate the magnitude of AT’ by factors of 4 or more. 
This is not too surprising because the magnitude of the anharmonic 
coupling is derived for linear chain models or simple cubic crystals with 
nearest neighbour interaction of a central force character. 

It seemed desirable to us to determine experimentally the heat con- 
ductivity of some of the solidified inert gases which crystallize in a face- 
centred cubic system; for these the interaction forces can be represented 
by a simple radial potential of short-range. Apart from solid helium 
which is close-packed hexagonal in structure and has been investigated 
extensively by Wilks and his collaborators (Wilkinson and Wilks 1951, 
Webb e¢ al. 1952, Webb and Wilks 1953), there appeared to be no available 
thermal conductivity data on the inert gas solids. A preliminary report 
of this present work was published (White and Woods 1956) when only 
some results on (solid) argon were available. Our experiments have not 
yet been extended to xenon, because the gas is very expensive. 


§ 2, EXPERIMENTAL METHOD 


The cryostat used by us has been described previously (White and 
Woods 1955). For these measurements a thin walled Inconel specimen 
tube, S, was mounted in the cryostat (see fig. 1) and could be filled with 
gas through the tube, T. The inert gas under investigation was supplied 
from a small cylinder through a reducing valve and needle valve to T. 
Swas sealed with copper end-pieces silver-soldered to it and had two copper 
wires (about 0-030 in. diameter) also silver-soldered through the side wall 
to act as connections for the thermometers, A and B. The temperature 
difference, AT’ =‘7', —7',, is produced by a measured amount of electrical 
power dissipated in the heater, H,, which is attached to one end of the 
specimen. The temperature of the other end of the specimen is controlled. 
together with chamber C and the inner shield, by pumping on a liquid 
refrigerant in C or by electrical heating in H, regulated automatically. 
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To fill S with solid, a suitable liquid refrigerant was placed in the dewar 
surrounding the cryostat, gas was fed to T at about atmospheric pressure, 
and a small pressure (~1 mm) of helium gas was kept in the cryostat to 
allow heat transfer from 8 to the refrigerant. When condensing argon 
into 8, the cryostat was covered with liquid nitrogen to a height indicated 
by the arrow 1. With neon, the dewar was filled with liquid hydrogen 


Fig. 1 


Schematic diagram of cryostat containing specimen tube. 


to a similar level (1). When condensing krypton, liquid oxygen covered. 
the cryostat to a level indicated by arrow 2. Since the filling tube T ee 
never in direct contact with the refrigerant and warm gas was led into § 
through T continually until S was filled, blocking of T by solid was usually 
avoided. At the conclusion of an experiment, the gas in 8 was evaporated 


3H2 
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and pumped away through a flow-meter so that a good estimate of the 
volume of gas liquefied (and thence solidified) could be made. 

We had hoped to reproduce the conditions of growth of solid used by 
Stansfield and Followell (1956) to grow solid argon. Although we were 
unable to watch the growth in the Inconel tube, some trial experiments 
were done with argon in an all-glass cryostat, i.e. a glass tube in a glass 
enclosure surrounded by a dewar of liquid nitrogen. In this case liquid 
gradually collected in 8, so that the level rose at about 1 or 2 mm/min ; 
the liquid soon began to solidify from the bottom so that the solid—liquid 
interface also rose at about the same rate but was a few mm below the 
free liquid surface. Although a few ‘vapour-snakes’ (Stansfield and 
Followell 1956) appeared in the solid just before it completely filled the 
tube, the solid mass soon became quite transparent. However, on 
replacing the liquid nitrogen by liquid helium and cooling the specimen 
by means of exchange gas, the solid became quite cloudy and opaque; it 
only regained transparency on rewarming and being left for a few minutes 
above about 60°k. Although the rate of cooling with liquid helium was 
slower and more easily controlled in the metal cryostat it seems quite 
likely that the behaviour of the specimens in an actual experiment was 
similar. At least the thermai conductivity results obtained below 10°xK, 
which are among those presented in §3, suggest strongly that the solids 
were far from perfect crystals at these low temperatures. In view of the 
large expansion coefficients of these solids, and the constraint imposed by 
the Inconel tube, this is probably not surprising. 

After performing a number of experiments on argon in a 4 in. diameter 
Inconel tube of about 3 in. total length with thermometer leads about 
2in. apart which were soldered right across the tube as shown in fig. 1, 
a 3 in. diameter Inconel tube about 2 in. long was put in its place. This 
provided greater economy in use of gas, and also the two copper leads 
were soldered to only one wall in an attempt to reduce the constraints 
somewhat. If the sample was indeed less constrained this fact is not 
reflected in the thermal conductivity results. 

The sources of gas used for the various specimens are listed below: 


Al, A2: Commercial purity argon from Linde Air Products 
Company (New York). 

A3, A4: High purity argon from Air Reduction Company (New 
Jersey). 

A5, A6: High purity argon from Matheson Company (New 
Jersey). 

Nel: High purity neon from Linde Air Products. 

Ne2, Ne3: High purity neon from Air Reduction Company. 


Krl, Kr2, Kr3: High purity krypton from Linde Air Products. 


The purity of Ne and Kr from Linde was quoted as 99-96% by the supplier, 
and we believe the other high purity gases which were supplied in small 
steel cylinders also contained less than 0-1° contaminants. 
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§ 3. REesuLts 


3.1. Corrections for the Thermal Conductivity of the Inconel Container 

A thin-walled Inconel tube (Inconel is an alloy of 14-17% Cr, 6-10% Fe, 
~0-1% C and remainder Ni) was used for the specimen tube, 8, because of 
its low thermal conductivity. The tubing, of about 0-006 in. wall thickness, 
was supplied in the hard-drawn condition but was partially annealed as a 
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Estermann and Zimmerman (1952). 


result of soldering operations. From the measurements of Estermann 
and Zimmerman (1952), data for its thermal conductivity in both the 
hard-drawn and annealed conditions were already available (see fig. 2). 
However to ensure that entirely reliable corrections could be made, we 
also measured the thermal conductivity of an Inconel specimen which was 
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in a similar physical condition to our sample tube. Because of the very 
low thermal conductivity of the material, experiments to try and obtain 
values directly for a thin-walled tube, particularly at the lowest tempera- 
tures, were both slow and rather unreliable. A specimen was therefore 
made by cutting about 20 strips, each about 1 cm wide x 7 cm long, from 
the thin-walled tube, placing them together to form a bar and silver- 
soldering end-contacts and potential leads to them. The thermal con- 
ductivity of this specimen (cf. fig. 2) was found to be generally intermediate 
between the values obtained on annealed and hard-drawn Inconel by 
Estermann and Zimmerman. 

When the tube, §, is filled with a solidified inert gas, the total heat 
conducted, per unit time, is given by 


. A 
doa), aton(t, 9 


where Ay, A, are the heat conductivities of the solidified gas and Inconel 
respectively, and (A/L), and (A/L), are their respective geometrical 
form factors (ratio of cross-sectional area to length). Thus in the 
absence of any correction for radiation loss, the true conductivity, A,, is 
related to the uncorrected value, 


ee 
7 ATA 


Ag =Ag* —A,(L/A)g/(L/ A), =Agt — Aa. 


by 


In our measurements A; and A, were generally comparable at temperatures 
of 10 to 20°K but (L/A)yg~4cem and (L/A);~70cm™ so that the 
correction is small below 20°K. At temperatures of 55° to 90°K, where a 
few rather laborious measurements on solid argon and solid krypton were 
made, the correction becomes large and a further small correction due to 
radiation has to be considered. This radiation correction had previously 
been determined in our cryostat, using a glass specimen of known thermal 
conductivity. As a further check on the accuracy of these corrections the 
effective thermal conductivity of the empty Inconel tube was determined 
at 77° and 90°K; the values obtained thus take into account both conduc- 
tion through the Inconel and the heat loss by radiation. 

For example for krypton at 78°K, \,*~ 10-7 mw cm- deg! and the 
correction is AA~7-1 mw em! deg-!, so that 


Ag~@10-7-—7-1 
= 3-6 mw cm~! deg-1. 
It may be concluded that results in this high temperature region are of 
only very limited accuracy, particularly as many hours were required 


for the temperature of the specimen and of the thermometers to reach 
equilibrium. 
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3.2. Argon 

The results shown in fig. 3 appear to be reproducible above 10°K and 
have the form to be expected if anharmonic interactions between lattice 
vibrations are the main cause of thermal resistance, that is \ is proportional 
to 7-'. Impurities might also simulate this 1/7 behaviour, but in view 
of the reproducibility between samples made from different cylinders of 
gas, it appears that the measured conductivity in this higher temperature 
tange is characteristic of pure argon. 


Fig. 3 
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At low temperatures, A falls approximately as 7? and is disappointingly 
small; application of the formula A= Cvl/3 shows that the mean free path 
of the lattice waves is only about 10-* cm at 2°K. If boundary scattering 
in large crystals was limiting the mean free path, then we should expect 
AccT?, and that the magnitude of the thermal conductivity would be 
much greater (cf. results of measurements on diamond, sapphire, Ge, Si, 
Te, Bi, solid helium, reviewed by Berman 1953, Klemens 1956). 5 

Since the density of solid argon increases by nearly 10% from 1-63 g cm™ 
at 80°K to about 1:77 g cm~? at 4°K (see Dobbs and Jones 1957) and the 
solid is constrained in an Inconel tube whose overall contraction is much 
less we might expect the solid to be in a highly disordered state at very 
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low temperatures. Sample A5 was cooled fairly slowly, over a period of 
many hours from 77°K to 10°K in the hope that some annealing would 
occur and thereby produce a less disordered solid, but the results, shown in 
fig. 3, are discouraging. 
3.3. Neon 

Neon solidifies at about 24°k, well below its Debye characteristic 
temperature of 0,~64°K. As a result we have fewer values for A (see 
fig. 4) in the temperature range where we expect that AocT’'. Because 


of the relatively low freezing temperature, the change in density during 
cooling from the solidification temperature to liquid helium temperatures. 
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should be less than in argon or krypton and a less disordered solid should 
be obtained. Certainly the increase in A with decreasing temperature is. 
more marked below 10°, although it falls far short of the sharp quasi- 
exponential increase observed by Wilks and others for solid helium. 

At 2°K, we find from our measurements that the mean free path of the 
lattice waves in neon is about 2 x 10-4 em. 


3.4. Krypton 
It appears from our results, which are shown in fig. 5, that above 


20 °K, Ax 7, but at temperatures below 10°K the conductivity is severely 
limited by imperfections of some type. Kr3 was annealed at 90°K and 
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Riek for some hours and at 55°K for about fifteen hours, but gave no higher 
conductivity than that shown by Krl which was cooled much more 
rapidly. Assuming \~ Ovl/3, then 1~ 10-5 em at 2:5°K. 
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§ 4. Discussion 


Apart from the general technical interest in the heat conductivity of a 
solid there appear to be some five fundamental interaction processes in 
dielectric solids which may be studied by means of thermal conductivity, 
namely the scattering of lattice waves by: 


(i) grain boundaries, which is usually important at very low tempera- 
tures (7' < @,,/30) ; 
(ii) dislocations, important also at very low temperatures ; 

(iii) point defects such as impurities or isotopes, a process which is 
particularly important at temperatures, 7'~ @,/20;. 

(iv) lattice waves in the temperature region just above the conductivity 
maximum where it often appears that A changes much more rapidly 
than with 1/7’; and 

(v) lattice waves in the higher temperature range where Aoc7'. 


It is apparent that in our experiments the method of making the samples 
produces severe disorder when they are cooled. This limits the con- 
ductivity markedly at the low temperatures where (i), (ii) (ili) and (iv) 
might be studied and unfortunately removes much interest in this tempera- 
ture range since we know little of the character of the disorder which is 
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introduced. Perhaps by growing inert gas solids under different con- 
ditionst, e.g. in a very thin plastic container which would contract or 
simply crack as the solid contracts, or by cooling under high pressure, this 
experimental limitation might possibly be overcome. Experimental diffi- 
culties would be very considerable, but it is of great fundamental interest 
to get further information about process (iv) in particular. Such informa- 
tion would supplement the available data concerning the validity of the 
theoretical expression, A~ f (7') exp (6/b7), in the temperature range below 
6/10; it might also show the effect of scattering by atoms of differing 
mass since krypton consists of a number of isotopes whereas argon and 
neon each consist principally of one isotope. 

We are thus forced to limit ourselves to considering the information 
obtained about process (v) from these experiments. That is, since the 
results seem to confirm that A\~ BT in this temperature region, we shall 
consider how experiment and theory agree in the evaluation of the constant 
B. Inthe table our experimental values of B are given (in mw cm~ from 
figs. 3, 4, 5) together with values of 6,, atomic weight A, melting tempera- 
ture 7',,, the lattice constant a, at 7’=0°K and approximate density 
pp at T=0°K; most of these values were obtained from the recent review 
by Dobbs and Jones (1957). Let us consider first the theoretical evaluation 
of the constant B. 


Physical Data for some Inert Gas Solids 


Te 4| On | oague| Wap | eee 
A (°K) | (RK) | (Ayo) (2 en?) ) (mw em-4 
Ne | 20-18 | 24-5 64 | 4:35 1-61 75 mn 
A 39-94 | 83:8 80 5:31 1-77 250 
Kr | 83:80 | 116 63 5-68 3-04 270 | 


Dugdale and MacDonald (1955) considered the properties of a linear 
chain with anharmonic coupling and proposed that the mean free path 
for lattice vibrations in a solid at temperatures 7’ $6, should be given by 


[~rolayT’ er Pierre Ven oct Ne, t18) 


where ry is the distance between nearest neighbours, « is the expansion 
voefficient, and y is the Gruneisen parameter. Thus they deduced that 


Ae VC [sayT Se a ee ane) 
or 


A=ryu/3yyT, since y=a/Cx. en aes ES 


_ + We are indebted to Professor G. O. Jones and Dr. R. H. Followell for 
information regarding their experiments on solid gases, and for suggestions as 
to crystal growing techniques. 
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Lawson (1957) has pointed out that eqn. (3) may be justified by the earlier 
work of Akhieser (1939) on the absorption of sound in solids; Akhieser 
showed that 


A=pv4/y*Twmax . . . ° ° . 5 . (4) 


where wmax=27vmax is the angular Debye frequency; Lawson then 
pointed out that since wmax ~ 7/19, y= a/Cy and v~ (py), eqn. (4) may 
be reduced to the form 

RESEND le ron, © oc etn rahe, backs. ald) 


in which the numerical constant is little different from that in eqn. (3). 
Another expression for A in terms of fundamental constants has been 
obtained by Leibfried and Schlomann (1954). They considered anhar- 
monicity in a linear chain generalized to three dimensions assuming central 
forces and only nearest neighbour interaction. With some approximations 


they deduced that 
12a 038 
ee Ee Bae a eS 
A rit yr Ae) 


where & and / are the Boltzmann and Planck constants respectively, NV is 
Avogardro’s number, A is the atomic weight and 5? is the atomic volume, 
i.e. for a face-centred cubic lattice with 4 atoms in the unit cell, a? = 483. 
Equation (6) may also be reduced to a similar form to eqn. (3) using the 
following substitutions : 


(k6/h) = ymax Lid v/[2ro 


GA) 2.25 
p=4M/a3 where atomic mass M = A/N 
and (esd Vay Qatar 
Thus from (6) we may obtain 
3 Ur 
\ eee a 
10 yy?! a 


where again the multiplicative constant is little different from that in (3). 
A further recent treatment of thermal conductivity in terms of anhar- 
monicity is due to Kontorova (1957) who shows that 


A=Cor,(f)?/3y2kT ee ee eee PS) 


where ¢ and 7 are the harmonic and anharmonic constants in the expres- 
sions for the force F or interaction potential U between neighbouring 
atoms : 

F=— he = ye" 

U = 472? — 472. 
Kontorova shows that by relating ¢ and 7 to y, « and 79, eqn. (8) becomes 
equivalent to the expression deduced by Dugdale and MacDonald : 


A=Coro/3ayT. 5. ope AM Mme rm (07) 
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It thus appears that experimental verification of any of these expressions: 
for constitutes an effective verification of them all. Because of our 
general lack of information concerning v, « and y in the inert solid gases, 
except in the case of solid argon, we choose to use eqn. (6) which may be 
written A=3-60,4 . a. 0,3/y27' w em deg-1; then using the values for 
A, a, 6, from the table we obtain A= 828/y?7' mw cm! deg! for Ne, 
3900/)27' mw cm-! deg for A and 4280/y?7' mw cm~! deg! for Kr. 

Concerning the value of y to be adopted, we have been led by the 
theoretical work of Barron (see review by Dobbs and Jones 1957) and by 
the experiments of Jones and his co-workers on argon to assume y~ 2-9. 
There is some evidence from this latter experimental work that for argon 
y may have a value as low as 2-0 at low temperatures but we are interested 
in the high temperature limit where both theory and experiment seem to 
lead to a value y~ 2-9. 

Thus we finally obtain from (6) that A~828/7?7’~98/T7 for neon, 
A~464/7 for argon and A~509/7 mw cm! deg! for krypton. These 
values are seen to be between 30 and 80% higher than the experimental 
values given in the table. 

In view of the uncertainties in the theoretical treatment and experimen- 
tal inaccuracies, the agreement seems tolerably good so that we might with 
some confidence apply a formula such as (2), (3) or (6) to calculate the 
approximate magnitude of the thermal conductivity in other cubic solids. 


§ 5. CONCLUSIONS 


These experiments have given values for the thermal conductivity of 
solid neon, argon and krypton that are severely limited at temperatures 
below about 10°K by physical disorder arising probably from severe 
strains imposed during cooling. However at temperatures approaching 
the melting temperature, the thermal conductivity appears to vary 
inversely as the temperature and has a magnitude somewhat smaller than,. 
but very comparable with, values predicted from theoretical considerations 
of the effect of anharmonicity in simple models of a lattice. 

It would be most interesting to obtain further results at low temperature 
on unstrained specimens of large grain size; in particular these might 
reveal details of the temperature dependence of the thermal conductivity 
in the range where anharmonic interactions of the Umklapp type are. 
becoming increasingly improbable; in the case of krypton the presence 


of more than one abundant isotopic species might also play a significant 
role. 
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Huceues anp Patevsky (1953) measured the total cross section o of cold 
neutrons in bismuth at different temperatures, and found that absorption 
plus the theoretical one-phonon scattering cross section o, were larger than 
the observed cross section at 300°K and 500°K. This discrepancy was 
attributed by them to the use of incoherent approximation in the calcula- 
tions. To test this point we have calculated, without making use of 
incoherent approximation, the one-phonon scattering cross section o, for 
cold neutrons at room temperature and find that the discrepancy is 
unresolved. 
Bismuth is a rhombohedral lattice with 


a=4:73564,. a=57° 14-2”. 
There are two atoms per unit cell, at wwu and aaa where uw=0-237. Using 


the Debye model of a solid and following Weinstock (1944), the one- 
phonon coherent scattering cross section is given by 


Sah ii | as q dq 
= M F?rexp(—2W,)A | ———+——___ 

* 16mck? J 2 P( ?) 1 exp (—fcq/kyT’)—17 
where S is the bound scattering cross section, W/, is the multiplicity of 
the reciprocal lattice vector 7. F* is the square of the geometrical 
structure factor, which for bismuth is 


F? = 4 cos? 27u(h+k+1). 
q is 27 times the phonon wave vector and has a maximum value q, 


given by q=k,.O,/hc, 0, being the Debye temperature and c the 
velocity ofsound. W, is the usual Debye Waller factor, given by 


672h2 (1 T\2 (a) 
i pati sel) 2 
a aor 4 3 (=) ,( 7 )he ; 


The constants used in the calculations are as follows: 


S = 9-37 barns B=35-17 x 10-*4 em? 
m= 209 a.m.u. T' = 300°K 
@, = 120°% o,=18m.b. per A 


c= 1-32 x 10° cm/sec 
Go= 1:19 x 108 em 
a ee eee 
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The limits A, and A, are calculated by the method described by Weinstock 
(1944) and integration over g is performed numerically. 

The sum of the calculated values of the one-phonon scattering cross 
section o, and the absorption cross section o, for different neutron 
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Total cross section (o,;+¢,) plotted against neutron wavelength. Curve A 
obtained by summing over the reciprocal lattice vectors. Curve B 
obtained by using the incoherent approximation. Circles indicate the 
experimental points. 


wavelengths is plotted in the figure, where for comparison, we also give the 
same obtained on the incoherent approximation. It will be seen that the 
one-phonon scattering cross section obtained by summing over the 
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various reciprocal vectors is larger than that calculated on the incoherent 
approximation. 

If we add to the one-phonon cross section the contribution coming 
from the higher phonon processes, which we estimate to be nearly 10% of 
the one-phonon cross section, the discrepancy between theory and experi- 
ment would further increase. It seems to us that the frequency spectrum 
of bismuth, on account of its peculiar structure, is very different from the 
Debye spectrum particularly in the high frequency region which contri- 
butes appreciably to the value of o,. This may well be a cause of the 
discrepancy. Unfortunately, no calculation of the vibrational spectrum 
of bismuth exists. 
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ABSTRACT 


Aluminium foils thinned by acid etching have been examined by trans- 
mission electron microscopy. Cellular patterns, 0-2 microns in spacing, 
have been observed superimposed upon, the normal subgrain structure of the 
metal. The cellular structure resembles patterns observed by other workers 
on replicas of chemically polished aluminium and interpreted as evidence of 
erystal misorientation. This communication shows that the cellular structure 
arises from a surface etching process and bears no relation to the distribution 
of dislocations in the metal. It is therefore most improbable that a further 
subdivision of the crystals is revealed. 


$1. INTRODUCTION 


HEIDENREICH (1949) and Hirsch et al. (1956) have shown that subgrain 
boundaries in thin aluminium foil can be revealed by transmission 
electron microscopy. Using the same technique, domains of similar 
character, 1-3 microns in diameter, have been observed by the present 
authors. However, superimposed upon this subgrain structure a fainter 
system of cellular markings, bounding regions only about 0-2 micron in 
spacing, has been detected. 

The secondary structure has not been reported in previous transmission 
studies, but in size and general appearance it closely resembles markings 
found on replicas of chemically polished (Hunter and Robinson 1953) and 
etched (Welsh 1955, 1956a) aluminium surfaces. Hunter and Robinson 
considered that the markings developed by chemical polishing represent 
the boundaries of unusually minute subgrain elements. This interpreta- 
tion has been questioned (Hirsch 1956, Perryman 1956, Welsh 1956 b) 
and doubt has remained whether the markings portray the internal 
structure of the metal or simply represent a special form of etching. Still 
smaller scale patterns (0-02—0-3 micron spacing), which frequently appear 
on anodically polished aluminium, are undoubtedly a surface phenomenon. 
associated with the electrode processes and not evidence of crystal 


imperfection (Welsh 1956 a). 
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The cellular markings described in this communication are shown to: 
be quite different in character to the subgrain structure and to bear no: 
obvious relationship to the arrangement of dislocations in the metal. 


§ 2. EXPERIMENTAL PROCEDURE 


Aluminium of 99% purity was supplied as beaten foil, 2-5 microns thick. 
The foil was annealed in vacuo at 350°C and then thinned to less than 
10004 by floating on dilute (1%) hydrofluoric acid. The etched foils 
were mounted on copper grids and examined in the experimental 
prototype Metropolitan-Vickers E.M.6 electron microscope. 


§ 3. STRUCTURE OF THE FoIL 


In fig. 1+ sharp polygonal boundaries delineate the substructure of the 
foil. That these truly represent boundaries between slightly misoriented 
blocks is confirmed by the fact that at higher magnifications individual 
dislocations are resolved, as at B in fig. 2. Moreover, as Hirsch ef al. 
(1956) first demonstrated, under stresses due to high electron beam 
currents the boundaries break up and the motion of dislocations can be 
observed. The misorientation across this type of boundary has been 
shown (Hirsch ef al. 1957) to agree approximately with values predicted 
from dislocation theory. 

Traversing the subgrains in fig. 1 a network of smaller cells is visible. 
In this region the cells are elongated and the major axes preserve a more: 
or less constant direction over large areas. In other regions the cells 
were equiaxed, as illustrated in fig. 3. By electron diffraction analysis the 
equiaxed cells were found to occur in regions where the surface of the 
foil had (100) orientation. The elongated cells appeared in regions of 
near (110) orientation, the major axis lying parallel to the [001] direction. 
The same orientation relationships were found in the structures described 
by Hunter and Robinson (1953). It will be noted that the larger polygons. 
in fig. 1 are equiaxed and do not change with the crystal orientation, 
and that the small cells pass without interruption over the subgrain 
boundaries (see also figs. 4 and 5). Since the latter undoubtedly comprise 
dislocation walls, the secondary structure apparently is quite unrelated 
to the main dislocation content of the metal. Moreover, individual 
dislocations visible within the subgrains were not preferentially sited at 
the cell boundaries and dislocations set in motion by thermal stresses 
passed through them without impedance. 

A further distinction can be drawn between the two types of structure. 
on the basis of contrast formation at the boundaries. Although 
dislocations could be resolved in the subgrain boundaries, the contrast 
often was due mainly to electron-optical fringes, as for instance at CC 
in fig. 2 ; this effect is well known and can be understood in terms of the 


yj All figures are plates. 
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dynamical theory of electron diffraction. These fringes did not occur at 
cell boundaries ; thus in fig. 5 the fringes at the subgrain boundaries are 
clearly visible, but the cell boundaries are revealed only by increased 
electron absorption. Concentric bands sometimes visible in the cells, as. 
at points A in fig. 1, are due simply to the variation in film thickness 
associated with the etching structure and are not to be confused with the 
fringes at misorientation boundaries. They are, in fact, contours of equal 
thickness similar to those observed by Heidenreich (1949) and are not 
confined to cell boundaries. 

To establish that the secondary structure does correspond to the 
markings described by Hunter and Robinson (1953), a specimen of the 
foil was thinned by chemical polishing in the reagent used by these 
authors (Alcoa Bright Dip R5—94 parts of orthophosphoric acid, 6 parts 
nitric acid). Foils treated in this way showed essentially the same 
features as those thinned by etching in hydrofluoric acid (fig. 4). It is 
important to note, however, that whereas Hunter and Robinson con- 
sidered the boundaries of the cells to be grooves, unequivocal evidence 
was obtained that they are actually ridges separating shallow pits. This 
conclusion is based mainly on the interpretation of shadowed carbon 
replicas of the foil surface, but is confirmed by the fact that in trans- 
mission micrographs the absorption of electrons was consistently greatest 
at the cell boundaries. Direct evidence of the differential surface attack 
was also available at punctured regions of the foil. Thus in some areas 
of fig. 6 (at D for example) the cell boundaries remain intact although 
the neighbouring foil has been completely etched away. Correlation of 
the outline of the sub-boundaries in fig. 5 with the etching structure also 
shows that the latter consists of shallow pits separated by ridges. It is 
significant that carbon replicas showed no change of profile at the 
boundaries of the larger polygons, emphasizing the difference in character 
between the two types of structure. 

Attention may be drawn to the minute dark spots of irregular shape 
which can be seen on many of the small cell boundaries in figs. 1 and 3. 
The spots appearing on films thinned by chemical polishing were more 
numerous, but much smaller than those found after hydrofluoric acid 
etching (cf. figs. 1 and 4) ; in neither instance did the spots appear upon 
the subgrain boundaries except where intersected by the boundaries of 
the cells. Carbon replicas showed the spots to be particles resting on the 
surface and not evidence of internal heterogeneity. Similar, though more 
uniformly distributed particles have been encountered while studying 
etched films of several other metals and apparently represent a deposit 
of the products of corrosion. They occurred even when using metals and 
reagents of the highest purity. 

It should be emphasized that the etching structures described in this 
section (except fig. 4) are typical of 2-5, aluminum foil of commercial 
purity further thinned by etching with dilute (1%) hydrofluoric acid. 
Spectrographically pure 0-5, foil thinned by similar etching did not 
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show. these structures though the sub-boundaries were clearly revealed 
(Hirsch et al. 1956, Bradley and Phillips 1957). Foils of other initial 
thickness and purity have not been examined. 


§ 4. Sepectep AREA DirrRACTION EXPERIMENTS 


The foregoing observations strongly indicate that the peripheries of 
the cells are not misorientation boundaries, but to substantiate this view, 
diffraction patterns were taken from selected areas of the foil. Mis- 
orientations within a selected area could be expected to cause splitting of 
the high index diffraction spots and splitting was in fact observed when 
the area contained a subgrain boundary. | However, when the area 
contained only cell boundaries even the highest index spots were invariably 
single, though often slightly distorted. This distortion was almost 
certainly simply the result of excessive film thickness at the ridges of the 
cells, a feature made manifest by a tendency to form Kikuchi lines in 
the patterns. 

Although the spot distortion makes it difficult to ascribe a definite 
limit to possible misorientation across the cell boundaries, any such 
misorientation must have been so slight that the dislocations accommo- 
dating it should have been resolvable in transmission micrographs. 


§ 5. Discussion 


Hunter and Robinson’s interpretation of the fine cell structure was 
first questioned by Hirsch (1956) who pointed out that if the cells truly 
represent misoriented crystal blocks their small spacing would imply 
a surprisingly high dislocation content for recrystallized aluminium. In 
the present study, direct observation of the internal structure of the 
metal has shown an absolute lack of correlation between the apparent 
dislocation distribution and the cell peripheries, the dislocations residing 
almost exclusively in polygonal boundaries of much larger spacing. The 
remote possibility that the cell boundaries did contain dislocations of 
such close stacking that they could not be resolved in the transmission 
micrographs has been negated by the selected area diffraction experiments. 
These show that any misorientation across the boundaries was so small 
that individual dislocations should have been easily resolved at the 
magnifications employed. } 

The question then arises, do the small cells in any way reflect the 
internal structure of the metal? The cell peripheries are now known to 
be ridges and not grooves ; there is indeed reason to believe that the 
cellular pattern develops by the merging of independent shallow etch pits, 
the cell boundaries representing the junction between pits (Welsh 1956 by: 
This interpretation has prompted the suggestion that the elongated 
character of the cells on surfaces of certain orientation derives from 
directional surface etching characteristics and does not portray the shape 
of any type of subgrain element (Welsh 1956b). The fact that the same 
orientation dependence is manifested by spurious patterns formed during 
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anodic polishing (Welsh 1956a), certainly lends force to this argument. 
However, the corrosion processes responsible for pattern formation 
during anodic polishing are definitely controlled by electrolytic factors, in 
particular the applied e.m.f. The structure discussed in this paper develops 
without electrolytic action other than any which may arise by local 
galvanic currents. In this respect the deposit of fine particles located on 
the boundaries of the cells (figs. 1, 3, 4, 5, 6) conceivably does indicate 
some type of galvanic process in which insoluble corrosion products 
deposit on cathodic sites, viz. the ridges. It may also be significant that 
Hunter and Robinson found the cell size to diminish as the purity of the 
metal diminished, an observation which could be interpreted as evidence 
that the number of anodic and cathodic sites increases with impurity 
content. 

The factors controlling etch processes on this scale are very imperfectly 
understood. Electron microscope replica studies have revealed domain- 
like etch patterns on several other metals (Delisle 1953, Barrett 1944), 
but the significance of these patterns is uncertain and any suggestion 
that they directly portray the shape of crystal elements must be viewed 
with caution. The present investigation certainly provides powerful 
evidence that the cellular markings found on aluminium are not the 
traces of misoriented subgrain blocks. 
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THE observation that fatigue failure in metals resulted from the spread 
of cracks from slip bands which had coarsened progressively during the 
application of fatigue stresses lead the present author to suggest (Feltham 
1957) that a close relation might exist between the fatigue limit o, in 
tension—compression and the critical tensile stress o, below which slip 
bands do not appear to form and coarsen in face-centred cubic metals in a 
tensile test. The object of the present work was to examine in how far 
existing experimental data and theoretical considerations provided 
evidence of the temperature-independence of the ratio o,/o,. which would 
be expected if o, were equal or proportional to o,. 

In face-centred cubic metals o, is generally regarded as the tensile 
stress at which Cottrell-Lomer barriers cease to be effective obstacles to 
the propagation of slip (Feltham and Meakin 1957). The temperature 
dependence of co, would be of approximately the same form whether the 
barriers became ineffective because (a) dislocations piled up against them 
could surmount them by cross-slip as suggested by Mott (1958) in his 
treatment of the spread of fatigue cracks, or (b) because the barriers failed 
by a break-through mechanism of the type considered by Stroh (1956); 
in either case a stress-induced dislocation reaction must take place. We 
shall consider both alternatives. 

Now, in relation to (a) the graphical representation of the temperature- 
dependence of the activation energy @ for cross-slip in copper (Schoeck 
and Seeger 1955) shows that 


Q~a+b log (a,/G), 


where @ and 6 are constants and G is the shear modulus. If the criterion 
for the occurrence of cross-slip is taken to be Q=mkT, with m a dimen- 
sionless constant, then the relation between log (c,/@) and 7 (and over 
limited temperature ranges between log c, and 7’) should be approximately 
linear. 

Figure 1 shows that the temperature-dependence of the fatigue limit 
of a number of metals} of purity approximately 99-99% does not deviate 
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t The fatigue limit has here been taken as the stress required to induce 
fracture after about one million cycles in tension—compression; the Woehler 
curves used were those of McCammon and Rosenberg (1957). The anomalous 
behaviour of aluminium may be due to oxidation. 
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greatly from linearity in the semi-logarithmic representation used. The 
temperature-dependence of o, therefore appears to have a functional 
form similar to that expected for o,. The variation of o, with temperature 
in the case of the break-through mechanisms (b) can be obtained from 


Fig. 1 


or > 


a xlO > dynjcnt 


T/T 


Temperature-dependence of the fatigue limit of metals of approximately 
99-99°/, purity, derived from the Woehler curves of McCammon and 
Rosenberg. > refers to Cottrell-Lomer barriers failing by recombi- 
nation. T',, is the melting point (°K). 


Stroh’s (1956) graphical representation of the ratio no,/2G as a function of 
T. Asn, the number of dislocations piled up against the Cottrell-Lomer 
barrier, is proportional to o,/G, the square root of that ratio, namely 


(no,/2G)2= x 
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is proportional to o,/G. The values of X for the case of the break down of 
the Cottrell-Lomer barrier in copper by ‘recombination’ (from Stroh’s 
fig. 8, curve ‘a’) is also shown in fig. 1; a similar curve (not drawn) was 
obtained for the case in which the barrier fails by ‘dissociation’ (Stroh’s 
fig. 4, curve ‘b’). 


Fig. 2 
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The relation between o;/G and = in copper for the case of failure of Cottrell— 
Lomer barriers by recombination (curve ‘a’, lower scale) and by 
dissociation (curve ‘ b ’, upper scale). 


The latter process requires a lower activation energy than break down 
by recombination. The values of © for both cases are plotted against 
o,/G in fig. 2; the temperature-dependence of G was assumed to be linear, 
with G = 5-18 x 10° kg/cm? at 77°K and 4-83 x 10° kg/em? at 293°K (Feltham 
and Meakin 1957). Line ‘b’, but not ‘a’, passes through the origin on 
extrapolation. 

The line drawn in fig. 3, correlating o, with o, (the latter variable taken 
from the measurements of Feltham and Meakin (1957) on polycrystalline 
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copper (99-99%)), has also been drawn to pass through the origin. Its 
equation is o,/o,= 2-0. 

The foregoing data and correlations show that the hypothesis of the 
temperature-independence of the ratio o/c, is highly plausible, but that 
further data of the type used in fig. 3 would be required to substantiate it 
fully. Two questions also require elucidation. First, the similarity of 
the temperature-dependence of the fatigue limit of face-centred cubic and 
of hexagonal metals (fig. 1), which suggests that in the hexagonal metals 
some form of thermally activated non-basal slip, i.e. cross-slip, may be the 
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Correlation of experimental values of of and o, for polycrystalline copper 
(99-997): 


pre-requisite for the formation of fatigue cracks, and, second, the 
significance of the unexpectedly high value of 2-0 of the ratio o,/o, in 
polycrystalline copper. 
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ABSTRACT 


The paper describes measurements of the ratio of z-mesons to electrons as 
a function of shower size and barometric pressure for showers containing about 
5x 10° particles at sea level. Measurements on the change of rate with 
atmospheric pressure and temperature are also described, and it is shown how, 
on the assumption that all showers suffer the same attenuation, these measure- 
ments yield the attenuation lengths of electrons and p-mesons. At each stage, 
the measurements are compared with the results of other experiments, and 
finally with various theories of the nucleon cascade. It is shown that none 
of the theories fits the data well, and that one may suspect the assumption of 
similar attenuation for all showers. Possible mechanisms giving rise to a 
distribution of attenuation lengths are briefly discussed, and a case made for 
an experimental investigation of this point. 


INTRODUCTION 


OnE of the main objects of research on extensive air showers is to achieve 
a satisfactory theory of the nucleon cascade. This is important for the 
bearing of air shower results on astrophysical problems and the problem 
of the origin of cosmic rays, because without it, we have only a rough 
knowledge of the energy and nature of the primary particles. But it is of 
out-standing importance for the light it may shed on the nature of nucleon 
collisions at energies far in excess of the energies available from machines. 

The early hypothesis that air showers were electromagnetic cascades 
from electron or y-ray primaries was discarded largely on the strength of 
arguments relating to the shower absorption length, A,= — (dln k/dt)™, 
where ¢ is the depth in the atmosphere, and F the rate of observing showers 
of given size. (The arguments are given in detail in Dobrotin et al. 1956, 
and Greisen 1956). A, can be determined from the measurement of 
altitude variation, of zenith angle distribution, or from the barometric 
effect. In general, each method requires that small corrections be made 
to the directly measured quantities before A, is deduced, and then by 
combining this with knowledge of the spectrum, A,, the attenuation lengtht 
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t Throughout this paper, we shall use the term * absorption length’ for 
quantities of the type d¢/d In R, referring to changes in the rate of counting 
showers with depth, and ‘attenuation length’ for quantities of the type 
t/ In N referring to the rate of change of number of particles with depth. 
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of particles in the showers can be deduced, and compared with the attenu- 
ation length of the E.M. cascade. 

With the introduction of the nucleon cascade hypothesis, the number 
of variables in the theory is greatly increased. For any choice of para- 
meters in the nucleon collision, the resulting cascade can be worked out, in 
some cases analytically, and in others by numerical methods, and the 
results compared with experiment. In order to make a choice between 
different models, it is important to determine as large a number as possible 
of experimental quantities which can be compared with the theory. 

In this paper, we describe some measurements of attenuation lengths of 
electrons and p-mesons, and of the ratio of ~-mesons to electrons as a 
function of shower size. The results are compared with the results of 
other experiments, and also with the predictions of various theories. 


A. MEASUREMENTS ON THE ELECTRON COMPONENT 
§ 1. APPARATUS 


The apparatus is the array of 91 G.M. units described in I (Cranshaw 
and Galbraith 1957 a). 


§ 2. THE MEASUREMENTS OF ),, THE ELECTRON ATTENUATION LENGTH 
2.1. Theory of Barometer Coefficient and Temperature Coefficient 
It has been found in many experiments that the rate of recording air 
showers varies exponentially with depth, i.e. 
Olnk 1 
ee Se el ee 
a OA ) 
It is also found that the size spectrum of showers can be represented 
approximately by a power law, i.e. 


J(N\GN =hN-Vol Gg Ni ON) ees eee en 
where JN is the number of particles in the shower. Hence we can write 
J(N)dN =k, exp (—t/A)N-7" dN. iJ ae eee eS) 


This equation can be explained if we suppose that the particles in all the 
showers are being attenuated exponentially, with an attenuation length X,, 


ein Nie ey sa 
oT deny 5A ea 
Thus A, can be determined by measuring y and A. We may remark that 
eqn. (4) merely defines a quantity A,, if there is a mixture of showers present. 
In Paper V, we have shown that y=1-88 for the showers with which 
we are concerned. We now show how A can be determined from the 
barometer and temperature coefficients. 
The barometer coefficient is defined as 


p=dinkidp) 2 fee ene) 


where in our case p is the pressure in cm Hg. The temperature coefficient 
is defined by 


6=dln R/oT oy ce oe ee) 
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where 7’ is the temperature in °c. In general, 8 is made up of two effects, 
the true mass absorption, Bmass, of the showers by the matter above the 
apparatus, and a variation of the probability of counting showers due to 
the change of density, and hence of the lateral spread of the showers. 
This density effect is most easily observed as a temperature coefficient. 

Hodson (1951) has calculated the effect for an array of counters which is 
small compared with the spread of a shower, and finds 

oink 2y—2-a, 
= a | ost ire ai 
oT 2:8 /o 7) 
where oo is the exponent in the relation Rad~*, giving the counting rate as 
a function of counter separation. In a particular case, Hodson has found 
6= —0-38+0-11% °c in good agreement with the theoretical value 
— 0-25. 

We can calculate the effect for a lattice array of counters, by assuming 
that the density distribution of particles in the shower is given by 


ip pNeE 
= eae oe bf ee ee ee LS 
(=) f: & (8) 


where ry is the radiation length, and a function of density, and the expo- 
nent a varies sufficiently slowly with r that for the present purpose it may 
be considered as constant. Then for a change in 7», the mean number of 
particles passing through a counter will change by 797? 


] 


olnk 
1.€ = —2 eon Hei ke Bee 
1.6. nia (a—2) (9) 
Oink y(a—2) 10 
‘a Oars eT ee 


when the change is due to a change of temperature at constant pressure. 
Thus 
_oeamRk- y(a—2) 
AE) aD 
We note that in the case of small apparatus, @ is always negative, but for 
a lattice arrangement, @ will always be small, and may be positive or 
negative, because a@ is near 2. 
The correction to 8 due to the density effect can now be deduced by 
using the gas laws. We have 


UE Oeil, see Rote tee e 


Oe eee eeentat a4 aie s0 (12) 


"9 we p 
where 7’ is measured at constant p, and p at constant 7. Then 


Bmass = B— 40 . : ° . ° . . . (13) 


_olnk topo lncinw 5 ol (8-40). ; (14) 
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Finally, a small correction can be made for the fact that not all the 
observed showers are incident vertically. In fact, if J( > NV) is the isotropic 
flux far from the earth of particles which would produce showers con- 
taining >N particles at sea level at vertical incidence, A is the rate 
absorption length for vertically incident showers and ¢t the height of the 
atmosphere, then the rate at sea level will be given by 


74/2 
R(>N)=J(>N) | In sin 0 cos Oexp [t/A(1 —sec 6) ]d8. (15) 
0 


The integral may be regarded as the solid angle subtended by the 
apparatus, and putting in the values of ¢ and A, we find Q=0-5 steradians. 
The value of dln R/dt is given by 
cues oh 1-3 yk ee anos 
ot A E,(t/A) 
where H,,(x) is the exponential integral. The expression in the bracket has 
the value 1-08. We therefore write finally 


II olnk OolnkR 
— = = Oy = —0:068 = — (7-:068(8 — 46 ‘i 17). 
A aI a (8 ) (17) 


where dln R/dp refers to constant air density. 


2.2. Huperimental Method 


In each hour interval, the apparatus (Cranshaw and Galbraith 1957 a) 
records the number of showers which have triggered >3, >5, >8 units on 
the array. In addition, the barometric pressure and the temperature 
inside a Stevenson screen on the roof of the building are also recorded. 
The data have been collected over a period of about three years, and a 
bivariate regression analysis has been used to extract the barometer and 
temperature coefficients. 


2.3. Results 


The results are shown in three rows of table 1, where the third 
column gives the barometer coefficient and the fourth the temperature 
coefficient. It is important to note that a solar diurnal time variation 
(Cranshaw and Galbraith 1957b) could give rise to an apparent tem- 
perature coefficient due to the diurnal temperature wave (Farley and 
Storey 1957). If there are sufficient counts available this can be avoided 
by averaging over 24 hour periods, or by using the change of rate with 
season (Hodson 1951). In the present case the statistical accuracy on 
the temperature coefficient is poor. We can calculate the expected tem- 
perature coefficient by noting that the distances at which the probability 
of a unit being discharged becomes one half for the three shower sizes 
chosen are 80, 110, and 130m, and the value of the exponent in the power 
law approximation to the density distribution at these distances is 2-18, 
2-28 and 2-34. Using eqn. (11), we find that the expected values for the 
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temperature coefficient are 0-12, 0-20 and 0-23% °c-!. The experimental 
values are clearly consistent with these values, and we have used the cal- 
culated values to obtain the mass absorption effect in column 5 from eqn. 
(13). 

In the sixth column we have calculated A using eqn. (17), and in the 
seventh, we tabulate the value of y from y=1-71+0-13log N/10°, which 
seems to fit the experimental data satisfactorily. Finally in the eighth 
column, we calculate A,. The figures are consistent with a constant value 
of A,, but a small increase of the order of 10% over the range of shower 
sizes considered cannot be excluded. 

At the end of table 1, we have included three measurements by the 
M.I.T. group, who have determined A by studying the dependence of the 
shower rate on zenith angle. The values are in good agreement with the 
values deduced from measurements of barometer coefficients. 


Table 1 
Barometer | Temperature " Mass : 
Experiment N coeff. coeff. BESOED EON A is y A 
6) emai OF Sq-1 effect gm em gm cm-* 
70 70 % em-1 
Citron (1958) | 2 x 108 — —0:12+0-03 | 10:6+0-4 138 1-36 188+6 
8 x 108 — —0-28+0-02 | 11:9+0-3 123 1-44 178+4 
104 _— —0-34+0-03 | 1214+0-4 121 1-45 Via 5 
2c LO = —0-37+0-04 | 12-9+0-6 114 1-49 170+7 
Miller (1951) |} 5x 104] 11-1+1-0 ~ —0:3 12*3 41-5 120 1-54 185+ 25 
Hodson 
(1951 — —_— —0-:38+0-11 — — — — 
(1953) HSeLOz 9-0+1-1 — 10-5+1-5 140 1-54 216+30 
Farley (1957) Oe 9-9+0-7 —0:89+0-24 | 13-0+1-5 114 1-45 164+ 23 
3x104 | 10-5+0-8 —0:37+0:26 | 12:0+1-7 122) 1-51 185 + 25 
9x104| 14:2+1-8 —0-37+0-60 | 15:6+2-5 94 1-57 148+ 30 
3x10°| 13°6+3-7 —2-141:3 14:-64+4-5 100 1-65 168+ 50 
Harwell 
(1954) tf 3x10°| 13:0+0-6 ~0 13:0+1 114 1-65 187+15 
8x 10°} 13:5+1-0 ~0 13-5+1-5 110 1-70 185 + 20 
2X 108) 13:34-2-0 ~0 13°3 + 2:5 111 1:75 195+ 30 
Culham 5x10*| 13:6+0-°6 0-08 + 0-04 | 13-3+0-7 110 1-80 199+11 
107] 14:5+1-3 0-24+0-17 | 13-8+1°5 106 1-84 195 + 20 
2x107| 14:2+2:-6 0-54+0-35 | 13:-4+2-8 106 1-88 200+ 34 
Clark (1957) 10° — = == LOT 11 = a 
Clark et al. 
(1957) 10° == oie = 98+5 = Ss 
10° a = — 90+5 — == 


+ Unpublished results obtained with apparatus described in Cranshaw and Galbraith (1954). 


§ 3. CONCLUSIONS 

The value of A, the rate absorption length of showers can be determined 
by the measurement of barometric and temperature coefficients. The 
value of X,, defined by A,=yA does not increase by more than about 10% 
when the shower size is increased from 2 x 10? particles to 10’ particles, 
but is nearly constant and equal to 190g cm. 

The solid angle subtended by typical shower apparatus is about 0°5 
steradians. 
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B. MBASUREMENTS ON THE PENETRATING COMPONENT 


§ 1. APPARATUS. 


The apparatus is the array of 91G.M. units described in I and the four 
penetrating particle sets described in IV (Porter et al. 1957). For the 
earlier part of the experiments, the p.p. sets were in the positions marked in 
fig. 1 of IV, but later all four were moved to position 3. Data from both 
dispositions were used in §3, using a suitable normalizing constant to 
combine the results, but only position 3 results in § 2. 

As in the earlier paper, either or both of the two interlaced trays in each 
unit may be triggered. We refer to events in which one bank only has 
been triggered as a 1 event, and one in which both banks have been triggered 
asa 2event. Thus, for example, an event 2+1 indicates that in one set 
both banks were triggered and one bank in another. In any particular 
case, we are unable to say whether a 2-event is caused by a nucleon or by 
a p-meson pair. However, we show that a statistical separation can be 
made. 


§ 2. Taz RATIO OF u-MESONS TO ELECTRONS AS A FUNCTION 
OF SHOWER SIZE 


2.1. Experimental Method 


The total number of w-mesons in a shower at sea level is rather difficult 
to determine, because over a large range of distances, the density falls off 
roughly inversely as the distance from the core, and thus a large fraction of 
the total number lies at a great distance, where it is not easily measured. 
Accordingly, we have decided to measure the density at one distance, 200 m, 
as accurately as possible, for two shower sizes. The distance 200m was 
chosen because it is far enough from the shower axis to be undisturbed by 
nucleons, and small enough for showers at that distance from the centre 
of the array to be substantially unaffected by edge effects. The two 
shower sizes chosen were those showers triggering just three units on the 
array, and those triggering between 8 and 12 units. 


2.2. Results 


The showers were divided into groups according to the distance of the 
centroid of the triggered units from the j.-meson detectors and the response 
of the u-meson detectors. The results on 2029 showers of the first size 
group and 258 of the second are shown in table 2 (a) and (b). In the first 
column are listed the possible responses from the ,.-meson detectors, and 
in the subsequent columns the number of times that response was obtained 
from a shower whose axis lay at a distance indicated at the top of the 
column. 

By inspection, it is clear that there is no evidence for any nucleons at a 
greater distance than 120 metres, i.e. the number of events of type 2 is not 
larger than one would expect from the number of events of type 1+1, 
taking into account their relative statistical weights. The two events 
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2+1+1 at 200-240 metres appear to be exceptional and these will be 
referred to later. The «-meson densities for 200m have been calculated 
from the data in the neighbourhood of 200m by a least squares method. 
The values found are 0-38 40-03 and 1:35+40-12 particles em-2. These 
values, and the slope of the density distribution are in good agreement 
with the values obtained in IV. 


Table 2 (a). 3-fold events 


Distance, metres 0-40 40-80 80-120 120-160 160-200 200-240 


0 10 97 163 278 370 619 


l 5) 61 60 88 68 94 


ee es es 0 0 1 0 0 


2+1+1 1 4 1 i 0 2 


242 0 0 0 0 0 


ae es | 


The mean numbers of electrons in the groups of showers have been 
estimated by using the data of V, (Cranshaw ef al. 1958). 


iN N-yP(N, f)dN 
N aa 

| N-y1P(N, f) dN 
J 0 


and we find V,=3-4 x 10° and N5_-12=1-9 x 10’, in satisfactory agreement 
with III and V. Accordingly, if we write 


N,=kN 
we find «=0-75 + 0:18. 
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2.3. Comparison with Other Experiments 

Other values which have been found for « are 0-87 at mountain altitudes 
(Cocconi et al. 1949), 0-89 at sea level (Milone 1952) and 0:75 + 0-09 
(Kraushaar 1957). (Data given by Dovzhenko et al. 1956 allow one to 
determine a value for « of about 0-60. However the experiments were 
not made with this measurement in mind, and it is not clear how much 
systematic errors may have contributed to this low value.) The values 
might be taken to suggest that « decreases slightly with shower size, but 
all results are consistent with the average value 0-81. It is interesting to 
note that if this value is correct, then the ratio of j-mesons to electrons 
at 200m from the core changes from 13° for showers containing 3 x 106 
electrons to 40% for showers containing 10 electrons. Thus beyond 200m, 
we must expect the majority of particles in showers containing 104 
electrons to be «-mesons and their secondary electrons. 


2.4. Hxceptional Events 
The two events referred to in § 2.2 are shown in fig. 1 A and B, where a 
solid circle indicates that the 20 em? ‘b’ counter has been triggered, and 
an open circle indicates that the two 200cm2 ‘a’ counters have been 
triggered. (A dot and circle indicates that both ‘a’ and ‘b’ counters 
have been triggered.) The position of the p.p. sets is indicated by a 


square. 
Fig. 1 
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Both events are slightly unusual even without the p.p. set results, 
because the triggered ‘a’ units are considerably spread out, over a distance 
of 300 metres in A, and 230 metresin B. In A, one might suppose that the 
axis has fallen somewhere between the p.p. sets and the two nearest trig- 
gered units, but then the fact that three counters, one of area only 20 cm? 
at a distance of 300 metres are triggered, suggests a local condensation of 
particles in that area. In B, two ‘b’ counters have been triggered in a 
region where the density is expected to be ~2 particles m~, i.e. the pro- 
bability of triggering such a counter is about 0-5%, In the 2000 showers 
examined, this is, of course, not highly improbable, but the two showers 
considered here were picked out for their peculiarity with respect to the 
p-p. sets. Taken together, the evidence suggests that a small fraction of 
the showers may have an ususual structure, as has been suggested before 
(Hidus et al. 1952, Cranshaw 1955). 

3K2 
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It is perhaps suggestive that the showers appear to be elongated in the 
North-South direction, as they would be if they had arrived at a large 
zenith angle. The counters are cylinders with their axes lying East-West, 
so that the projected area is unchanged for showers inclined in the North— 
South plane. We might suppose that because of the large zenith angle, 
or for other reasons, for example a heavy primary nucleus, the showers 
may consist almost entirely of .-mesons, the density of which falls off only 
slowly with distance from the axis. The observed patterns would then 
be easily understood. 


§ 3. THE NuMBER OF NUCLEONS 


The data of table 2 may be used to deduce the number of nucleons in the 
showers, using the methods of IV §4.2, and V §4. By considering the 
relative numbers of events such as 1+1+1 and 2+1, we see that almost 
all the events with double bank coincidences at distances less than 120 
metres are due to nucleons. A small correction for the number due to 
pairs of »-meson can be made. Then using the formula of V, we find that 
the numbers of nucleons detected in the two categories are 9-3 x 108, and 
8:5 x 104. As described in IV, the probability of detecting a nucleon in 
the detectors is 0-77, and the energy required is about 3x 10%ev. We 
conclude that showers containing 3-4 x 10® and 1-9 x 107 electrons contain 
respectively 1-2 x 104 and 1-1 x 10° nuclear interacting particles of energy 
greater than 3 x 10% ev. 

Alternatively, we may take the data without reference to distance of the 
axis from the detectors, and apply the same methods, correcting the 
“nucleon events’ like 2+1 by the number of u-meson events 1+1+1, 
with a factor for their relative statistical weights. In this way, the 
figures are 1-2 x 104 for showers containing 3-4 x 10° electrons and 2-5 x 10° 
for showers containing 2-2 x 107 electrons. These figures are in reasonable 
agreement with those given in the preceding paragraph. The error is 
difficult to estimate, and depends on the structure of the nucleonic com- 
ponent of the showers. 


§ 4. THe ATTENUATION OF u-MESONS 
4.1. Hxperimental Method . 

For this measurement, we have counted the number of showers observed 
for one ».-meson detected as a function of pressure. The analysis simply 
consists in sorting the cards recording the showers into categories of 
different pressures, and then finding the fraction in which ,.-mesons are 
observed. When p-mesons were observed in two or more units, eon er 
events, these were counted twice. However as has been shown, an appre- 
ciable number of 2 events are due to nucleons, and so such events were 
excluded from the analysis. 


4.2. Results 
We have used a least squares analysis to find the rate of change of pu :e 
ratio with pressure. The value found is 4:5+2°%,em- with the ratio 
increasing as the pressure increases. This is the slope for constant J,. 


82} 


Observations on Extensive Air Showers: VI 


‘[eAjNou oq OF powmMsse ov suOSOUT-4 fo PATYyy-9UQ *4/,(q4)o=wu Aq UeAIs st Uotonpoad uosout-x yo Aprordyjnur oy, t 
‘Sou0d PLBMYOVG PUL PIVALOF UT poyesyuoUOD ADIOU “JT “([EG{) ‘WIeYyDg-1eqey] 07 Sutpaioooe Adosyostuy *T 4 


= e-OLX@- ¢80-0 18-0 O00FI~ OST OST | [equcuttiedx ay 
zs 0-1 0-o-G@1 z-O1 — = a = CST — | (Zg6T) [equezoxy 
AIOIY} 8,5.10qUISTOF{ 
09 oy suryeurxoidde c-OLX9 60-1 £9-0 0S9 009 OLT OSI | (SG6T) 2401107 
[Pepoul [VOVUTEyye IL pure turpng 
$1 = = II 24-0 
09 O-T a a = a 
€1 = — 4109-0 (LG61) ®PO 
06 €-0 Z, s-O1 X &-h Z10-0 GLO |201X¥-9— eOLXL-1-— | — OLZ . 
06 G-0) z e-OLX LT LGO-0 OL:0 | CL XF3— Oe | OLT | (LS61) 83150 
06 OT Ye c-O1 XG 90-0 19-0 Ole cOl ~e— | = OSI pue vpoy) 
= — a 901 901 — 901 »01 901 rOL sopoqaed 
OZIS TAaMOYG 
ae bs tp "NIN ° NIN ” g-0 3 "Y g—Uld 3 °y Aqryuenty 


id MICA 


$22 T. E. Cranshaw et al. on 


If we want the slope in an individual shower, or the average shower, 
we note that an increase of pressure of 1 cm decreases the number of 
electrons in the showers by 7-4°% and therefore decreases the p :¢ ratio by 
7x 0-:19=1-4%, i.e. in an average shower, the  :e ratio is decreasing with 
height by 5-9°/, em Hg. 

4.3. Comparison with Other Eaperiments 


Both Cocconi (1949), Milone (1952), agree that the fraction of ~-mesons 
in showers of density ~ 60 particles m~? at sea level and mountain altitudes 
is 2-65% and 1-46% respectively when the showers are detected solely by 
counters distributed over asmall area. According to Greisen, the showers 
which trigger a close array of counters at mountain altitudes would, at sea 
level contain 0-29 times the number of particles in the showers triggering 
the same apparatus at sea level. Therefore for these latter showers, the 
fraction at mountain altitude is 2-65 x 0-29-919=3-24%. Moreover, the 
mean distance from the axis at which showers are detected is inversely 
proportional to the densities about two radiation lengths above the 
observation point, i.e. the distance at sea level is (644)/(954)°°? ~ 0-73 of 
the distance at mountain altitude. Since the p :e ratio is almost directly 
proportional to the distance from the axis, the p :e ratio at a comparable 
part of the shower is found to be 4.5%. The decrease in ratio with altitude 
in one shower is thus 4:5/1-48=3-0 times. This corresponds to change of 
pe ratio of 5-0% cm Hg, with an error which may be about 2%. 

By combining a number of observations at mountain altitude and sea 
level, Greisen (1956) finds that the total number of x.-mesons in showers at 
sea level in twice the number at mountain altitudes, in agreement with 
Dobrotin et al. (1956). This corresponds to a change of » :e ratio in the 
average shower of 4:5°%, with an error again of about 2%. 

As in the determination of A, one would expect to obtain the same co- 
efficient for the change of p :e ratio by measuring the pz :e ratio at different 
zenith angles. Measurements of this type have been reported by the 
M.1.T. group (Kraushaar 1957) and they find that yw :e ratio is independent 
of zenith angle, in apparent disagreement with the above results. These 
experiments are referred to again in the discussion. 


§ 5. CONCLUSIONS 


Except in the calculation based on the total numbers of ,.-mesons at 
mountain altitudes and at sea level, which may not be quantitatively very 
accurate, we have neglected the fact that the ~-mesons are probably 
expanding from the shower axis as the shower proceeds down through the 
atmosphere. Then measurements which are made at a fixed distance 
from the axis, or which integrate the number of j-mesons up to a fixed 
distance will show an apparent absorption, due to the fact that particles 
are diffusing out of the measured region. In the absence of other infor- 
mation, we have to make assumptions if we wish to estimate the magnitude 
of this effect. If all the ,.-mesons were produced at an infinite height, and 
diffused from the shower core by scattering, it can be shown that they 
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would be expanding as though they diverged from a point twice the scale 
height of the atmosphere, i.e. 16km, above the observation point. A 
minimum value for an apparent absorption due to this geometrical effect 
is then 0-65% cm, but if the height of origin of the x-mesons is quite low, 
and the production angles play a significant part in producing divergence of 
the u-mesons from the shower core, then this value could be about 2% cm. 

To sum up, we have found from three experiments a decrease of the p :¢ 
ratio of 5+ 2% cm, and it would be unwise to suppose that the final result 
is known to be much higher accuracy than this, because the same calcu- 
lations subject to the same errors appear in each result. Since the electron 
component is being attenuated by 7-4°% cm, we find an apparent j.-meson 
attenuation of 2-4°% cm. About 1-4% may be allowed for the geometrical 
effect mentioned above. We are thus left with a true .-meson attenu- 
ation of 1+3%. We may express this result as 1/Au=7-3+15 x 104g—1 
em*, or the observed value of A, is 1400 gem, with a 60% probability of 
lying between 340 and — 680 gem-?. 


§ 6. Discussion 


In table 3, we have collected the results of theoretical calculations of 
the nucleon cascade, and also the experimental data as discussed pre- 
viously. We wish now to see how far the experiments enable us to make 
a choice between possible models, and how far the theoretical work en- 
ables us to design future experiments that will be significant. 

First we consider the values of A,. Ueda and Ogita (1957) have 
remarked on the great difficulty of finding a theory which gives a constant 
or nearly constant value of 4, over four orders of magnitude of V,. They 
suggest that the inelasticity, 7, in nuclear collisions may change from 0-3 
at low energies to near 1 at high energies. Apart from the fact that it 
seems unlikely that a large enough range of shower sizes could be covered 
in this way, it seems slightly artificial to suppose that 7 changes in just such 
away astomake,,constant. Wesuggest therefore that another mechanism 
which operates to equalize the value of A, should be examined. This 
mechanism is the possible existence of a mixture of showers each being 
absorbed with a different value of X,. Then the showers for which A, is 
large survive to a lower level, and tend to raise the value of yA which is 
measured there. On this picture, the value of A,=yA does not necessarily 
refer to the attenuation of any particular shower, and the apparent con- 
stancy of A, with altitude does not mean that any particular shower is 
attenuated with a constant length. Such behaviour could seriously 
falsify attempts to estimate the energy of showers by evaluating a track 
length integral, in the manner of Greisen (1956). ae 

The mixture could arise in several ways. We would here mention first 
the possibility of a mixture of primaries, so that showers of the same energy 
may result from interactions of widely different energy per nucleon. 
Secondly, if the primaries have an interaction length comparable with 
A~ 100g cm~, showers of a given size can result from primary interactions 
distributed through quite a considerable depth of the atmosphere, and with 
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therefore, a wide distribution of age. The fact that in the M.I.T. experi- 
ment (Clark et al. 1957) there is no distribution of age apparent from the 
structure may mean either that the changes are too small, or that the 7° 
particles which contribute most to the shower at sea level are produced not 
in the first interaction, but much later. An extreme result of this effect 
would arise if the interaction length of the primary were in fact 100g cm~. 
Then A would always be measured as 100g em~*, and the average charac- 
teristics of showers measured at all altitudes would be the same. A third 
possibility is that there may be large fluctuations in the development of 
the cascade. If, forinstance, 10°/ of high energy nucleon—nucleon collisions 
have 7=1, instead of 0-3, then we may expect one such collision to occur 
in each shower. But the development of the cascade will be very different 
if this collision occurs at the beginning, or near sea level. 

Next we consider the value of A,. Here the experimental uncertainties 
are rather large, and in general, the measured values could be in reasonable 
agreement with the theories. The result of the M.I.T. measurement 
(Kraushaar 1957) however is anomalous. They find that the attenuation 
of -mesons measured by studying showers at different zenith angles is the 
same as that of the electrons. This is, of course, just the result which 
would be obtained on the extreme hypothesis mentioned above, but then 
the measurements of long attenuations for the .-mesons would be hard to 
explain. It is also true that when a shower is produced in the atmosphere 
at a large zenith angle it spends more of its development in air of low density, 
so that the probability of a »-meson decaying becomes appreciable at higher 
energies. Thus at sea level, we may expect fewer »-mesons each with 
higher energy in showers at large zenith angles. However, approximate 
calculations show that unless appreciable w-meson production occurs by 
other mechanisms, this effect would only account for an apparent attenua- 
tion of about 2% em. One may expect that the further investigation of 
variations with zenith angle would be profitable. 

The value of « is as yet, not well determined. A lower value than 0-81 
is certainly possible, as would be deduced from the figures published by 
Dovzhenko (1956). However, in general, the measured values seem higher 
than calculated. This might indicate that the energy per nucleon is not 
increasing with shower size as rapidly as is assumed in the theory, as would 
be predicted by the hypothesis of a mixture of showers. Alternatively, 
it might indicate that there is a progressive shift in the spectrum of primaries 
towards heavier nuclei, as would indeed occur if the energy spectra of the 
individual nuclear species were approximately power laws with an ex- 
ponent increasing with energy. Then as the energy increases, heavy 
nuclei are relatively favoured. On the other hand, the fact that « is less 
than one is an argument against the hypothesis of an energy cut-off and 
the production of larger showers by larger primaries (Peters 1957) or dust 
grains (Alfven 1954), because such a hypothesis would indicate a constant 
energy per nucleon above some critical shower size. Further, in such a 
theory, larger showers would be initiated on the average higher in the 
atmosphere than small showers, leading to a value of « greater than unity. 
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The calculation of Oda (1957) shows that the effect of introducing an 
anisotropy in the nucleon-nucleon collisions concentrating the energy 
in the forward and backward directions, which is required by considerations 
of angular divergence, makes the calculated value of « further from the 
experimental value. 

The ratio of «-mesons to electrons is probably not a very significant 
quantity in testing the theories because it depends rather strongly on the 
lower energy limit. However Rozental (1952) shows that the observed 
values are inconsistent with multiplicities in the nucleon-nucleon collision 
which increase more rapidly than #, with v=, and Ueda and Ogita have 
shown that the height of origin of the «-mesons is sensitive to the value of 
7. More careful measurements of ;.-mesons may be able to fix v and 
more precisely. : 

We conclude that at the present time, it is important to perform experi- 
ments designed to detect the presence of a mixture of showers in the flux 
at sea level. Experiments on the j:e ratio in individual showers are in 
preparation. 
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ABSTRACT 


The energy density of the electron—photon component has been measured 
at distances of 60-500 metres from the axis of showers containing 3 x 10° 
particles at sea level. The mean electron energy appears to be constant, 
and close to the critical value of 84 Mev, over this region. A lower limit for 
the energy appearing in the soft component of these showers at sea level 
is set at 6 x 1014 ev. 


§ 1. INTRODUCTION 


THE estimation of primary energy in extensive air showers is usually 
based on the total number of particles observed at a given altitude, with a 
calculation from some model of the cascade in the atmosphere. The 
energy derived is rather insensitive to the form of the theoretical model 
used, provided that fluctuations in the shower development are taken into 
account, but depends on the assumption that the particles observed are 
true cascade particles, with energies greater than about 84 Mev, the 
critical energy in air. Since a Geiger counter array will detect electrons 
with energies greater than 5 Mev, a verification of this assumption is 
desirable. 

In the experiment described here, the distribution of energy in the 
electron—photon component has been measured by absorption in a deep 
Cerenkov counter, at distances of 60-500 metres from the axes of showers 
containing about 3 x 10° particles. This region contains about two-thirds 
of the total number of particles in the showers. Previous measurements 
have shown mean electron energies close to the critical energy, but refer 
to showers of undefined size, and averaged over unknown distances from 
the axis (Mitra and Rosser 1949, McCusker and Nevin 1951, Rosser 
and Swift 1951). Recently, Matano et al. (1957) have found energies 
considerably larger than the critical energy close to the shower core, « 
but falling off below the critical energy at distances of about 50 metres. 
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§ 2. APPARATUS 


A diagram of the Cerenkov counter is shown in fig. 1. It was square in 
horizontal section, with a sensitive area of 14 400 cm”, 92 cm deep, and 
filled with de-mineralized water. The water was contained in a sealed box 
made of the commercial material ‘ Darvic’, which has a white diffusive 
surface. A single 5-inch photomultiplier was mounted at the top in 
optical contact with the water. The lead shield over the counter was 
0-6 em thick, and acted as a photon-converter, as well as improving the 
response slightly at high electron energies. The counter was calibrated 
continuously with «-mesons, which were selected by a small Geiger 
counter telescope. A single relativistic particle traversing the system 
vertically produced 7-10 photoelectrons at the photomultiplier, and the 
mean response varied by about 30°% over the sensitive area. 


Lead 
shield 


Darvic box 


A Monte Carlo calculation indicated that more than 80% of the energy 
of an electron or photon would be absorbed within the counter at incident 
energies up to 300 mev, but that at energies above this the response would 
become appreciably non-linear. The calculation was used only to 
determine corrections, for edge effects, non-linearity of response, and for 
the lead shield ; which totalled about 20°, of the signal in the energy 
range 100-300 Mev per particle. The p-meson calibration, therefore, 
gave a direct lower limit for the energy associated with each incident 
electron, but the upper limit depended to some extent on cascade calcula- 
tions. 
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The counter was situated at the centre of the array described in I, 
(Cranshaw and Galbraith 1957) and signals from it were measured in 
coincidence with showers. Logarithmic amplification was used to cover 
the wide range of pulse heights encountered ; pulses were measured with 
a mechanical kicksorter, and recorded on punched cards. 

Measurements were carried out in three independent runs, and approxi- 
mately 20000 showers were analysed. In the first two runs the system 
was set to cover the range 10 Mev/metre?—25 000 Mev/metre”, and in the 
third run, 100 Mev/metre?—100 000 Mev/metre?. 


§ 3. RESULTS 
The experimental points are plotted in fig. 2, for showers in which 3 
units only were struck on the array. These showers contained on average 
3x 108 particles (see Paper III, Cranshaw et al. 1957). The line is the 
particle density distribution measured in III, and plotted on ordinates 


Energy density (MeV m °) 


1) 100 200 300 400 500 600: 
Distance from axis (metres) 


such that the mean energy associated with each electron is 210 Mev. 
The points are consistent with a distribution of this form, and show no. 
significant deviations from it. Major contributions to the errors shown 
in fig. 2 arose from uncertainties in the positions of the shower axes 
which were determined from the centroids of the units struck on one 
array ; and from the ~-meson calibration, which contained an appreciable 
contamination of knock-on particles and small air showers. No correc- 
tion has been applied for variation in the zenith angle of the showers, or 
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for penetrating particles. From previous results (III ; and IV, Porter 
et al. 1957) these effects are known to be small. 

The results are consistent with a constant mean energy of 210 + 60 Mey, 
associated with each electron observed by the array in showers. Since 
the energies carried by electrons and photons are probably almost equal at 
sea level, the mean energy of each electron appears to be about 100 Mey, 
close to the critical value of 84 Mev. This energy is within the range of 
linearity of the Cerenkov counter for both electrons and photons, but 
must be regarded as a lower limit for particles within 100 metres of the 
axis. It is unlikely that the mean energy of each electron at 60-500 metres 
exceeds 200 Mev, but it may be considerably higher at distances of 0-50 
metres from the axis, as was reported by Matano ef ad. 

A lower limit for the energy appearing at sea level in the electron— 
photon component of 3-fold showers may be set at 6x10!¢ev. The 
distribution of energy for showers striking eight or more units on the 
array and containing on average 2-7 x10’ particles, appears, although 
statistically poor, to be similar to that for 3-fold showers, and the mean 
energy of the electrons is approximately the same. 


§ 4. Discussion 


The theory of Nishimura (quoted in Matano ef al.) for electron—photon 
cascades with s=1-0, predicts a mean electron energy close to the critical 
value, at distances greater than one shower unit from the axis. It would 
be expected that this would apply also for showers with s=1-4, the age 
parameter derived from the particle density distribution in III, and 
probably for shorter distances. To this extent, the experimental results 
are consistent with the general conclusion that the electron—photon 
component is probably in equilibrium with the nucleon component, 
and has a structure characteristic of old showers. 

The estimation of the primary energy required to produce a shower of 
3 x 108 critical energy electrons at sea level, is not strongly dependent on 
the form of the model assumed for the cascade. Pure electron—photon 
cascade theory, under approximation A or B (Rossi and Greisen 1941) 
predicts an energy of 2-5 x 101° ey, provided that the number of initiating 
y-rays is less than about 10. Theories which take into account the nuclear 
cascade in the atmosphere (Rozental 1952, Oda 1957, Ueda and Ogita 1957) 
would also demand energies of the same order of magnitude. 

The effect of fluctuations in the early development of the shower must 
be important, however, if the primary particles have interaction cross 
sections comparable with the geometric value for protons (Kraushaar 
1957). A reliable estimate of the mean primary energy could be made 
only if the interaction lengths of the primaries were known. At present 
it must be assumed that appreciable contributions to the rate of showers 
of a given size at sea level, may be made by primaries varying over an 
order of magnitude in energy. 
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§ 5, CONCLUSIONS 


The distribution of energy in the electron—photon component of showers 
containing 3 x 10 particles at sea level has been measured at distances of 
60-500 metres from the axis. It has a form similar to that of the particle 
density distribution, and the mean energy of the electrons in this region 
appears to be close to the critical value for air. The majority of particles 
observed by a Geiger counter in showers of this size appear, therefore, 
to be true cascade particles, and not low energy products of collision 
processes. The estimation of the primary energy corresponding to the 
showers cannot be made reliably without knowledge of the interaction 
cross section of the primaries. 
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ABSTRAOT 


It is shown that the Griineisen constant derived from the elastic properties 
of a solid can be negative in certain cases. The lattice models examined are 
(a) ionie lattices with the zine blende structure, (b) ionic lattices with the 
rocksalt structure. In the first case it is found that a negative value is 
obtained over the whole range of parameters examined. In the second, 
negative values are obtained when the shear modulus of the lattice is suffi- 
ciently small. In all these cases the volume expansion coefficient of the 
cubic lattice would be negative at sufficiently low temperatures. 


§ 1. [yTRODUCTION 


Ir has been commonly assumed that a crystal will have a positive 
volume expansion coefficient at all temperatures, if only vibrational terms 
in the free energy are considered, and the interaction between the particles 
of the lattice has a reasonable form. This belief is based on calculations 
for simple models, as in the Griineisen theory, and on the generally good 
agreement of this theory with experiment. It is, of course, conceivable 
that a negative volume expansion can result if the free energy contains 
configurational terms, and an electronic component has been suggested as 
responsible for negative values in uranium and plutonium. 

On the experimental side the evidence is meagre. The only relevant 
cases in which negative volume expansion coefficients have been reported 
are silver iodide, zinc blende and silicon. Of these the best investigated 
instance is that of silicon (Erfling 1942) where a region below —130°c 
shows this negative value. There is so far only one investigation for zinc 
blende (Adenstedt 1936) and the silver iodide measurement, that due to 
Fizeau, is subject to some uncertainty as to the nature of the material. 

It is the purpose of this paper to show theoretically that cubic lattices 
can have a negative volume expansion coefficient at low temperatures 
when a normal type of potential interaction is employed. This arises 
from the fact that in certain types of transverse waves, the quantity 
y, = —d log v/d log V (where v is the frequency of a normal vibration and 
V the volume) can be negative. This has been noted by Dayal (1944) 
and by Barron (1955). 

The starting point of the investigation was that for an alkali halide 
of the rocksalt type, the quantity dc,,/dr is positive (see Blackman 1957) 
where ¢,, is the shear modulus in Voigt’s notation. Since at sufficiently 
low temperatures, the heat capacity and the free energy are controlled 
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by the c,, value, if this is sufficiently small compared with the volume 
modulus (Blackman 1935), it seemed likely that for a suitable model of 
this type a negative volume expansion would be found. The results of the 
investigation (§ 4) show that this is the case. On the other hand, such 
examples are rather special in that these do not fall in the experimentally 
observed region of alkali halides. The experimental results quoted above 
have one common factor, namely the coordination number four. A model 
which has this in common with the above crystals, is a purely ionic lattice 
of the zinc blende type. This model is considered below, and it will be 
shown that a negative volume expansion is found in that region where 
only the elastic constants are important, i.e. at very low temperatures. 
The character of the results is such that it is expected that the region will 
extend to moderately high temperatures. Furthermore, the results do 
not appear to be dependent on the particular force constants used for the 
repulsive forces. 


§ 2. THe ELastTic CONSTANTS OF THE ZINC BLENDE MODEL 


The elastic constants of the zinc blende crystal have been investigated 
by Born and Bollmann (1920) using an ionic model. As the aim of the 
calculations was the elastic constants in the equilibrium state, these were 
not obtained in the form required here. It is essential that the equilibrium 
conditions as such should not be used, if one wishes to investigate the 
changes of the elastic constants with volume—which is what is required 
in the calculation of the volume expansion coefficient at very low tempera- 
tures. Furthermore, the calculation of Born and Bollmann contains 
numerical errors, as was pointed out by Born (1923). 

Expressions for the elastic constants of the zine blende lattice in terms 
of lattice sums were obtained by Born (1923) by specializing his more 
general relations for the static elastic constants of a crystal. The same 
formulae were also obtained in the course of this work by starting from the 
general equations of motion of a crystal and expanding these for vibrations 
of long wavelength. The three elastic constants ¢,1, Cy, C4, can be obtained 


from the following relation (cf. also Born 1923): 
, 
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Here /= (1,, J,, 13) numbers the lattice cells, k and k’ (equal to 1 and 2) 
number the two particles in the rhombohedral unit cell of volume A, 
The quantities 2!,,,, ete. are the Cartesian components of the vector Pn 
joining particle k in cell / to particle k’ in the zero cell; P¥,,, and Q! kh ate 
given by the relations: 


Lddbny 
PEF et te ee 
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Qn = ee = ip I | 
radr\(r dy Pek J 


where ¢,,,, is the potential interaction between two particles k, k’. 

The calculation of the elastic constants} is carried out by taking the 
potential interaction as the sum of two parts, the electrostatic and the 
repulsive. The method adopted for the electrostatic part, labelled by 
the index (e), follows that due to Born and Bollmann for c,,®. That for 
€44*, which was not evaluated there, was determined by a similar method 
as a check on the calculation. It can be shown by symmetry arguments 
that c4,* = —4c,,%. The other constant c,. can then be obtained by 
relating the expression containing (c,,“ + 2c,,) to that for the Madelung 
constant of the zinc blende lattice. The calculation for C also followed 
the same lines as in the calculation of Born and Bollmann. In addition 
numerical methods were used to obtain values of c,,? and C, namely to 
serve as a general check of the method employed. The results obtained 
are listed below. In these r is half the lattice spacing of the cubic cell 
(A= 2r°) and e is the magnitude of the charge on each of the ions. 


C46? =0-1255 a 6, *O = — 09-0630 
lng) Aechmea" Te Dy 
CF 
Oye, = Tey a ol 8 
C19 1-324 (3) 
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The contribution to the elastic constants from the repulsive part of the 
potential was calculated under the assumption of an nth power law. For 
simplicity only nearest neighbours were taken into account. The re- 
maining constant can be related to the Madelung constant «(a= 1-891) 
by introducing the equilibrium distance r, of the lattice. 


+ There is a slip in notation in Born’s expressions (Born 1923, p. 548, eqn. 
(32c)) where two indices have been interchanged; the quantity defined as 
[wy&7] should be [x@y7]. This does not really affect the expressions for the 
elastic constants given by Born for cubic crystals, as the interchange will not 
alter the magnitude of the elastic constants once the equilibrium conditions 
are used. It is, however, important in our case where the elastic constants 
are subject to differentiation with respect to the lattice distance. 
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In this notation the elastic constants can be written in the form : 


1 e2 — we’7,”-1 (n—1) 7} 
[et aera 
1 e = axe*rg"—1 (n+ 5) 
cam q[ 1825 + m2 
; 1 e «aer,” 7 (n—l1) 
C* y= x | — 00885 aie yn a6 r . . (4): 


§ 3. y VALUES OBTAINED FROM THE ELasTIc CONSTANTS 


The expressions obtained for the elastic constants of the zine blende 
model allow one to obtain the quantities y, =d log v/log V for the vibrations. 
of sufficiently low frequencies. In this region of the vibration spectrum 
the y, values will depend on the direction of propagation of the elastic 
waves of a given type, as distinct from the frequency itself. The 
expression for y, can be transformed into one depending on the velocity 
of the elastic waves (Barron 1955) : 


1 9 fide 1 [ d(pc?r) 
es = al SS sae 5 
UE Ges ( dr i 6pc?|_ dr {c (6) 


where c is the velocity, p the density of the solid, and r the lattice spacing 
defined above. 

The quantity pc? can be obtained, from the equations of motion of the 
elastic continuum, in terms of the elastic constants and the direction of 
propagation of the wave considered. For a given direction of propagation 
there are three solutions for pc? and these solutions are of a particularly 
simple form when the direction of propagation coincides with particular 
directions in the solid associated with a high degree of symmetry. 

To obtain numerical values one has to choose a value for the 
parameter n. To illustrate the results n=6 has been taken. The 
numerical quantities used in the calculation are quoted below: 


eyo LITO, - (cyr) = —8 fae | 

Caa'g = 0-644 - , = (Cyy" — yp") = 5-78 ma > | 

Cap? = 0-987 = = (Cag) = 2-97 o po ake 
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The derivatives are, of course, taken at r= ro, the equilibrium value. 
It is of interest to consider the special cases for which pc? = ¢41, C44, $(¢1, — C42) 
since these are associated with directions of special symmetry. 

The three y values are 1:24, — 0-766, and — 5-09 respectively. Two of 
these, which have negative values, are associated with transverse vibrations 
which have also relatively low velocities. In the weighting of the y values 
to form the correct average value, this point may be important. It is the 
fact that the relatively highly weighted vibrations contribute negative 
terms to the average y value which can result in an overall negative contri- 
bution, even if the magnitude of the contributions were not as large as 
those obtained above. 


g t v1 Per'%o Y2 pls"To Ys pls lo 
1:0 | 0-0 | 0-841 | 1-724 | —5-10 0:0945 | —0-769 | 0-644 | 

0-2 | 0-823 | 1-751 | —4-47 0-107 —0-801 | 0-608 

0-4 | 0-798 | 1-808 | —3-33 0-135 —0:-965 | 0-521 

0-6 | 0-774 | 1-857 | —2-72 0-178 —1-132 | 0-423 

0-8 | 0-759 | 1-898 | —2-05 0-233 —1-503 | 0-333 

1-0 | 0-757 | 1-909 | —1-75 0-278 —1-75 0-278 

0-8 | 0-0 | 0-843 | 1-708 | —4-21 0-112 —0-769 | 0-644 

0-2 | 0-823 | 1-740 | —3-82 0-124 —0-808 | 0-601 

0-4 | 0-788 | 1-806 | —3-18 0-156 —0-951 | 0-502 

0-6 | 0-769 | 1-864 | —2-47 0-199 —1-173 | 0-401 

0-8 | 0-759 | 1-895 | —2-05 0-234 —1-504 | 0-325 

0-6 | 0-0 | 0-909 | 1-640 | —2-91 0-180 —0-769 | 0-644 

0-2 | 0-822 | 1-682 | —2-40 | 0-189 —0-869 | 0-592 

0-4 | 0-788 | 1-765 | —2-11 0-208 —0-987 | 0-490 

0-6 | 0-772 | 1-883 | —2-08 0:217 —1-148 | 0-414 

0-4 | 0-0 | 0-988 | 1-558 | —1-192 | 0-316 —0-769 | 0-644 

0-2 | 0-855 | 1-567 | —1-212 | 0-309 —0-830 | 0-589 

0-4 | 0-804 | 1-679 | —1-358 | 0-274 —0-991 | 0-511 

0:2 OV ae 1-02 1-308 | —0-715 | 0-512 —0-769 | 0-644 

0-2 | 0-927 | 1-405 | —0-738 | 0-456 —0-859 | 0-603 

0-0 | 0-0 | 1-244 | 1-176 | —0-769 | 0-644 —0-769 | 0-644 


To obtain a survey of the y values for various directions in the crystal, 
a system of coordinates introduced by Barron and Domb (1955) is em- 
ployed, in which the components q, 72, 73 of the wave vector q are related 
in the form 7, =sq>, 73=1q2. Because of the cubic symmetry of the crystal, 
it is necessary to consider only one-sixth of an octant in q space. The 
results are set out above. Instead of the velocities, the value of c?pr, is 
given, the units of this being e?/2r,°. 
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It will be seen that two of the solutions of the equation give negative y 
values in all cases. The average value of »(y=yo) is obtained from the 


expression : 
3 


“yf dQ 
10a = en i= 


3 1 dQ (7) 
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where the index j numbers the three solutions and the integrals are taken 
over the solid angle. In the notation used dQ =(1 +s? +4)? ds dt. 
The numerical integration gives : 


Yo= ery) s P A - 5 5 : (8) 


This negative value is not confined to the special value of n=6. In fact 
it is always negative for all values of n which have been considered, e.g. 
for n= 10, yo is larger in magnitude than for n= 6, and this latter case was 
chosen because the y, value was relatively small. 

The constant y,) derived from the elastic constants is the Griineisen 
constant at low temperature, i.e. in the 7? region of the heat capacity of 
the solid. The only other region for which a y value can be obtained with- 
out a detailéd examination of the vibrational properties of a crystal is at 
relatively high temperatures. Here the y values depend on the evaluation 
of moments of the frequency spectrum. The auxiliary quantities y(p) are 
calculated (Barron 1955) where 


1 d log p, 
- p (a eens MP aig et 
y(p) Lys / 2 p dlog V S 


and , is the pth moment. The summation is taken over all frequencies 
of the lattice. The limiting value of y at high temperature, y,,, is equal 
to (0). 

The quantities derivable from the secular equations for the frequencies 
of vibration of a lattice are y(2), y(4), ete. while yy=y(—3). In order to 
obtain y(0) an extrapolation is needed, which is uncertain where the 
y(s) values do not lie on a monotonic curve. 

The value for y(2) can be calculated very simply for the zinc blende 
model used above, as the quantity $v? can be obtained directly from the 
determinant for the frequencies of the lattice. 'The coulomb terms in the 
expréssion for }v? vanish for a cubic lattice and the contribution of the 
nearest neighbour repulsion terms is easily evaluated. The result is 
y(2)=$(n+2). For the case n=6 it follows that »(2)=1-33. The extra- 
polation to y(0) would involve the evaluation of y(4) ete., and this has not 
been thought worth while at this stage. It seems, however, clear that 


the y value is positive at high temperatures as contrasted with the negative 
value at very low temperatures. 


§4. Tue ALKALI Hatipe MopEn 


As has been pointed out in §1, the formal expressions for the elastic 
constants of an ionic alkali halide lattice suggest that the y, value can be 
negative in extreme cases. As has been shown by Barron (1955) and 
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Blackman (1957), there are negative y values associated with transverse 
vibrations whose velocity depends on the elastic constant c,,. Typical 
values are yO = —0-23 for a model with the repulsive constant n=7 and 
y= —0-74 for n= 10 (these values have been unfortunately interchanged 
in the paper by the present author quoted above). The corresponding 
quantities for other important cases are all, however, positive. It appears 
that the negative value is larger in magnitude as n is increased, but this in 
itself is not sufficient to suggest a negative value for the average value. 
It seems clear that one would have to rely on the weighting factor which 
occurs in the expression for yp (eqn. (7)). Since the square of the velocity 
is proportional to an elastic constant, it follows that if c,,/c,, is sufficiently 
small, the negative y values associated with the shear waves characterized 
by the modulus c,, would give such a weighting. The ratio c,,/c,, can be 
made as small as one wishes, by increasing the value of n. This leads to an 
increase in ¢,,, leaving ¢,, unaltered in the model used. While a value of 1 
greater than 10 would not be in the range found with the use of experi- 
mental data, there is no reason why such values should not be considered 
as long as the lattice remains stable. 

The value actually chosen for the calculation was n=21. Here 
C44/Cy,= 9-06. With such a small value of ¢,,/c,, it is to be expected that 
the y values would vary markedly with direction and that a numerical 
method of averaging, as in § 3, would be unsatisfactory. It is, however, 
possible to use an approximation discussed elsewhere (Blackman 1935) 
in which the non-diagonal terms in the determinant for the velocity of 
elastic waves are neglected. The velocity c can be written in the form 


pe? = C1, p* + Cy,(G? +177) Oa We ae LO) 


together with two similar solutions in which the directions cosines p, g, 7 of 
the wave vector are interchanged in cyclic order. This approximation 
allows the integrals for y, to be worked out numerically. The solutions (10) 
are, however, not sufficiently reliable to be used in the calculation of y 
values at those points where they yield the same values of pc?. These 
have been studied separately and are found to produce a negligible effect 
on the value of the integral. 
The value of y) obtained with n= 21 is: 


ied OL0 AU debe? Weis. <2 vi nso LL) 


The calculation of the value of y,, was carried out using the known expres- 
sions for y(2) and y(4); and these lead to y(2)=4-00 and y(4)=46. An 
extrapolation taking into account the yp) value yields y,,=3-4. 

The variation of the average y value for this alkali halide model is 
unusually large, as is to be expected in view of the very small ¢4,/¢,,._ The 
point at which the y value will change from negative to positive values as 
the temperature is raised cannot be deduced from the above calculations. 
The general behaviour of the contour surfaces for the frequencies of 
vibration of a lattice with a highly anisotropic elastic behaviour suggests 
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that the negative value which holds in the 7° region, will be replaced 
rapidly by positive values as the temperature is raised. 

The alkali halide model forms, however, the simplest case in which the 
existence of a negative volume expansion coefficient can be demonstrated, 
and is one in which this effect can also be seen to follow from qualitative 
arguments once the formal expressions for the elastic constants are taken 
as given. 


§ 5. GENERAL COMMENTS 


The calculations described above are restricted to certain ionic lattice 
models and do not in themselves apply to specific practical cases. There is 
no direct relation, for instance, to the case of silicon, and even for zinc 
blende itself a purely ionic model would not necessarily be a good approxi- 
mation. It would, however, be of particular interest to have additional 
confirmation of the results on zinc blende as well as measurements on 
other substances at low temperatures. The calculations suggest that 
negative volume expansion coefficients are to be expected from open 
structures rather than close packed ones and that those with relatively 
low shear moduli would be favoured. One of the interesting theoretical 
questions, which is still open, is whether the negative volume expansion 
coefficients are restricted solely to temperatures which are small compared 
with the Debye temperature—in particular whether y,, can be negative. 
This does not appear to be excluded. While it is true that the quantity 
y(2) will always be positive, it is possible to conceive of cases where this is 
smaller than unity and where »(4) is appreciably larger than y(2). Ifthe 
y value is negative and sufficiently large in magnitude, the y,, value could 
well be negative. This argument is based on an extrapolation procedure. 
If one considers the definition of y,, as the average value of y for all vibra- 
tions, then it would be necessary for a sufficiently large number of these 


to be negative. This again does not appear to be beyond the realm of 
possibility. 
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ABSTRACT 


It has been found that when one of the alkali metals is brought into 
electrical contact with some other metals like aluminium, cadmium, 
magnesium, ete., the contact resistance of the junction undergoes remarkable 
variations under the influence of an electrical stress. The exact nature of the 
variation depends primarily upon the contact-combination of the metals. 
However, the nature and magnitude of the electrical pulse and polarizing 
voltage, the pressure applied to the surface junction, etc., play important 
roles in determining the finer characteristics of the variation. 

The present paper gives a brief account of our investigations on the 
potassium—aluminium surface combination only. A tentative explanation 
of these phenomena has been suggested on the basis of formation and disrup- 
tion of an intermediate oxide layer under appropriate electrical conditions. 


$1. INTRODUCTION 


THE influence of Hertzian waves on the contact resistance of metallic 
particles has long been recognized. Branly (1890, 1892) found that fine 
metallic filings contained in a tube of glass or ebonite between two 
metallic electrodes offered an enormous resistance of about 10’7ohms to 
the flow of current from a single Daniell Cell ; the resistance diminished 
to some hundred ohms or even less, by exciting at a distance, a spark 
from a Leyden jar, or simply by putting a high voltage momentarily in 
the circuit. The original resistivity was immediately restored by a very 
slight mechanical shock on the tube containing the filings. Minchin 
(1894) and Lodge (1894) also reported on the sudden acquisition of 
conducting power by a series of discrete metallic particles by the action 
of Hertzian waves. Lodge (1899) called this instrument a ‘ coherer’ 
and attributed its property to the electrostatic attraction between 
particles. He found however that a few metals (e.g. silver) showed an 
increase of resistance under similar electric influence. Bose (1899, 1900) 
showed that alkali metals, notably potassium, also showed an increase 
of resistance under the influence of Hertzian waves, followed by an 
automatic recovery. These facts could not be interpreted satisfactorily 
by any contemporary theory. However, owing to the rapidity with 
which wireless science has progressed, new methods of detection super- 
seded the old coherer, before the implications of the latter were fully 
understood. 
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In recent years, the importance of electrical contacts in power and 
communication engineering has led to considerable interest in the physics. 
involved in their action. Holm (1958) and Llewellyn-Jones (1957) have 
discussed the fundamental processes of electrical contact phenomena 
in great detail. Bowden and Tabor (1950) have given a general account 
of the nature of contact between two solid surfaces ; while Williamson 
et al. (1956) have investigated the influence of dust particles on the 
contact of solids. Bowden and Williamson (1954) have investigated the 
effect on the nature of contact between metals of the passage of short 
duration electric pulses of sufficient intensity to cause softening in the 
region of contact. They pressed two pieces of gold together with a force 
of a few grams, and obtained typical values of the true area of contact 
and constriction resistance of 10->cm? and 10-3ohm. This work shows 
that oxide films, if any, were thinner than 10 A, and good conduction took 
place according to quantum mechanical tunnel effect. Wilson (1955) has 
studied the contact resistance and mechanical properties of surface films 
on metals, by measuring the electrical contact resistance and the 
coefficient of friction. His study was confined to the behaviour of metals 
when in sliding contact but covered a wide range of conditions. 

This paper describes experiments, carried out rather interruptedly 
since 1952, in which the conductivity of the surface of contact of 
dissimilar metals has been studied under the influence of transient 
electrical stress. They were undertaken with the idea of obtaining fuller 
information about the ‘ cohering’ action of a surface combination of 
metals exposed to Hertzian waves. Subsequently, however, it was found 
that a single voltage pulse from a discharging condenser across the 
contact was essentially similar to that of a group of Hertzian waves. 

The metal potassium was selected because of its low melting point 
(63°C), high electro-positive character, and its abnormal behaviour as a 
coherer (Bose 1900). All the other metals, Al, Cd, Mg, etc., were chosen 
for their affinity for K and also their readiness to form oxide layers under 
normal atmospheric conditions. Since each combination of metals has 
its own intrinsic characteristics, we shall for the sake of simplicity restrict 
ourselves in the present paper to the discussion of the essential features 
of the contact between potassium and aluminium surfaces. It may be 
added that the observations described hereunder were consistent and 
reproducible. 

The properties of other combinations will be taken up in proper 
sequence. 


§ 2. EXPERIMENTAL INVESTIGATION 
2.1. Haperimental Arrangement 
The contact arrangement usually consists of a rectangular block of 
aluminium and a truncated cone of potassium. <A sphere of aluminium 
pressed against a rectangular block of potassium provided an equally 
efficient junction. The holder and leads are made of silver. Analytically 
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pure metals are generally used. The junction is kept immersed in pure 
paraffin oil in a suitable glass container. It is equally satisfactory if the 
paraffin oil is replaced by argon gas. 

The aluminium surface is cleaned by fine grade carborundum paper, 
while the potassium surface is placed inside the oil and cut by means of 
a sharp blade. A mechanical device is constructed by means of which 
the potassium wedge can be lowered vertically on to the aluminium block. 
The micro-manipulator also provides minute adjustments of pressure at 
the contact of two surfaces without any sliding or rotational movement. 


Fig. | 
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The junction is placed in an electrical circuit as shown in fig. 1. The 
magnitude of the polarizing voltage at the contact can be varied by 
means of a potentiometer, while its sign can be altered by a reversing 
switch. An automatic recording of the current through the junction is 
made by a pen recorder. 

The Hertzian waves are produced by sparking between two platinum 
beads and an interposed platinum sphere 6mm in diameter. By a single 
sudden break of the primary of an induction coil, a flash of radiation can 
be emitted consisting of a Hertzian wave packet. The spark-gap is 
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placed at the end of a brass tube 5cm in diameter having longitudinal 
metallic grids at the outer end while the coil and its battery are enclosed 
in a copper box. The unit serves as a portable generator of damped 
polarized Hertzian waves of about 5mm wavelength. The intensity of 
the incident wave-packet can be varied by placing the radiator at different 
distances from the junction. 

The single voltage pulses from charged condensers are injected across 
the junction by means of conventional switchgear. 


2.2. Rectification 


When a freshly cut potassium surface is brought into electrical contact 
with a lightly ‘ aged ’ aluminium surface (an hour’s standing of a freshly 
scraped surface in paraffin oil, preferably in contact with potassium, 
usually suffices), the contact, on proper adjustment of pressure, acquires 
strongly rectifying properties for polarizing voltages ranging from about 
0-1 to lv. The resistance in the forward (or easy) direction is of the 
order of 20ohms, while in the reverse (or blocking) direction it is about 
100000 ohms, in the equilibrium state. The latter situation corresponds 
to potassium having a negative potential with respect to aluminium. 
The resistance in both directions are, more or less, voltage dependent. 


2.3. Induced Conductivity in the Blocking Direction 
(K negative, Al positive) 

If exposed to the electromagnetic influence of a spark-gap oscillator 
the contact suddenly loses most of its resistance and, with it its rectifying 
properties. The current through the junction, even for the blocking 
direction, jumps up to almost the maximum value allowed by the 
associated circuit at that polarizing voltage. The magnitude of the 
polarizing voltage (denoted in this paper by P.V.), seems to have little 
bearing on this sudden loss of resistance, except in special circumstances 
to be discussed later. While the recovery progresses even when the 
P.V. is reduced to zero, the rate of recovery is much enhanced if a higher 
P.V. is applied, or if the original P.V. is stepped up during the course 
of recovery as shown in fig. 2. 

The effect is more clearly revealed by the current vs time curves ; 
sample curves are shown in fig. 3. 

These response-cum-recovery curves are obtained with the junction 
maintained at different polarizing voltages (with K negative) and 
subject to weak irradiation. The enhancement of recovery rate at higher 
P.V. becomes more clearly manifest, if the P.V. is increased during the 
course of recovery, as shown in fig. 4. 

However, if the voltage injection at the terminals of the contact is 
increased by bringing the radiator nearer the junction, the recovery 
curves become more flat-topped and sometimes show signs of permanence. 
Figure 5 shows the recovery curves at various P.V. when the injection 
is strong. If, however, the flow of current is temporarily interrupted by 
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switching off the battery, there is a marked enhancement of recovery rate 
as shown in fig. 6, in which the ‘ switch-off ’ period is shown to scale. 
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2.4. Induced Conductivity in the Forward Direction 
(IX positive, Al negative) 

If now the terminal of K is made positive with respect to Al, a basic 
current flows through the junction, the rectifying contact being in the 
forward direction. Injection of voltage by electromagnetic radiation 
usually leads to a further fall of the junction resistance, followed by a 
rather jerky recovery to its original value. It is interesting to note, 
however, that full recovery of current in the forward direction does not 
necessarily indicate corresponding recovery in the blocking direction. 
This is exemplified in fig. 7. 

The equilibrium resistance in the forward direction increases as the 
P.V. in this direction is gradually lowered. Finally, at a value of P.V. 
of about 0-05v., the response-cum-recovery curves become identical for 
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both directions. Under such circumstances, the recovery is slow but 
regular. 

All the observations given above are related to the usual type of 
contact, with a pressure of a few tenths of a gram across a thin oxide 
layer. The nature of the response-cum-recovery curve radically changes 
if the aluminium surface is coated with a thick oxide layer by special 
treatment. The resistance for such a junction is usually higher. Injection 
of voltage by electromagnetic waves now causes a steep fall of conductivity, 
practically to zero. The following recovery is equally sharp, and can be 
almost instantaneous if the P.V. (in the forward direction now) is 


Fig. 7 
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sufficiently high (about 1v). At lower P.V. the recovery-jump is some- 
what delayed and at P.V. of about 0-2 or less, there may be no recovery 
at all. The normal recovery curve rises sharply to about four-fifths of 
the original value, and then climbs up rather slowly. A few typical 
negative response curves are shown in fig. 8. 

This type of contact is highly rectifying but shows no response to 
irradiation if the P.V. is reversed to the blocking direction. The normal 
positive response (elastic increase of conductivity) can be obtained by 
increasing the pressure of the contact, but this again leads to the 
disappearance of the negative characteristics referred to above. In fact, 
the positive and negative response characteristics do not appear to 
coexist under identical contact conditions. 

It is found that the observed induced conductivity cannot be ascribed 
to any absorption of 5mm electromagnetic wave quanta like photo- 
electric effect. Rather, it is caused by the voltage induced across the 
surface junction. In order to test this point, pulses have been injected by 
discharging a condenser across the junction. It is found that a pulse of 
—3v or +1:5v produces similar increase of conductivity in the blocking 
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direction, followed by usual recovery, though the rate of recovery is 


slower in the former case. 
In fig. 9, the effect of small pulses injected across the junction during 


the recovery periods are depicted. A positive pulse causes a small rise 
Fig. 8 
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of conductivity, while an equivalent negative pulse produces a correspond- 
ing diminution. The saturation effect due to stronger voltage injection 
by the negative pulse of —3v may also be noted. In all the cases, the 
basic P.V. has been kept at — 0-4 v. 


2.5. Voltage/Current Characteristic Curves 


A variety of voltage/current relations have been observed on potassium— 
aluminium junction. They depend on the previous history of each 
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particular contact. The feeble current (about 10-¢amp) flowing through 
the junction in the blocking direction usually obeys Ohm’s law within the 
voltage limits zero to 0-5v. However, if this limit is exceeded, the 
current increases in short jumps. If the P.V. is now gradually reduced 
to zero, the return path takes place along a route of higher conductivity. 
Resistance loops occur in the forward direction also. Figure 10 shows 


Fig. 10 
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voltage/current characteristics of a particular specimen in the reverse 
direction. The measurements have been made under static conditions. 
The specimen maintains constant resistance up to —0-5v and then 
suddenly diminishes if the voltage is increased. However, if the voltage 
is gradually diminished, the new resistance remains stable. During next 
outward excursion, it goes up to —0-7v and subsequently undergoes a 
cross-over before assuming another value. Such abnormal hysteresis 
phenomenon has been reported by Henisch and Granville (1951) under 
dynamic conditions only. More pronounced loops appear if the measure- 
ments are made before the specimen has recovered from its previous 
breakdown. 


§ 3. Discussion oF EXPERIMENTAL RESULTS 
3.1. Oxidation of Metals 

It is well known that with the exception of gold, no pure metal 
(including platinum) and no alloy is stable at room temperature. All 
metals tend to form oxides under normal atmospheric conditions. Recent 
experimental work indicates that oxidation is initially extremely rapid, 
but that the rate soon drops to a very low or negligible value, a stable 
film 20-1004 thick being formed. An explanation of this behaviour was 
first given by Cabrera and Mott (1949). It is based on the hypothesis 
that a strong field is set up in the oxide film due to a contact potential 
difference between metal and absorbed oxygen, which enables the metal 
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ions to move through it without much help from the temperature. The 
theory gives a logarithmic growth law of the type : 


i 

gd = Bint» ite ene ieante 
X being the thickness in time t. Thus a stable film will grow until it 
reaches a limiting thickness if the temperature is low enough. The 
influence of temperature, gas pressure, texture and orientation of surface 
layers on further growth of oxidation layers have been discussed in 
details by Kubaschewski and Hopkins (1953). 

Some of these oxide films are insulators and some are semi-conductors ; 
in the first case the contact ceases to function, while in the second case 
a rectifying action may be produced. The formation of a rectifying 
barrier between two oxide layers of dissimilar metals of different work 
function has been shown by Kittel (1953). A non-rectifying high resist - 
ance barrier is usually presented by two oxide layers of the same metal 
in loose contact with each other. 


3.2. Action of Coherers 


A normal coherer consists of oxidized metallic powders of the same 
kind lightly pressed together between two metallic electrodes. Ordinarily, 
the combination behaves as a high resistance unit. The incident Hertzian 
wave induces a potential difference between the electrodes and the 
individual particles, resulting in an electrostatic attraction between them 
given by 


F' =(V2/87l?) x 9 x 104d ynes em PORE ee se | 


where F/ is the mechanical force per unit area of attraction due to a 
potential difference V volts between conductors /cm apart. 

If V=1, and /=24A, this force exceeds the yield stress of a typical 
metal (=10*°dynescm-). 

Thus the metallic powders are further compressed to admit current by 
tunnel effect. Once the current has been established, the various micro- 
contacts are heated up to yield plastically, forming larger areas of contact 
with consequent lowering of resistances in accordance with the theory 
of Bowden and Tabor (1950). The effect of the mechanical shock is to 
break the closeness of contact and restore the originally high resistivity. 


The maximum temperature @,, in a monometallic contact is related to 
V by the relation 


BOn?=0,°)= V3j4) eee ee 


where @» is the maximum temperature at the contact and 6, is the 
general temperature outside the constriction region (both measured in A) 
and B is a constant=2-3 x 10-§ ohm watt (°c)-. Slepian (1926) 
applied the above results to metallic contacts and showed that contact 
voltages of one volt should produce contact temperatures of about 
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2000°c to 3000°c. Such high temperatures can produce thermionic 
emission, thermal ionization, evaporation, melting and boiling. This fact 
lends support to Lodge’s theory of electrical welding which explains in a 
simple manner the diminution of contact resistances of various metallic 
filings when subjected to strong electric field variation. 


3.3. Induced Conductivity : Physical Theory 


Many rectifiers are based on the Mott-—Schottky rectifying barrier 
formed between a metal and a semi-conductor. In fact, in the case of 
K—AI junction, the barrier seems to correspond with the electron energy 
diagram showing rectification at a metal and a n-type semiconductor 
junction, K being a metal and Al,O, a semiconductor. In the blocking 
direction (K negative), the height of the barrier prevents electron flow 
from metal to the semiconductor. However, because of the extreme 
thinness of the barrier layer (about 10-*cm), the injected pulse creates a 
field of the order of 10° volts/cm, resulting in field emission as suggested 
by Mott and Gurney (1948), It also causes the breakdown of the covalent 
bonds of Al,O, and the formation of holes. The positive holes will be 
attracted to the barrier region, and if enough of them are trapped there, 
the summit of the potential hill would be pulled down. There will now 
be a relatively free flow of electrons from potassium to aluminium, which 
will continue until the abnormal concentration of positive holes disappears 
through recombination. Taking into account the case of a. barrier at the 
junction of a metal and a p-type semiconductor whose height should 
increase due to hole injection, Chatterjee and Sen (1955) have tried to 
explain the behaviour of * negative ’ coherers of the type introduced by 
Bose. 

The lifetime of free electrons and holes in solids has been analysed 
by Rose (1957). The enormous range of lifetimes encountered in solids 
extends from 10~!%sec to several seconds. The lifetime of induced con- 
ductivity of K—Al junction is generally of the order of a few seconds, 
although under strong voltage injection, it may extend to several minutes. 

The nature of the response and recovery curves, however, suggests that 
the phenomenon of induced conductivity and its consequences may 
perhaps be more satisfactorily explained on the basis of chemical trans- 
formation of the intermediate oxide layer under suitable circumstances. 


3.4. Chemical Theory 


Since the metal K and Al are quickly covered with thin oxide layers, 
the contact combination consists of the following disposition :— 


KO. AON rr) 


Thus, an intermediate layer of potassium aluminate, KAIO, is formed, 
which partially ionizes as |K*+A10,~|, giving rise to basic conductivity. 


P.M, 3M 
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On the application of polarizing voltage of either sign to the electrodes, 
the following chemical changes are likely to occur :— 
(A) Blocking direction (K negative, Al positive) 

Electrons are available to the KAIO,, and facilitate the following 
reactions : 


2K A110, = DT DAS + Al,Os + Ome | 
(5) 


30---+ Al>2A1,0,+ 6e- 


(anode) J 
Thus Al,O, forms as a high resistance layer (of about 10°ohms) in the 
blocking direction. 
(B) Forward direction (K positive, Al negative) 
The cations and anions K+ and AlO,~ migrate to their corresponding 
electrodes to combine into : 


K++ Al+e—>K. Al 


(cathode) (Alloy) ( 6 ) 
and K+ Al0,->KAI0,+ €- 
(anode) (Potassium Aluminate) J 


which are conducting films. The actual resistance of about 20ohms may 
be due to traces of unconverted A1,O, films. 

Since the current has been found to be equal in either direction at the 
P.V. of about 0-05 v, it seems to be the ionization potential of potassium 
aluminate salt. 

The situation may perhaps be compared with the physical theory of 
barrier rectification which necessitates a barrier voltage greater than 
kT'/e (greater than 0-026 v at ordinary temperature) to enable the junction 
to function as a rectifier. 

The various features of induced conductivity in either direction may 
now be interpreted in the following manner. 


Case A. Blocking Direction (K negative, Al positive) 

(a) The injected pulse causes field emission, resulting in microarcs and 
high temperature at the junction. Under such appropriate conditions, 
Al,O; may be reduced to metallic Al (Thorpe 1937). The reformation of 
of Al,O, layers is indicated by recovery curves shown in figs. 2, 3 and 4. 

(6) Increase of P.V. during the course of recovery accelerates the 
formation of Al,O; as shown in figs. 3 and 4. 

(c) At high P.V. (about 0-8v), a weak pulse injection leads to positive 
response with sharper recovery (fig. 3), as expected. But stronger pulse 
injection leads to a ‘ saturation effect ’ as shown in figs. 5 and 9. In this 
case the intense pulse of the current causes a conducting path of aluminium 
to be formed between the metals, possibly by the mechanism discussed in 
case A(a). This bridge will be relatively large because the amount of 
charge passing between the surfaces was large ; the surfaces may even 
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be ‘ welded ’ together in the sense used by Bowden and Williamson (1954). 
The temperature then falls because the large bridge passes the small 
polarizing current easily. The very slow recovery then results from the 
slow oxidation of the aluminium bridge in the manner discussed by 
Fairweather (1953). 


Case B. Forward Direction (K positive, Al negative) 


(a) Partial positive response and recovery curves (fig. 7) may be due 
to the combined effect of the reduction of residual oxide layer and 
increased formation of KAIO, under suitable conditions. 

(6) The negative response and recovery curves shown in fig. 8 are 
rather interesting. It may be noted that in this case the oxide layer is 
too thick (in the blocking direction) for the injected pulse to trigger field 
emission. Consequently, the function of a sinusoidal pulse would be to 
facilitate chemical reactions (5) and (6) in corresponding half cycles. 
However, on account of the greater mobility of electrons, the rate of 
reaction (5) involving the formation of Al,O, is quicker and prevails over 
reaction (6) forming KAIO,. This explains the steep decline of conduc- 
tivity caused by the induced a.c. voltage. The subsequent rise of current 
is due to the exponential formation of KAIO, in the forward direction. 
As the P.V. is lowered, the formation of K AIO, is delayed. The straggling 
curve R in fig. 8 probably represents the initial stage in the formation 
of K AIO, under lower P.V. 


Case C. Mixed Type of Response 


The effect of single voltage pulses of both signs and different magnitudes 
on K—AI junction has been shown in fig. 9. A P.V. of —0-4v has been 
used in these experiments. However, the effective P.V. across the 
junction changes depending upon its conductivity, the remainder being 
distributed across the recording instrument. 

It may be noted that a positive pulse of +3v gives rise to a positive 
response, followed by a sloping recovery. In this case, the transient 
pulse acts in the forward direction, assisting the formation of K—A] alloy 
and KAIO, accompanied by field emission. The recovery indicates 
gradual reformation of Al,O;. The second pulse +1-5v causes a small 
rise of conductivity. The fourth and fifth pulses of —1-5v show that 
small negative pulses produce a decline of conductivity. The sixth pulse 
of —3Vv causes an intense field emission followed by a slow formation 
of Al,O, as suggested in Case A (c). It has also been found that a positive 
pulse is more effective in producing positive response than an equal 
negative pulse. However, a detailed analysis of these phenomena is 
reserved for a subsequent communication. 


§ 4. CONCLUSION 
The behaviour of a surface combination of metals like K and Al has 
been studied by simultaneous measurements of current and resistance 
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under a wide range of conditions. The observation of the induced 
conductivity due to injected pulses seems to afford a powerful method 
for studying the surface states of metals. The results are in accord with 
the hypothesis of the formation and disruption of an intermediate oxide 
layer under appropriate electrical conditions. The term oxidation has 
been used in the general sense of any surface process that is accompanied 
by a delivery of electrons to the metal. 
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ABSTRACT 


A summary is given of the published data of electrical resistivity of semi- 
conductors which exhibit ordered spin arrangements below a Curie or Néel 
temperature. It is found generally that electrical anomalies accompany the 
advent of magnetic order only if the specific resistivity at the transformation 
temperature is less than a critical value. The implications of this rule 
are discussed in the light of recent theories of electrical conduction in 
spontaneously magnetized materials. 


§ 1. INTRODUCTION 


THE relationship between resistivity and magnetic order in ferromagnetic 
and antiferromagnetic metals and alloys has recently been discussed by 
Mott and Stephens (1957) and by Coles (1958). The most remarkable 
conclusion reached by these authors is that the anomalous reduction in 
resistivity below the magnetic transformation temperature is caused by 
distinctly different mechanisms in two groups of materials. On one hand 
the resistivity of nickel and alloys based on this metal is described in 
terms of conduction (4s) electrons which may be scattered into a partly 
filled 3d band under the influence of crystal lattice vibrations. On the 
other hand the anomalous behaviour of other metals is described in terms 
of the disorder of the exchange potential experienced by conduction 
electrons interacting with those electrons which are localized at lattice 
points and are mainly responsible for the magnetic properties (Kasuya 
1956). 

There remains a large family of compounds known to be mainly ionic 
in nature and to have ordered spin arrangements below a Curie tem- 
perature. These compounds contain at least one constituent element 
belonging to the first transition group and one non-magnetic element. 
They exhibit a wide range of magnetic and electrical properties. The 
former properties include examples of ferro-, ferri- and antiferromagnetic 
ordering ; electrically these materials range from semiconductors with 
relatively high activation energies to good conductors with practically 
metallic behaviour. The present paper discusses the available data with 
special reference to the influence of magnetic order on electrical 
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conductivity ; a novel, and hitherto overlooked general relationship between 
electrical properties and magnetic ordering emerges from these considera - 
tions. The new law is then discussed in the light of the theory put forward 
by Mott and Stephens and elaborated by Coles. 


§ 2. SuRveY or Data 


In order to determine the influence of magnetic order on electrical 
resistivity, p, one requires data of p as a function of temperature extending 
over an appreciable range on either side of the Curie point. Such data 
are available for compounds based on the chemical formulae NiO, CoO, 
Fe,03, Laq_»Sr,MnO;, MnTe, MnSe, Mn,Crq_,Sb, FeTiO, and MeFe,0,, 
where a denotes a proper fraction and Me a divalent metal or a mixture 
of divalent metals. 


pOnnaee Approximat 
Miah Magnetic | Conduction | transformation he igtivi = " 
interaction type temperature ie y 
(°K) ohm em 
NiO anti I 520 107 
CoO anti i 291 4x10? 
HesOs anti i! 950 108-107 
Lavo.98t(o.y MnO; ferro II 46 4 
Lay <o.9)8t(> 0-1) Mn0, ferro III 480-530 10-1_10-2 
MnSe anti I 160 10° 
MnTe anti IMME 307 1 
MnSb ferro Ill 578 2:6 10=4 
CrSb anti Til 690 aq On 
FeTi03(0.g)Fe.0 30-2) ferri x Pad 224. 3 x 10 
FeTiO 3(0:5)H'e20 30-5) pa t I os a 
Fe,0, ferri Il 858 Sie Wie 
NikFe,O, ferri I 858 102 
CuFe,O, ferri I 728 10 
MnFe,0, ferri I 573 104 
MgFe,0, ferri [ 713 10 


The normal oxides of the three ferromagnetic metals are themselves 
anti-ferromagnetic, their transformation temperatures being given in 
the table. The nature of the magnetic interaction is believed to be due to 
‘super-exchange ’ in which the oxygen ion plays a dominant part. This 
has been discussed by Anderson (1950) and van Vleck (1951). Careful 
electrical measurements have been reported by Morin (1951, 1954) and by 
Verwey et al. (1950). The stoichiometric oxides in each case behave like 
simple semiconductors and are accurately described by the usual equation 
p=Aexp (</kT), with A and € constant over a wide range of temperatures. 
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There is thus no detectable anomaly at the Néel point. A theoretical 
discussion of the conduction process has been given by de Boer and 
Verwey (1937) and by Mott (1949). 

The mixed lanthanum-strontium manganites have been found to be 
ferromagnetic by Jonker and van Santen (1950). The same authors have 
also reported electrical measurements (van Santen and Jonker 1950). 
The main features of their results are reproduced in fig. 1. It will be noted 


Fig. 1 


(la,_y3 Sr) Mn 0; 


1000 
a 


log p vs 1/T for lanthanum-strontium manganites. Curie temperatures 
"are indicated by arrows. (After van Santen and Jonker 1950.) 


With acknowledgment to the Editor of Physica. 
that there is a gradual transition from semiconducting LaMnO, to ne 
metallic properties with x=0-3 to 0-5. A small anomaly at the C Hate 
point is found for «=0-1 in the form of a change in slope of the log p vs 1/T 
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characteristic. This anomaly becomes more pronounced with increasing 
strontium content. For all the mixed La,Sr manganite specimens 
reported by Jonker and van Santen the advent of magnetic order is 
accompanied by a reduction in the resistivity from that extrapolated 
from the paramagnetic region. A theoretical discussion of the nature 
of the magnetic interaction and its bearing on the electrical properties 
has been given by Zener (1951). 

Manganese telluride and manganese selenide are both antiferromagnetic 
and resistivity measurements show no anomaly at the Néel temperature 
for the latter material (Palmer 1954) while the former shows the distinct 
anomalous reduction of resistivity with magnetic order (Uchida et al. 1956) 
like that described for the mixed La,Sr manganite system. 

The system Cr,Mn,_,Sb has been investigated by Suzuoka (1957) and 
is of particular interest as it undergoes a transformation from ferro- 
magnetism to antiferromagnetism as v exceeds 0-75; the crystal structure 
remains that of the NiAs type for all compositions. A marked anomaly 
at the magnetic transformation temperature exists for all values of 2. 
The electrical resistivity of ferromagnetic manganese antimonide is similar 
to that of the manganites discussed above while the antiferromagnetic 
chromium antimonide resembles the conduction behaviour of « manganese 
(White and Woods 1957) in shape and in magnitude. 

Most ferrites exhibit ferrimagnetism below their magnetic transforma- 
tion temperature the nature of the magnetic interactions being discussed 
very fully by Néel (1948) and Gorter (1954). 

An electrical resistance anomaly at the transformation temperature has 
been observed so far without doubt only for the well-conducting ferrous 
ferrite, Fe,0, (see for example Smith 1952). Bochirol (1951) and Parker 
(1957) have reported a linear relationship between logp and 1/7 for a 
large number of ferrites including those of Zn, Me, Co, Cu and mixtures 
of Ni-Zn, over a temperature range which includes the Curie temperature. 
On the other hand, Komar and Kliushin (1954) have found very slight 
anomalies in manganese ferrite and in nickel ferrite. The author has 
recently made careful measurements on a number of specimens of these 
materials both with d.c. and a.c. up to several Me/s and is unable to 
confirm the observations of Komar and Kliushin. Typical results 
obtained for nickel ferrite and a mixed nickel-zine ferrite are shown in 
fig. 2. The Curie temperature is indicated by arrows. 

Measurements have recently been reported for the system 


vHeTiO; — (1—x)Fe,O, 


which gives rise to both ferromagnetic and antiferromagnetic compounds 
depending upon composition and thermal treatment of specimens. As in 
the case of ferrites no conductivity anomalies accompany the ordering of 
magnetic spins. 

The results described above are summarized in the table. I refers to 
a linear log pvs 1/7 relationship throughout a temperature range 
including the Curie or Néel point as illustrated in fig. 2. II refers to a 
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small change in slope in this relationship. III denotes a large anomaly 
with the material exhibiting semiconducting behaviour above and pseudo- 
metallic behaviour below the transformation temperature. (Examples of 
modes II and III are shown in fig. 1.) The last column shows the 
approximate value of resistivity at this temperature. 


Fig. 2 
(B) 
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(A) 

Q 
2) 
a 
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log p vs 1/T for (A) Ni ferrite fired in air at 1350°c, (B) Ni,Zn,Fe,O,, fired 
in air at 1150°c. 


§ 3. THEORETICAL CONSIDERATIONS 


There is at present no satisfactory explanation even of a qualitative 
nature of the results summarized in §2. A theoretical treatment of the 
conductivity of magnetic compounds has been attempted by Heikes 
(1955). Described in its briefest outline it is as follows. In materials 
coupled by super-exchange interaction a given electron can drift down 
the crystal lattice under the influence of an externally applied field only if 
suitable unoccupied energy levels are available at each lattice point along 
a row of atoms. This will only be the case along rows of atoms where the 
3d electrons are coupled ferromagnetically. According to this view 
ferromagnetic compounds should all exhibit metallic conduction, with 
anomalies in resistivity at the Curie point similar to ferromagnetic metals. 
Antiferromagnetic compounds should all be insulators, conduction being 
possible only as a result of thermal activation to excited energy states. 
This model is unsatisfactory on several counts. Firstly, even in an 
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antiferromagnetic lattice rows of atoms can be defined along which the 
spin alignment is parallel. Electrons could travel along such directions 
without temperature activation so long as some overlap in the ionic wave 
functions exists. Secondly the disappearance of short-range order above 
the Néel temperature should be accompanied by a resistivity anomaly 
as pointed out by Heikes himself. This is not in accord with the published 
data on pure antiferromagnetic oxides. Thirdly, since the appearance of 
Heikes’ paper ionic antiferromagnetic compounds have been discovered 
(CrSb for example) which are good conductors and thus violate his rule. 

Heikes’ theory has also been criticized by Mott (1956). Mott has shown 
that the electrical conductivity of compounds of transition metals is 
more convincingly described in terms of the overlap of 3d electronic wave 
functions of the transition metal ion with the p electrons of the negative 
ions than in terms of the long-range magnetic order. A 3d shell less than 
half full is associated with a larger overlap integral (and hence with a 
higher conductivity) than a shell containing 5 or more 3d electrons. Mott’s 
theory is in accord with the resistivity values listed in the table, except 
for the case of manganese telluride and for those compounds known to 
have positive ions of different valencies. 

Neither of the two theories described in this section appears, at the 
present stage, to account satisfactorily for the resistivity anomalies which 
accompany the magnetic transformation in some of the compounds listed 
in §2. The present author is unable, at this stage, to offer any explanation 
for these properties. He wishes however to draw attention to the 
relationship between the type of temperature dependence of conduction 
and the numerical value of the resistivity at the magnetic transformation 
temperature. All materials for which this resistivity exceeds 4ohm cm 
exhibit type I conduction while a value of p less than 0-6 ohm em is always 
associated with type III. Between these two values apparently either 
type is possible and one example of type Il is recorded (which may be 
considered as the transition between I and III) with a Curie point 
resistivity of 4ohmcm. It appears therefore as a general rule that a 
resistivity anomaly at the magnetic transformation temperature exists 
only if the resistivity of the material is already below a well-defined 
threshhold value. It is remarkable that this rule is obeyed by compounds 
differing widely in crystal structure, the nature of the magnetic order 
and even in the currently postulated origin of the exchange interaction. 

It may well be that the true transition region is defined even more 
sharply than indicated by the results in the table. It is known that 
measured values of specific resistivity of these materials may differ from 
the true values in the crystal lattice by a factor of 2 or 3 unless the most 
stringent precautions are taken to eliminate errors due to contact 
resistance (Hlaschen and van Uitert 1956) and due to high resistivity 
surface layers between the boundaries of crystal grains (Volger 1954). 

A further remarkable confirmation of our empirical law is obtained 
from the results of Morin for impure or non-stoichiometric oxides of 
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nickel and iron. It has already been pointed out that the pure oxides 
are highly resistive at the Néel temperature and conform strictly to type I 
conduction. The resistivity of nickel oxide can, however, be greatly 
reduced by firing the material in an atmosphere of oxygen and thus making 
the compound deficient in nickel. A similar reduction is obtained by 
‘ doping’ the oxide with a small quantity of lithium. In either case it 
can be noted in fig. 3 that an anomaly appears in the log 1/p vs 1/T' relation - 
ship when the resistivity at the Néel temperature is less than 5ohm cm; 
this anomaly corresponds to the conduction type I. The resistivity of 
ferric oxide can be reduced in a similar manner by doping with titanium. 
Unfortunately the measurements for ferric oxide do not extend over a 
large region above the Néel point (950°K) but nevertheless the results 
show clearly a strictly linear relationship between logp and 1/7’ for 
specimens which have a Néel point resistivity in excess of 5ohm cm, and 
a marked anomaly of type LI for specimens of lower resistivity. 

The results on ferric oxide and nickel oxide just described suggest 
strongly that the cause for the anomalous resistivity behaviour at the 
Curie temperature must be thought in the increase in concentration of 
the conduction electrons rather than in an increase in their mobility. 
One may conclude therefore that magnetic ordering affects the con- 
ductivity of a material only when the concentration of conduction 
electrons rises above a sharply defined value. 


§ 4, Discussion or Resistiviry ANOMALIES IN 
MaGnetic Merats anp ALLOYS 


The considerations of §3 have an important bearing upon the theory 
recently put forward by Mott and Stephens and by Coles. According to 
these authors the anomalous reduction of resistivity below the magnetic 
transformation temperature is explained in the following manner. In 
nickel and alloys based on nickel half the conduction electrons cannot be 
scattered into the 3d band and their mobility is therefore enhanced to a 
value approximating to that of conduction electrons in copper. In other 
metals and alloys the conduction electrons are not scattered by a spacially 
disordered exchange potential field when the material has become 
spontaneously magnetized. Neither of these mechanisms appears, to a 
first approximation, to depend upon the concentration of electrons which 
partake in the conduction process. It appears therefore that the current 
theories of electrical conductivity in magnetic metals and alloys are quite 
inadequate in accounting for the observed properties of a large number 
of magnetic materials. 
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ABSTRACT 


The theory of small angle scattering from dislocations (Atkinson and 
Hirsch 1958) has been extended to take into account the splitting of disloca- 
tions into partial dislocations. At very small angles the splitting does not 
change the scattered intensity, but its effect becomes important as the angle 
increases till at large angles the intensity is reduced to 3 of its value at 
zero angle. It is suggested that it may be possible to determine the widths of 
extended dislocations from the shape of experimental scattering curves for 
deformed metals. 


§ 1. INTRODUCTION 


In a previous paper (Atkinson and Hirsch 1958, to be referred to as I) 
a theoretical investigation was made of the small angle scattering of 
X-rays or neutrons from cold-worked polycrystals containing random 
networks of unextended dislocations. However, in face-centred cubic 
metals, dislocations are split into two partial dislocations bounding a 
ribbon of stacking fault (e.g. Seeger and Schoeck 1953), and it is of 
interest to know in what way the splitting of the dislocations might affect 
the intensity of small angle scattering. Seeger (1956) has calculated 
the scattering from a single extended dislocation assuming that the 
density change is due entirely to a second-order elasticity effect ; however, 
this term, and the term due to the core of the dislocation, have been shown 
in I to result in much less scattering in the range of scattering angles for 
which experimental data were available than that due to first-order 
elastic effects from edge dislocations. The purpose of the present paper is 
to calculate the effect of the splitting of dislocations in terms of the 
first-order theory. The notation used will be the same as that given in I. 


§ 2. THEORY 


Consider an extended dislocation with total Burgers vector b at an 
angle % to the dislocation (fig. 1). The Burgers vectors b,, by, of the 
partial dislocations have magnitude 6,=b//3 and make angles (jb + 47) 
and (%—¢7) with their respective dislocation lines. Hence the two partial 
dislocations have edge components of magnitude b,sin (+4) and 
b, sin (s— 477) respectively. Assuming that all the scattering is due to the 
As order elastic effect, then in the notation of the previous papel see 
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eqn. (I, (19)), the amplitude scattered by one extended dislocation is 


crystal b 


A,f 6... ors 
A(k)= a | A iz sin (b+ $a) + — sin (xh — bar) exp (ik w)| exp («k.P) dv 
(1) 
where A is the dilatation produced by an unextended edge dislocation 
given by eqn. (I, (2)) and w is the vector distance between the partial 
dislocations. (It is assumed that w<R, the radius of the crystal.) Note 


Fig. 1 


partial dislocations 
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“1 stacking fault 


An extended dislocation. 


that the ribbon of stacking fault itself results in no density change and 
therefore does not contribute to the small angle scattering. Since w is 
not a function of z, we can proceed as before and write, corresponding to 
(I, (21)) and (I, (22)) 
v5 +H 
A=— | exp (izk,) dzX ,’ Sal napeesee he) 
J_H 
where 


DERE DGS iz sin (+ $77) + 4 sin (xb — dar) exp (th,w cos a} eh tees) 
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After integration with respect to z, we obtain 


A N\2 sin Hk,\? 
oa. ae 2 A He GS Sn ee een 
[A| (=) ue Hk, ) |X." (4) 


Now consider a random array of dislocations in a polycrystalline 
specimen : The dislocations are supposed to have no preferred direction 
of orientation in their slip planes; hence any particular set of dislocations 
parallel to a chosen direction in the specimen will contain dislocations 
with all possible Burgers vectors; the average intensity from such a set 
will therefore be given by an expression similar to (4), where |X,’ is 
replaced by its value averaged over all possible Burgers vectors, i.e. 


PP=2f Pa 
i) 


Substituting (3) in this equation we obtain 


2 7 
BOA Be (3) 7+ | 2sin (s+ ¢77) sin (xb — 477) cos (k,w cos €) i . (6) 
w 0 
The intensity per unit volume, /,,, scattered from the whole random 
array of dislocations is then found by averaging the scattering from the 
sets of parallel dislocations (mentioned above) over all possible directions 
with respect to s: 7, is then given by the expression, which corresponds 
to eqn. (I, (30)) 
[oe eee 
= — Cea Pa ee et 230 Gi 
: 47s QQ? rae e Pads (7) 
where L is the total length of dislocation line per unit volume (i.e. the 
dislocation density). 
In the particular case for unextended dislocations, (6) reduces to 


|X’? =5X 2. In I the dislocations were arbitrarily divided into edges 
and screws, and the symbol Z in eqn. (I, (35)) or (I, (36)) represents the 
density of edge dislocations only. The present more sophisticated 
analysis leads, for unextended dislocations, to the same result as (I, (35)) 
or (I, (86)) except that LZ should be taken as half the total dislocation 
density ; i.e. if the dislocations take up random orientations in their slip 
planes, the effective edge dislocation density is equal to half the total 
dislocation density. 


|X,’ and hence J, are now calculated for extended dislocations. Two 
cases are considered: In the first w is taken to be constant, equal to w,,, 


independent of , and in the second w is assumed to vary as expected 
from the elastic theory approximation, i.e. 


Wad == 7 COS 2us) ve aa nS) 


where ¢=2v/(2—v) (see Read 1953, p. 131 eqn. (9.8)) and Wy is the width 
of a dislocation making an angle of $7 with its Burgers vector, 
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In the first case 
[X,'P=3X.%{1 +4008 (k,w,cosé)}. . . . . . (9) 
In thé second case 


|X. |? =4X .2{1 + 4 cos (kw, cos £) J 9(qk,wy cos €) 
— sin (k,w cos €) J 1(qk,wy cos €)t. =. (10) 
In the first case we obtain, after substituting (9) in (7), putting 
X,.=27iF sin €/k,, see eqn. (1, (34)), and integrating 
L F? 
eae imsqp te 8 5 Lee ee eet, Pea 11) 


where #, =(1+J,(a,,)/a,) and a,,=2asw 


m m* 
In the second case the integral involved in (7) is 


-2 1 = 2m : 
be |X,’P?dé= —— a | {1+ 4 cos (k,wy cos €)J 9(qk,w, cos €) 
0 


de 


—sin ae cos £)J,(gk,wycos€)}sin®Eédé. . . . (12) 


The integral containing J, can be solved by using the Jacobi expansion 
(e.g. Gray and Mathews 1895, p. 18) for the cos (k,w, cos €) term and an 
expansion for J, in terms of squares of Bessel functions (e.g. Gray and 
Mathews, p.28); the term involving J,(= —J,’) is first integrated by 
parts to give an integral containing J); part of this latter integral is 
similar in form to the above integral; the other part can be solved in a 
similar manner again using expansions for the terms J, and sin (k,w, cos €) 
appearing in the integrand. The resultant integral is found to be 


L 1M 


Dip Fr 
J. |X. Pdé= 37 Bs, and therefore J, = —— A,"47 4 —=H, . (13) 


drs Q? 
where 
By= 1+ (5-2) {Jalal Jarchag) + Sede} 
J (a) [J 9°(20q) + 25 17(Saq) + J 2*(309) | 
+ J 4(a)[J4?($aq) + 2F9°(5aq) + J3°(30q) ] 
rae} 
+ = {(0 a) Fo°($aq) — J*($0q) 1+ F(a) [417(3a9) — J2*(3aq) | 
J5(a)[Jo2(4aq) — Je2($ag)J+...} . . . (14) 


and a= 27s. 

E, and E, have been evaluated as a function of a, and are shown in 
fig. 2. Curve (a) represents #, for a value of v=4; curves (b) and (c) 
a. EB, for w,,=Wp, and w,,=1-38w, respectively. The latter value of 

m Das been chosen so that E, fits #, closely. It is clear that for practical 
ee the simpler expression of eqn. (11) is adequate, Seah i 
the value w,, =1-38w) is used. The curves show that for values of aS 3 
the intensity decreases to about } its value for a=0. The effect of the 
splitting of the dislocations on the intensity is therefore a small one; 
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the neutron scattering experiments on copper (Atkinson 1958) were 
carried out at a mean value of s~0-014—, and taking wy as ~ 8-56 
(Seeger and Schoeck 1953), a ~ 1-4 and hence the effect due to splitting 
is to decrease the intensity by about 13%. If the scattered intensity 
could be measured as a function of a in the range in which H, or H, is 
varying, it might be possible to determine the widths of dislocations in 


this way. 


Lene [ee eee 
O 5 1O 


Q= 21TSWo 


Curves of H,, #, as a function of a. (a) H, for v=% (solid line). (6) EH, for- 
Wm=W, (long dashes). (c) H, for w,=1-38w, (short dashes). 
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ABSTRACT 


The conditions are described under which thin metal specimens can, be 
produced by electro-polishing for examination by transmission electron 
microscopy. The method has been applied to Mg, Co, Cu, Ni, Al, «-brass 
and Fe, and representative micrographs showing the arrangement and 
movement of the dislocations are presented. Some practical hints are 
given relating to the use of the transmission electron-microscope technique 
for the study of dislocations in metals. 


§ 1. INTRODUCTION 


Recent work has shown that it is possible to see dislocations directly 
in thin metal foils by transmission electron microscopy. A number of 
different techniques have been used to prepare specimens : by beating 
and etching (Hirsch et al. 1956), by electro-polishing (Bollmann 1956, and 
Whelan et al. 1957), by evaporation (Pashley et al. 1957), and by electro- 
deposition (Jackson, unpublished). The electro-polishing technique is 
particularly useful for the study of dislocation arrangements in deformed 
bulk specimens, which are thinned after deformation. This technique 
has now been applied to a variety of metals, and it is the purpose of the 
present paper to describe the details of the method and to record the 
optimum polishing conditions applicable to the various metals. Repre- 
sentative transmission electron micrographs are included to illustrate 
the types of photographs which might be expected. 


§2. PoLisHinc TECHNIQUE 


The technique described by Bollmann (1956) was found to be suitable 
for stainless steels, where a very viscous polishing solution is used. For 
all the other metals examined so far the following method is employed. 

The initial material in the form of a strip of metal sheet or foil, 0-001 in. 
to 0-0lin. thick, is coated with a varnish leaving a half inch square 
‘window’ on either side of the foil. The varnish is a solution of poly- 
styrene in trichlorethylene which can be used without danger in a 
perchloric acid polishing bath. 

The polishing bath itself is very simple, consisting of a 150 ml glass 
beaker surrounded with ice and two-thirds filled with polishing solution. 
A one inch wide strip of the appropriate material (see the table) placed 
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in this beaker is used as cathode. A variable voltage is supplied to the 
bath through a potentiometer circuit across a 110 volts d.c. mains supply 

The specimen serving as the anode is suspended vertically in ihe 
bath from a spring-clip and is perpendicular to and 4in. to 1 in. away 
from the cathode. The foil is not mounted rigidly, as a slight amount of 
agitation (by hand) helps to produce a uniformly polished surface, and 
facilitates the rapid removal of the specimen from the bath after polishing 
before etching takes place. , 

The foil is polished at a voltage corresponding to the high voltage 
end of the plateau on the current-voltage curve. At this point the 
current just begins to rise rapidly for a slightly increased value of the 
resistance in parallel with the cell, while the voltage across the cell 
remains almost constant. This point can be quickly determined for 
each specimen. It is found that polishing under these conditions gives 
specimens with the least amount of surface structure, although this 
appears to vary with the orientation of the grains. The time taken for 
polishing is from 30sec to 10 min depending on the thickness of the 
specimen, and on the current density employed ; the polishing is usually 
stopped when about half the * window ’ has been polished away. (fig. 1). 


Fig. 1 
Cut microscope 
specimen from here 
Varnished metal 
surface 
ae) 


Appearance of ‘ window’ in specimen after polishing. 


The current is then switched off, and the foil is very quickly placed 
in a dish of methyl alcohol to remove the polishing film. In the case 
of nickel and aluminium the foil is then dried on filter paper ; the other 
metals, which readily oxidize, are left in clean alcohol until required, 
but microscope specimens should be made within an hour after polishing. 
These may then be kept under vacuum conditions if necessary. 

It is found that if the current through the polishing bath is switched 
on and off rapidly for the last few seconds of polishing, specimens are 
produced which are uniformly thin over rather larger areas ; further, the 
remaining polished part of the foil is left with jagged edges, from which 
it is very easy to cut a microscope sample. 
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The specimens are mounted for examination in the electron micro- 
scope as follows. The polished foil is placed between two sheets of 
tissue paper and a suitably sized area, approximately one square milli- 
metre is cut away with a razor-blade from the edge of the foil remaining in 
the ‘window’. The specimen is mounted on a 200 mesh copper erid 
coated with a 1°, solution of polybutene in xylol, and the ‘ polished ’ edge, 
as distinct from the ‘ cut’ edges of the foil, is arranged at the centre of 
the grid. Sometimes it is convenient to sandwich the foil between the 
200 mesh grid and a 100 mesh grid, as this prevents the specimen from 
becoming detached from the support during insertion through the airlock 
of the microscope. 

It is usually possible to cut three or four specimens from the ‘ window ’. 
If necessary, the remaining part of the ‘window’ can be polished to 
make further specimens. The extent of the thin regions is rather variable, 
sometimes the specimen is uniformly thin over the whole area of a grid 
square (fig. 2 (a), Pl. 47), or more, but more generally, useful areas are 
found around the polished edge or around holes, and these areas may be 
10-100 » wide. On the average about 50°, of the specimens are of this 
type and suitable for study of dislocation arrangements and movements 
over considerable areas. 

Little difficulty has been encountered in finding suitable conditions 
for thinning the metals so far examined by transmission microscopy, and 
this suggests that the method could readily be extended to other metals 
and alloys. 

The metals thinned by electro-polishing, the composition of the baths, 
and the polishing conditions are summarized in the table. It will be 
noticed that the solutions used are not always those generally recom- 
mended, (see for example McG. Tegart 1956) which may give a better 
surface for optical examination but do not produce such uniformly 
thinned specimens. 

Typical micrographs from metals other than stainless steel (see Whelan 
et al. 1957) are shown in figs. 2-8, Pls. 47-49. These photographs show 
the type and quality of the micrographs which should be obtained. No 
attempt is made here to discuss the details of the dislocation arrange- 
ments and movements observed. Systematic studies of the dislocation 


arrangement as a function of deformation of a number of these metals 
will be published elsewhere. 


§ 3. Hints on THE MicroscopE TECHNIQUE 
The following points may be of interest to readers wishing to use the 
transmission technique. 


(1) The electron microscope should be operated at the maximum 
possible voltage (i.e. 100 kv for the Siemens Elmiskop I) when trans- 
mission is optimum. 

(2) It is essential to be able to tilt the specimen through a few degrees. 
In the Siemens Elmiskop I this can be achieved with the aid of the 
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stereoholder. The contrast at a dislocation is critically orientation- 
dependent, and it will be found that maximum contrast is often obtained 
when a strong low-order extinction contour is near the dislocation. 
Sometimes dislocations can be observed in apparently dislocation-free 
areas after tilting the specimen. The significance of this orientation- 
dependence of the dislocation contrast in relation to the contrast mecha- 
nism will be discussed elsewhere. 

(3) In areas of large curvature of the foil the extinction contours 
(Heidenreich 1949) appear as narrow bands. These must not be confused 
with dislocation lines and can always be recognized by their movement 
on tilting the specimen. 
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ABSTRACT 


Previous investigations of the effect of temperature on the static and 
dynamic strength of steel are briefly reviewed. A description is given of 
tests in which small specimens of a medium-carbon steel were compressed at 
high rates of strain and at temperatures ranging from 15°c to —121°c. 
Dynamic stress-strain curves obtained from these tests are compared with 
stress-strain curves obtained from tests at normal rates of strain, at tempera- 
tures from 15°c to —183°c. The results indicate that the effect on the 
stress-strain curve of increasing the average strain rate to 600 sec is 
comparable to that of lowering the test temperature to —183°c, the upper 
yield stress being multiplied some 24 times and the lower yield strain being 
increased considerably. The effect of testing at high rates and low tempera- 
tures is to raise the curve still further. 

Metallographic examination shows that the slip occurring in the ferrite at 
low temperatures and normal rates is finer than that obtained at room 
temperature ; also there is less multiple slip. These characteristics have 
also been observed at high strain rates. No twinning took place in the 
static tests, even at the lowest temperature. Im the dynamic tests, twins 
were not formed at room temperature, but were formed at the lower test 
temperatures. 


§ 1. IyTRODUCTION 


THE strength properties of materials at low temperatures are of interest 
not only because such conditions may be encountered in practice, but 
also because a study of the variation of properties with temperature may 
yield information concerning the laws governing the processes of 
deformation in general. Much of the work done in this field has been 
limited to slow loading rates ; but in recent years increased attention has 
been given to dynamic properties, and to the possibility of correlating 
reduction in temperature with increase in strain rate. 

Among researchers primarily concerned with finding static properties 
of iron and steel at low temperatures have been Hadfield (1905), McAdam 
and Mebs (1943), Eldin and Collins (1951) and Cottrell and Churchman 
(1951). The results obtained by these investigators indicate that in 
general the static stress-strain curves of iron and iron alloys are raised by 
lowering the test temperature. 

One of the first investigations into the effect of temperature on the 
dynamic properties of steel was performed by Zener and Hollomon (1944), 
a ee ee 
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who used temperatures down to that of liquid nitrogen. They observed 
that at constant strain rates the yield stress of a medium-carbon steel 
increased as the test temperature became lower. Clark and Wood (1949), 
Krafft et al. (1953), and Krafft (1956) have also investigated the variation 
in the static and dynamic upper yield stresses of iron and steel at 
temperatures down to —196°c. Their results indicate that the time of 
loading and temperature have a marked effect on the strength of iron and 
the softer ferritic steels. 

Although the above researchers, and others, have investigated the 
stresses required to cause yield under conditions of slow or rapid loading 
and low temperatures, it appears that only Baron (1956) has obtained 
dynamic stress-strain curves at low temperatures. He carried out tensile 
tests on a number of metals and alloys at temperatures of 20°c, —78°c, 
and —196°c, using strain rates of 10?sec-! and 10-%sec!. His results, 
like those of other investigators, not only indicate a considerable increase 
in the upper yield stress at low temperatures and high strain rates, but 
also show that under these conditions twinning is an important 
mechanism of deformation. 

The dynamic loading machine used in Baron’s tests was designed to 
permit measurement of tensile stress-strain curves up to fracture. It 
was, therefore, not suitable for an accurate determination of yield stress 
or of the first part of the curve where the plastic strains are small. No 
account was taken of the effects of wave propagation in the specimens or 
dynamometer. It is well known that at high rates these effects give rise 
to considerable uncertainties in the stress measurements, particularly 
when the material exhibits a sharp yield point. 

The tests described in the present paper were made in order to obtain 
accurate information concerning the yield stress and early post-yield 
behaviour of steel at low temperatures and high strain rates. The 
specimens were subjected to rapid compression, and stress-strain curves 
derived from a consideration of wave propagation in the apparatus. 
These curves are compared with those obtained in conventional tests at 
ordinary ‘ static ’ rates of strain. 

Since the stress-strain behaviour of metals is dependent on the 
mechanism of deformation, metallographic examination of the specimens 
was made and micrographs obtained ; these permit a comparison of the 
modes of deformation occurring under various conditions of temperature 
and strain rate. 

§ 2, APPARATUS 
2.1. Static 

The static testing machine used was a 10-ton compression machine 
accurate to 1%. To minimize eccentricity of loading the load was 
applied through ball bearings and high-tensile steel plungers (fig. 1). 
The specimens, }in. in diameter and 4in. in length, were lapped in contact 
with the plunger faces and centred by means of a brass ring. 
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At room temperature, a cylinder gauge reading to 10-4in. was used 
to measure the strain. This gauge gave the relative displacement of the 
plunger faces, i.e. numerically one-half of the specimen strain. As this 
gauge could not be used for the low-temperature tests a special gauge 
was designed ; this gauge is also shown in fig. 1. The plungers fitted in 
accurately turned recesses that had holes in the centre to permit the 
specimen to bear against the plunger faces. As the specimen was deformed 
the rings of width $in. were squeezed in, the greatest strain in each ring 
occurring at A and B. The strain at these points was measured by 
resistance strain gauges placed in a Wheatstone bridge circuit so as to 
give maximum sensitivity, i.e. the strain gauges on the outside of the 
rings were in series in one arm of the bridge, and those on the inside of the 
rings were in series in the other arm. 


Fig. 1 
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This low-temperature gauge was calibrated at room temperature and 
the calibration assumed to hold at the lower test temperatures. As the 
Wheatstone bridge circuit described above is temperature-compensated 
the room-temperature calibration is valid at the low temperatures 
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provided the resistance strain gauge factor remains constant. The 
Appendix shows that at least down to —121°c such an assumption is 
very nearly correct. Further evidence in favour of this assumption is 
that in the low-temperature tests the permanent deformation, as 
calculated from the room-temperature calibration of the gauge, compared 
well with that obtained from micrometer measurements before and after 
test. 

To calibrate the gauge, known specimen strains were simulated by 
squeezing the plunger bearing faces together with two micrometers. 
A set of readings of strain and changes in the bridge balancing resistance 
was thus obtained and hence a calibration curve was plotted. This curve 
was approximately a straight line and did not differ by more than 6% 
from the calculated line. Although the strain at points A and B was 
considerably less than the specimen strain, the accuracy with which this 
strain could be measured by the Wheatstone bridge circuit was such that 
the low-temperature strain gauge effectively measured the relative move- 
ment of the plunger faces to approximately 10~‘in., i.e. to the same 
order of accuracy as that given by the cylinder gauge. 

After the gauge had been placed in position, the specimen centralized 
and a copper—constantan thermocouple bound around the specimen, two 
tufnol boxes were bolted to the steel block (fig. 1). The container so 
formed was partially filled with liquid oxygen, and the specimen tem- 
perature was kept constant by maintaining a slow feed of liquid oxygen 
to the container from a constant head reservoir. This method was found 
to work satisfactorily; if the temperature did alter slightly it could be 
brought back to its original value by altering the flow of liquid oxygen 
to the box. 

2.2. Dynamic 

The impact machine used in the present investigation was that described 
in a previous paper (Campbell and Maiden 1957). The stress caused by a 
falling weight was transmitted to the specimen through a high-tensile 
steel weighbar. The specimen temperature was lowered by using liquid 
oxygen, the arrangement being similar to that used in the static tests. 

As the specimens, }in. in length and Zin. in diameter, were cooled by 
conduction, there was a temperature gradient along the axis of each 
specimen. These longitudinal temperature gradients were investigated 
and found to be small : at two arbitrary average test temperatures the 
variations from the average were + 5° at —107°c and + Tati 41°C. 1 
was therefore assumed that the variations in temperature throughout the 
specimen in a given test would be small enough not to cause any 
significant error. 

§ 3. SPECIMENS 

The steel used was a medium-carbon steel and had the following 
composition : 

C Si Mn S P 
032% 0-235% 0:-62% 0:037% 0:025% 
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The method of machining the specimens from the original 2}-in. 
diameter bar and the heat treatments to which the steel was subjected 
have been described in an earlier paper (Campbell and Maiden 1957). 
The 4-in. diameter specimens used in the static tests were the remainder 
of a batch prepared for an earlier investigation (Campbell and Duby 1957), 
whereas the 2-in. diameter specimens used in the dynamic tests formed 
part of a new batch of specimens. 


§ 4, EXPERIMENTAL RESULTS 


4.1. Static 
4.1.1. Stress—strain curves 
Static stress-strain curves were determined at temperatures of 
+15°c, —41°o, —84°o, —121°c and —183°c. At all these temperatures 
the rate of straining before yield was of the order of 10~* sec. The 
room-temperature stress-strain curve was found by using the cylinder 
gauge to measure the strains, whereas the low-temperature stress-strain 
curves were obtained by using the specially designed gauge. The stress— 
strain curves obtained are shown in fig. 2 
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Static stress-strain curves. 


Figure 2 shows that the whole of the stress-strain curve is raised as the 
test temperature is lowered. At —183°c the upper yield stress has 
increased to 106000 |b/in?, 2-4 times the room-temperature value. After 
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yield the ratio of the stresses, at the same strain, decreases slightly as the 
strain increases. Also the stress-strain curves show that not only does 
the lower yield stress increase as the test temperature is lowered but also 
the yield-point strain becomes greater. 


4.1.2. Micrographic examination 


A number of specimens with prepared flats were tested statically at 
various temperatures. By examining the prepared flats before and after 
test it was possible to compare the mechanisms of deformation at the 
various test temperatures. Figure 3 (a), Pl. 50, shows the microstructure 
of a prepared surface after 4% static deformation at +15°o. It is seen 
that deformation occurred by coarse slip in the ferrite, the slip lines often 
being forked; some multiple slip is also evident. Figure 3(b), Pl. 50, 
shows the microstructure of another specimen after 4°% static deformation 
at —183°c. The magnification in both micrographs is x1400. A 
comparison of figs. 3(a) and 3(6) shows that at —183°c the specimen 
deformed by a much finer slip than at room temperature. In addition 
there is much less multiple slip at —183°c than at room temperature. 
Examination of specimens tested at —41°c, — 84°c and — 121°c indicated 
that the lower the test temperature the finer the slip by which the steel 
deformed. 

Examination of the prepared flat after deformation at — 183°c failed 
to show the formation of any twins. As these might have been disguised 
by slip, two specimens were polished and etched after they had been 
deformed 4% at —183°c. When examined, these specimens also failed 
to show the formation of twins. Similarly the specimens used to find the 
static stress-strain curves at —41°c, —84°c and — 121°C were examined. 
Twins were not observed in any of these specimens, and it was therefore 
concluded that twinning was not a mechanism of static deformation even 
at the lowest test temperature used. 


4.2. Dynamic 


4.2.1. Dynamic stress—strain curves 


To obtain a dynamic stress-strain curve corresponding to a given impact 
velocity and test temperature, the relevant elastic, weighbar and anvil 
records were analysed in the manner described by Campbell and Duby 
(1956). The elastic records were obtained from the weighbar strain 
gauges which were centred 3in. above the specimen upper face, the elastic 
record being obtained with a high-tensile steel dummy in place of the 
specimen. The anvil bar records were obtained from strain gauges on the 
anvil bar, centred 3in. below the specimen lower face. 

Dynamic stress-strain curves were found for up to three impact 
velocities at each of four different temperatures of test. The impact 
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velocities used were 430in/sec, 470in/sec and 503 in/sec, and the tem- 
peratures +15°c, —41°c, —84°c and — 121°c respectively. At —121°c 
yield was found to occur only at an impact velocity of 503in/sec and 
hence only the dynamic stress-strain curve corresponding to this 
particular impact velocity is presented. For the same reason, at — 84°C 
only the stress-strain curves at 470in/sec and 503 in/sec were obtained. 

In general two or more weighbar and anvil bar records were obtained 
at each temperature and impact velocity. Under identical conditions the 
records obtained at each gauge position were found to be consistent to 
+1% at room temperature and approximately +2% at lower test 
temperatures. The average curve for each gauge position was used to 
obtain the dynamic stress-strain curve. 

Figure 4, Pl, 51, shows two typical stresstime oscillograms, each of 
which has been superimposed during enlarging on the relevant recorded 
calibration lines. The traces are interrupted every two microseconds in 
order to obtain an accurate time base. 

It was assumed that for a given impact velocity the elastic stress—time 
waves at the different test temperatures would differ very little from one 
another ; justification for this is given in the Appendix. Hence in finding 
the dynamic stress-strain curve at a particular temperature the elastic 
record used was the average of those found at room temperature (15°C) 
at the impact velocity under consideration. 

The derived dynamic stress-strain curves, together with the static 
stress-strain curves at the temperatures concerned, are shown in 
figs. 5(a)-(d). In these figures, the stress plotted is the mean of the values 
derived for the upper and lower specimen faces; these values were found 
to differ by not more than 8% at any given time. The strain plotted is 
the mean strain in the specimen at the time concerned (the time, in micro- 
seconds, being indicated on each curve). The values of permanent strains 
shown were obtained from micrometer measurements of the specimens 
before and after impact. 

The table on p. 881 has been drawn up to compare the dynamic stress— 
strain curves shown in fig. 5. In this table the following quantities are 
tabulated: test temperature, impact velocity, time to yield, dynamic 
upper and lower yield stresses, ratio of dynamic to static upper yield 
stress at the temperature concerned, maximum and average strain rates, 
and permanent deformation produced. 


4.2.2. Micrographic examination 


Several specimens were prepared for micrographic examination before 
impact and examined after test. The predominant mechanisms of deforma- 
tion at room temperature and at —41°c appeared to be very fine almost 
indiscernible slip and a certain amount of localized shear near some of 
the grain boundaries. These characteristics have previously been 
observed in steel impacted at room temperature (Campbell et al. 1957) 
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Fig. 5 (continued) 
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Stress-strain curves obtained at various temperatures. Test temperatures 
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Neither in specimens impacted at room temperature, nor in speci- 
mens impacted at velocities of 430in/sec and 470in/sec at — 41°C, were- 
twins observed. However, in specimens impacted at 503in/sec at — 41°C 
and in all impacts at —84°c and —121°c, twins were formed. The 
number of twins evident in specimens impacted at 503in/sec and —41°c, 
and in specimens impacted at — 84°c was very small; but at —121°c the 
number of twins formed was appreciably greater. 


§ 5. Discussion 

The results given in the table show that for equal impact velocities (and 
hence approximately the same elastic stress pulse) the time to yield and 
the upper yield stress increase steadily as the temperature is lowered. 
These two quantities are not independent, but are related by the form 
of the elastic stress pulse applied to the specimen. The dynamic lower 
yield stress also increases steadily as the temperature is lowered. Its value 
is presumably affected by the post-yield strain rate, which decreases with 
decreasing temperature. However, the values obtained at any one 
temperature show very little variation with strain rate; it is therefore 
considered that, within the range of strain rates and temperatures referred 
to in the table, the effect of variation in strain rate on the dynamic lower 
yield stress is negligible compared to that of change of temperature. The 
ratio of dynamic to static upper yield stress, on the other hand, shows a 
consistent fall with strain rate, at any one temperature ; this is indicative 
of the much greater sensitivity of the upper yield stress to changes in 
strain rate. 

Comparison of figs. 2 and 5 (a) shows that the effect on the stress-strain. 
curve of increasing the average strain rate to 600 sec is similar to that of 
lowering the test temperature to —183°c. In each case the upper yield 
stress is multiplied by about 2-5 and the lower yield strain is considerably 
increased. 

The static stress-strain curves (fig. 2) show that the drop of stress at 
yield and the lower yield strain both increase as the temperature is 
lowered. It seems likely that these effects are related to the fact that the 
mechanism of deformation changes from coarse slip to fine slip as the 
temperature is reduced. This change indicates that during low- 
temperature yielding the number of dislocation sources activated is. 
considerably greater than the number activated at room temperature. 
At any given strain, therefore, the average number of dislocations 
generated per source will be correspondingly smaller. Thus the rate at 
which back stresses build up at the sources with increasing strain will 
be reduced, and consequently the strain which can occur before the onset. 
of work-hardening will be increased. The fact that relatively little- 
multiple slip occurs during low-temperature yielding suggests that 
hardening due to interaction between intersecting dislocations is reduced ; 
this would also tend to increase the lower yield strain. 
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Similar characteristics in the mechanism of deformation have previously 
been observed at high rates of strain and room temperature, and the 
_ present investigation has shown a further refinement of slip when high 
rates of strain are combined with low temperatures. It appears, therefore, 
that the number of sources activated during yielding increases with the 
stress at which yielding takes place. When this stress is large, the rate 
at which energy is being dissipated in the lattice by dislocation movement 
must be large; it is suggested that this energy assists thermal fluctuations 
in freeing dislocation sources which would otherwise have remained 
inactive. 

Figures 2 and 5(a)-(d) show that in those tests in which twins were 
formed (all dynamic tests at —121°c and —84°c, and tests at 503in/sec 
at —41°c) the upper yield stress was greater than about 120000 Ib/in?; 
whereas in the remaining dynamic tests and in all the static tests the value 
of the upper yield stress was less than this. These observations are 
consistent with the hypothesis (Krafft et al. 1953) of a critical stress for 
the formation of twins. The hypothesis has received further support from. 
tensile impact tests (Wood, private communication) in which specimens 
of the same steel as that used in the present investigation were yielded at 
room temperature at stresses up to 160000Ilb/in?; in these tests twins 
were only observed when the yield stress exceeded about 125 000 lb/in?. 

Figures 5 (b)—(d) show that the formation of twins has very little effect 
on the general shape of the stress-strain curve. The fact that the curves 
do not show saw-tooth features, such as have been obtained by some 
investigators, is probably due to the relatively small number of twins 
formed in the present tests. 


§ 6. CONCLUSIONS 


(i) The effect on the compression stress-strain curve of a medium- 
carbon steel of increasing the strain rate to 600 sec“ is comparable to that 
of lowering the test temperature to — 183°C; in each case the upper yield 
stress is multiplied some 2} times and the lower yield strain is considerably 
increased. The effect of testing at high strain rates and low temperatures 
is to raise the curve still further. 

(ii) At a given velocity of impact, the time to yield and the dynamic 
upper and lower yield stresses are increased by lowering the test tempera- 
ture; the ratio of dynamic upper yield stress to static upper yield stress, 
however, is decreased. 

(iii) At low temperatures and normal strain rates the slip occurring 
in the ferrite is finer than that obtained at room temperature; the slip 
lines are straighter and less forked ; also less multiple slip occurs. 

(iv) At low temperatures and high strain rates the slip is still finer, the 
slip lines being almost indiscernible. Twin lamellae are also observed 
when the temperature is low enough and the rate high enough ; the results 
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obtained are consistent with the existence of a critical stress for twinning 
of between 2:7 and 3-0 times the static upper yield stress, for the steel 


tested. 


APPENDIX 


The Effect of Temperature on the Shape of the Elastic Wave and on 
the Wire Resistance Strain Gauge Factor 


It was thought that for a given velocity of impact the elastic stress waves 
at the different test temperatures would differ very little from one another. 
A theoretical estimate of the difference in the maximum stresses at room 
temperature and —121°c may be made as follows: 

Ideally the maximum stress o applied to the high-tensile steel dummy 
by an impact at velocity V is given by 

o=2-10V (Hp)? > g: aa0 i eee. ae eee 


where H ig the value of Young’s modulus for the weighbar steel and p is 
the density of the weighbar. Thus, provided # and p are known at 
+15°c and —121°c, the variation in o may be found. 

Koster’s (1948) expression for the variation of the shear modulus G 
with the absolute temperature 7’, 

G = (13-30 — 0-003’) x 108 Ib/in2, 
indicates that at — 121°c G has increased 3-:1% from the value at + 15°c. 
Hence, assuming that Poisson’s ratio remains constant, H may be 
taken to have increased 3-1%. The increase in p is negligible, taking 
the temperature coefficient of expansion of the weighbar steel as 
1-1x10-° per °c. Thus from eqn. (1) the maximum stress should be 
increased about 14%. 

In deriving this estimate it has been assumed that the entire weighbar 
is at a uniform temperature, whereas in fact only the lower part of the 
weighbar is reduced in temperature by the liquid oxygen. It may be 
shown that the effect of this is to make the increase in the stress applied 
to the specimen even smaller. It is therefore reasonable to assume that 
the shape and magnitude of the incident stress wave is sensibly independent 
of temperature. 

Making this assumption, it is possible to determine the variation in the 
resistance strain gauge factor with temperature. To permit this, elastic 
waves were recorded at temperatures of 15°c, —41°c, —84°c and — 121°c, 
for an impact velocity of 238in/sec. The oscillograms and calibration lines 
obtained at 15°c and — 41°o are shown in fig. 6, Pl. 52, and it can be seen 
that the waves are almost identical in shape. By assuming that their 
magnitudes are equal, and taking the makers’ value of 1-94 for the gauge 
factor at 15°c, the value at —41°o was found to be 1-97. Similarly the 
values at —84°c and — 121°c were found to be 1:93 and 1-98 respectively. 
This variation is of the same order as the accuracy of the measurements 
and hence the gauge factor was taken as sensibly constant over the range 
of temperature investigated. 
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ABSTRACT 


When a strong current is flowing axially in a plasma, with cylindrical 
geometry, the charged particles of the plasma will be prevented from striking 
the cylindrical wall both by the self-magnetic field of the current and, in the 
case of the positive ions, by the electric field due to the negative charge 
in the region of the current. Calculations are made of these effects, which 
are relevant to the loss of heat by conduction from a thermonuclear device 
of the zeTa type. It is shown that in order to reduce the loss of heat by con- 
duction to a value less than the energy generated in the discharge by a thermo- 
nuclear reaction in deuterium, it is necessary either to have a magnetic field 
corresponding to a current of the order of a thousand million amps, or to have 
a substantial electric field, giving a radial potential difference between the axis 
and wall of the order of 250-500 kv, together with currents of a few million 
amps. Reasons are given for supposing that such an electric field would in 
fact be set up as a concommittant of the current and the reaction. 


§ 1. INTRODUCTION 


IF a completely ionized gas in an infinite cylindrical container carries a 
strong current along the axis of the cylinder, the self attraction of the 
moving electrons will produce a negative charge near the axis. This 
negative charge will hold the positive nuclei in towards the axis creating a 
radial gradient in their density. Since however the electrostatic coupling 
between the electrons and nuclei is very large the charge densities of 
electrons and nuclei can only differ infinitesimally, so the density of the 
whole gas diminishes towards the walls of the cylinder. If the gas is hot 
this diminution in density will reduce the loss of heat to the walls by 
conduction and the loss is further reduced by the bending of the paths 
of both nuclei and electrons, but especially the latter, by the magnetic 
field of the current into trochoidal tracks. Apart from a change from 
cylindrical to toroidal geometry, this is the principle on which zeTa and 
SCEPTRE depend for attaining much higher temperatures than could be 
reached if conductivity were normal. 

The object of the present paper is to examine in more detail how the 
conductivity is in fact affected by the two factors of diminished density 
and magnetic field near the wall. For mathematical simplicity cylindrical 
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geometry will be taken and the radial electric field outside the ‘pinch’ 
will be taken as due to a charge on the pinch. 


§ 2 

Spitzer (Physics of Fully Ionized Gases, pp. 42-44) has pointed out that 
the problem of the pinch discharge is not definite unless some assumption 
is made as to how it is set up, but that on a plausible assumption the mass 
velocity of the nuclei—which we shall take as deuterons—is likely to be 
small. In some early unpublished work the present author pointed out 
that, ifa pinch discharge is fed by letting in gas to take the place of particles 
lost to the walls, a quite appreciable axial velocity of the deuterons will be 
set up; but for the present rough calculations this effect will be ignored 
and the mass velocity of the deuterons taken as zero. In that case, assuming 
always that the Boltzmann equation of transfer is valid (which may not 
be strictly true with long range forces but seems the best we can do), the 
equations of radial equilibrium for a steady state become 


dn 
cp Ny inde 
neL,. Pap? 
nev .H = —neH..—kT dn 
oe a ° dr’ 


where n is the particle density of electrons or deuterons, 7’, and 7’, are 
their temperatures, v, is the mean axial velocity of the electrons, H is the 
circumferential magnetic field in gauss, H, is the radial electric field in 
-e.m.u. and is a negative quantity, and e¢ is the charge on a deuteron. 

If T,=T,=T we can rewrite these equations as 


kT d 
f= cape) 


kT d 
v= —H,— me (log 7). 


The first equation leads to 
—eV, (t 
np =n exp ( I i) where V,= — hs E,dr. 

‘This is of course the usual Boltzmann relation which is valid in this case 
in spite of the presence of the magnetic field. The corresponding relation 
for the electrons does not hold because of the magnetic force due to their 
motion, in fact the electrons are Jeast numerous where their electrostatic 
potential energy is greatest. 

Clearly we have v,H=—2E,, showing that the pinch force on the 
moving electrons is twice the radial electric force and balances both the 
-outward electric force and the gas pressure of the electrons. 

If one assumes that v,, produced by an applied axial electric field, is 
inversely as n, the current density becomes independent of n and if the 
axial field is independent of r one finds, after taking account of the Ampere 
relation for H in terms of the current, that the pinch has a sharp edge for 
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a certain radius 7,, outside which there is no gas so that the current density 
falls suddenly to zero. This is obviously at best an approximation. It. 
would mean an infinite potential difference between the pinch and the 
region outside. 

In view of its generality is seems probable that the Boltzmann relation. 
is at least approximately true over the whole volume and assuming its. 
validity we shall proceed to calculate the number of collisions between the 
deuterons and the wall for a given V, and a given current. The calcu- 
lation will be limited to the region outside the pinch which will be assumed 
free of axial magnetic field and V, will be taken of the form 2c log (7/79), 
where r, is the radius of the pinch, which is equivalent to assuming the 
region outside the pinch to be electrically neutral. This no doubt is in- 
accurate, but it seems unlikely that the loss to the walls will depend 
critically on the form of the electric field for fixed J. 


3 

The stages of the calculation are as follows : 

(a) Each deuteron will be moving in a trochoidal path in an overall 
axial direction. Assuming that one of them makes a collision which gives 
it a velocity random in direction but of constant magnitude, we find the 
chance that its new trochoidal motion will make it hit the wall. This 
chance will be zero if it is too far away and will vary with distance if it is 
not. We thus find the number of particles striking the wall per second in 
terms of the density, the electric and magnetic fields and the frequency of 
collision for random scattering. To give an idea of the orders of magni- 
tude in a workable thermonuclear device, the magnetic field from a current 
of 10°amps will give a deuteron whose energy is 50kv a Larmor radius 
about 4 of its distance from the axis. 


(6) The next process is to calculate the frequency of the scattering 
collisions. The most important collisions with Coulomb forces are 
unfortunately the small ones. We evaluate the effect of these by a modi- 
fication of Spitzer’s method (loc. cit. pp. 68-75). In effect we calculate what 
the frequency of large collisions would have to be (if these alone existed) 
to give the same mean-square sideways energy to a particle moving through 
the gas as Chandrasekhar has shown it would actually acquire. This we 
assume to be the effective collision frequency. 


(c) Combining (a) and (6) we find the energy loss to the wall per second 


and compare this with the yield of energy from thermonuclear action in 
the pinch. 


$4 
To calculate the effect of a collision at the point (r, 6, z) on a deuteron. 
of mass m, assume cylindrical symmetry. with a total current J in the 
pinch, where the density of deuterons per cm’ is n). Take the z-axis. 
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along the direction of motion of the electrons giving the current J, the 
equations of motion are 

mé=—eHft; H=2I/r 

m(# — eS) =eHz+ eH, 


whence 


where C is a constant. 
The first equation can be integrated to give 


2=%)— | (eH /m)i dt = 2 — (2eI/m) log (r’/r), 


where 2, is the initial value of 2 and r is the initial radius of the particle’s: 
position. To find the maximum value of r’ we use the equation of energy, 
since at this point *=0, giving 
3m (2? + O7/r'?) = dw? — eV, + mC?/2r?, 
where w is the velocity in the axial plane. 
Hence 
{2 — (2eL/m) log (r’/r)}? = w? — (4ce/m) log (r’/r) + C2(1/r? — 1/r’?). 

We are concerned with the chance that the larger of the two roots of 
this equation in 7’ shall exceed R, the radius of the cylindrical container. 
If it does so the deuteron will hit the wall, since we suppose the free path 
so long that it is unlikely to make another collision before making one 
complete period of its orbit. Clearly the term C? (1/r?—1/r’*) is positive 
and less than the square of the initial circumferential velocity. For ease 
of calculation we shall neglect this term and take w as a constant given 
by 4mw?=kT, i.e. take w? as the mean for two degrees of freedom. The 
error in neglecting O?(1/r? — 1/r’?) is less than the difference between taking 
two and three degrees of freedom for w”. 

Our equation can now be written 

3 = (2Le/m) log (7r'/r) + \/{w? — (4ce/m) log (r’/r)}. meee (i) 
For every value of log (r’/r) less than mw?/4ce there will be two real values 
of 2, satisfying eqn. (1), but these do not always represent possible paths 
since |2)|<@. Within the range —w<z2)<@ the chance of a deuteron 
having a value of 2, between 2, and 2, + 52, is 62)/2w, provided the direction 
of the initial velocity of the deuteron is random. 


$5 
We shall now calculate the chance that a deuteron which makes a 
collision at (r, 6, z) shall hit the wall. This is the chance that 2) satisfies 
eqn. (1) with r’>R. We calculate this probability as a function of 7. 
We write eqn. (1) as follows: z=az + 1/(w?—bx) where z= %, x= log (r'/r), 
a=2lelm. b=4ce/m. Clearly x<w?/b. We write log (R/r) = &, 
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b/w? =2ce/kT'= x. If zis plotted as a function of x we obtain a parabola. 
Three cases arise : 

Case 1. Thelinez=w cuts the curve only once fora>0. The condition 
for this is 2aw <b, which means that the magnetic field is comparatively 
weak. In this case, if z,, 2. are the two values of z corresponding to 
By %—%_ = 24/(w?—b&) and any value of z between z, and 2, corresponds 
to a value of «>, i.e. to a deuteron hitting the wall. Hence the chance 
of a particle hitting the wall is 


24/(w?— be) _ Mba tk 
Ba - fi Sloss} 


provided this quantity is real, which is the case ifr >7, where 6 log (R/7)=w?. 
Now by Boltzmann’s law 
n=n, exp { — 2ce log (7/79) }/kT = n(7/79)~* 
where a= 2ce/kT’. 


Fig. 1 


‘The total number hitting the wall is 
-R ——s 
(C n® ow s/{1—alog (R/r)} 2ar dr. 


After each collision either or both the particles may hit the wall at some 
stage of the trochoid, since in general a particle will describe many loops 
before making a second collision. Here a,@ denotes the mean over the 
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scattering cross section and relative velocity. Hence the number is 


el R 
2 . Og). My? . N** r—*e+14/{1 —alog (R/r)} 


= 2mro wn? ,(R2/x) { exp {t(2 —2/a)}. 4/(1—#) dt 


where Np =Ny(R/r9)-* and t=alog (R/r) 

the integral though a function of « is of order unity for a>1. It is 1:19 
for a=3. Hence the number striking the wall is approximately 
270 ,0N 2 R2kT |2ce which does not contain I. 


Case 2. 2aw>b but aw<b. The parabola touches x =w?2/b at C whose 
ordinate is less than w. The line z=w. cuts the curve at D. Here 
again the number hitting the wall depends only within rather narrow 
limits on the parameters. The numerical values are not much different 
from Case 1 and we shall therefore use whenever aw <b, the expression 


Qro wn RT /2ce for the number striking the wall per cm per sec. 


Case 3. Strong magnetic field aw >b. A calculation on the same lines as 
‘before gives the number per em per sec reaching the wall as approximately 


2ro wn, RL —y+y2(exp yt — 1)]2w/a where y= a/4e(«—1). For aw>b 
this becomes 270,w n,2R2w/a. This expression varies inversely as It 


$6 
It remains to determine c,w. In doing so we shall assume the 
deuterons to be moving in epicyclic orbits parallel to the walls, and we 
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shall make the usual assumptions about distant collisions (which with 
long-range forces are apt to lead to an infinity), namely that collisions in 
which the shortest distance of approach is more than a Debye radius do 
not count because the forces are neutralized by other charges. This, as 
is well known, leads to multiplying the cross sections for ‘large’ collisions 
by a logarithmic factor which has been tabulated by Spitzer and is of the 
order 20 in the cases of interest. This gives a cross section o, ~ 6 x 10~19/W? 
where W is the relative kinetic energy of the colliding particles. If o, is 
averaged over the relative kinetic energy of the particles in the simplest 
way the result contains another logarithmic infinity. I have therefore. 


gone into the matter in somewhat greater detail and the results shows that 
the logarithm occurs only as a simple factor. However, it is better to 
use the results of Chandrasekar as quoted by Spitzer (loc. cit. p. 73). These 
give an expression for the mean value of the square of the velocity at right 
angles to its original direction acquired per second by a test particle in 
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moving through the gas, written as ((Aw ,)*). If this quantity is multi- 
plied by the number of particles per em’, each of which can be regarded as 
a test particle, and by two, to take account of the particles that do the 
scattering, one gets the sum of the squares of all the velocities perpendi- 
cular each to its original direction acquired per second by all the particles 
in unit volume of the gas. We call this quantity S. 

The scattering is in fact caused mostly by a great number of small 
collisions, but we shall not go far wrong if we replace this by that number 
of ‘large’ collisions, involving a sideways velocity of «/(kT'/m) which 
would give the same value of 8S. Ifo, is the cross section for this process, 
g,w (the mean of o, multiplied by the relative velocity w) is given by 

S=n2o,w kT /m. 
But S=2n<(Aw , )?) by definition. Hence 
_of m\ ((Ao,)*) — 2 8aetlog A / m \32 

rao=2( 7) ee (mm) {O(w) —G(yo)} 
in Spitzer’s notation, writing 4/(3k7'/m), the r.m.s. velocity, for Spitzer’s 
‘w’, the velocity of the test particles. From his table 5-2 

(1,0) — G(L,w) = 0-629 

for the case when the test particle and those of the gas have equal velocities. 
Hence 


_ 18-2etlog A 


OW 
m2 


s 


m 3/2 
(ey =0-34log A 1-32 
(i) Be ee : 


with 7 in degrees. 


$7 
The same argument should apply to the electrons for which the particle 
density is the same. In fact o, is the same for deuteron—deuteron as for 


electron—electron collisions, but o,w is larger for the latter for equal 
temperature in the ratio of the square roots of the masses. However 
w/a is smaller in the same ratio so o,w.w/a is the same for both. Thus 
for very strong magnetic fields the first and most important term in 
the expression is the same for both for equal temperatures. The later 
terms give more deuterons sent to the wall than electrons. On the other 
hand the scattering of deuterons by electrons is negligible, but that of 
electrons by deuterons is substantial. Even in a strong magnetic field 
there will be a slight tendency to form a sheath near the wall due to excess 
loss of electrons. In a weaker field (cases 1 and 2) this tendency will be 
even more marked. 


§8 
It is of interest to compare the loss of energy by conduction to the walls 
with the energy produced by nuclear action and to calculate the condition 
that the former shall be a moderate fraction f of the former. In an actual 
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D-D power producer the temperature must exceed 35kv for the energy 
produced to exceed the radiation not counting neutrons, and so keep the 
reactor hot. We will calculate results for 50kv and 100kv. 
We double the expressions found for the loss of deuterons to allow for 
a loss of electrons equal to that of deuterons; each particle will on average 
carry energy 2k7'. We assume constant density m) in the pinch where 
alone the nuclear reaction is important and denote by o,,w the mean product 
of cross section and relative velocity for this process; then 
a0 49x 10-33 
Fae belo 
? -18 
At lOOky Ce! ere 
OW 3:0 x 10 


The following expressions hold for the different cases : 
Cases 1 and 2, 


At 50kv =3°5 x 104 


= 5:8 x 10%. 


aw <b, i.e. [Tw < 20 


fas zen 2 R22kT kT 


where £ is the energy released and retained in the vessel per D—D reaction 
including secondary D-T reaction; H=4-16 Mev. 
Case 3. Iw>2c 


Play Gist JRE 
—— Se om _y tyX(expyt— 1)]. 


f=8 


Potentials in kv 


Temp (kv) Cases 1, 2 Case 3 


Current (amps) 3 x 10° | 5x 108 10? 12% 10° 


50 257 233 210 193 — 
100 470 = 409 343 295 


Cases 1 and 2 apply up to a limit of 3-2 108 ; : 
amps for 100 ky. a amps for 50kv and 4:2x10 


We now examine which of cases 1-3 is most likely to apply to some 
special examples. The deciding parameter is b/aw=2c/Iw, where 
20 log B/rg is the potential difference between pinch and wall. 

For zeta: I~ 2x 104 (200000 amps) 

w ~ 6 x 10° corresponding to 5 million degrees k 
o has not been measured but log (R/r,) ~ 1. 


The p.d. might be 100-1000 volts, ie. c=10!°~ 10"; thus 
10 


1 
Odor tO 
PS rs 
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For a possible power producer with D-T': 
I~2x 10°, i.e. 2 million amps 
w~2x 10? for 50 million degrees k 
2c log (R/ry) of order 5kT' ~ 25kv 
log (R/r9) ~ 1 giving b/aw = 25 x 10/4 x 10!2~ 1. 
Using the Bennett relation J?=2Nk7’,, this would give 
(2 x 105)? 
2x 1:36 x 10-16 x 6 x 10’ 
or with 100cm? area of pinch, 1-5 x 10!° particles of either sign per cm3,. 
a reasonable figure. 
For a reactor aiming at energy from D-D: 
I~5x 10° 
w~8x 107 for a temperature of 50kv 
2c log (R/79) ~ 200 kv 
blaw = 1/2. 
For this the Bennett relation gives 
N=1-6 x 108 perem. 

We conclude that we may expect conditions somewhere near the border-. 
line between cases 2 and 3. Thus the containment is due about as much 
to the radial electric field as to the magnetic field; the two help one 
another. The electric field which produces the radial change in density 
can only exist in a good conductor like a plasma in the presence of the 
magnetic field. A magnetic field is possible alone, if it can be established 
before the plasma conducts too well, but a very large field is required to 


produce good containment. It is worth examining our results in the 
limiting case when the electric field vanishes. 


=3 x 1018 percm, 


§9 
Case of zero electric field, i. b=0, ~=0, n is constant. We now have 
for the number of deuterons reaching wall per cm per second 


me a Rk a R 
5, cated AY ae cae —\)=1 
3 (2 = log =) r dr where on log (=) 


= na, ngtP| 1+ = {exp(— ~)- 1} | : 
163) 


The ratio of the loss by collision at the walls to the energy produced is, 
allowing for electrons. 


ow R22kT a — 4y 
Eat i fae i ag id 
ate E [ +e {e( a ) a 


which gives at 50 kv 


a ( —4w 
2:5 x 108 | 14 as (Ca (=) = 1} | : 


Qo gw . Ny" | 
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‘To make this small a/4@ must be large; the expression is then approxi- 


mately equal to 


4-2 x 108 2x 108 


at 50 kv and at 100 kv. 


Hence we should need more than 109 amps to reduce the loss by con- 
duction to no more than the nuclear energy produced, if no electric field 
were present. 


§ 10 

We have not considered how the radial field will be produced, or whether 
in fact it will be sufficient. A state of zero field, though unlikely, could 
conceivably arise if outside the pinch the electrons and deuterons had equal 
but opposite axial velocities. However, the fact that the simplest relation 
between particle density and current density leads to an infinite potential 
difference suggests that the field will be adequate. The central negative 
charge is produced and sustained against loss of electrons to the walls by 
the action of the axial electric field which, accelerating the electrons in the 
outer parts of the discharge, drives them in towards the centre by inter- 
action with the field H. Because of their greater mass the deuterons are 
-accelerated less. 

The electric field is a link in the chain by which the deuterons are 
anchored, through the electrons to the magnetic field. 

Once the reaction has started the fact that the high energy products 
have positive charges is an independent reason for a negative charge to 
build up on the axis. 


§ 11. ConcLUSION AND SUMMARY 


When completely ionized deuterium is contained in an infinite cylinder 
by a strong electric current down the axis there will usually be a radial 
electric field set up. We have calculated the effects of this electric field 
and of the magnetic field of the current acting together in diminishing the 
loss of thermal energy from particles hitting the wall. In the absence of 
any electric field this loss is very large and in this case it would require a 
current of several thousand million amps to reduce the energy loss 
below the rate of production of thermonuclear energy even at temperatures 
of 100kv (about 1000 million degrees). However, a radial potential 
difference of 250-500kv between the pinched current and the wall will 
reduce the energy loss to acceptable values even with currents of the order 
of a few million amps. Reasons are given for supposing that these 
potential differences are likely to arise automatically. 
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ABSTRACT 


Electron transmission micrographs of polycrystalline aluminium quenched 
from ~ 600°C into iced brine show that the specimens contain many disloca- 
tion loops with diameters of the order of several hundred A. Some of the 
loops are hexagonal in shape, and occur on (111) planes, and have a Burgers 
vector } [110] not in the plane of the loop. The loops are thought to be 
produced by the collapse of dises of vacancies, and the corresponding vacancy 
concentration is 10-4. Some of the loops disappear, while under observation 
in the microscope, by a mechanism of prismatic slip. In the neighbourhood 
of grain boundaries there is a zone depleted of loops; similar zones are 
observed near dislocations, whose complex configurations suggest that climb 
may have taken place. The observations are discussed in relation to the 
problems of ‘ annealing out’ of vacancies, of quench-hardening, and of the 
origin of the dislocation network. 


§ 1. INTRODUCTION 


WHEN a pure metal is rapidly cooled from a high temperature, vacancies in 
concentrations of up to 1 to 104 can be quenched in. These vacancies 
produce an increase in electrical resistivity and measurements of this 
physical property can be used to determine the concentration of vacancies 
and to follow the ‘annealing out’ of the vacancies during ageing at suitable 
temperatures (for review see Cottrell 1957, Broom and Ham 1957). For 
aluminium ageing effects are observed at room temperature (Wintenberger 
1956, 1957, Bradshaw and Pearson 1957). Apart from the resistivity 
effects, quenching also leads to a decrease of the damping dueto dislocations, 
as shown by internal friction experiments (Levy and Metzger 1955), and to 
large changes in the mechanical properties of metals. The shape of the 
stress-strain curve is very dependent on the rate of cooling and large 
increases in the yield stress may be introduced by quenching (Maddin and 
Cottrell 1955). This increase in yield stress is associated with a marked 
change in the appearance of the slip lines (Maddin and Cottrell 1955), an 
increased tendency to ‘overshoot’ (Kimura and Maddin, reported by 
Cottrell 1957), and a different preferred orientation texture (Smallman and 
Westmacott, unpublished) on subsequent deformation. (For review see 
Cottrell 1957.) 


+ Communicated by the Authors. 
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While much valuable information has been gained about the energies of 
formation and movement of vacancies for example from electrical resistivity 
effects, it is not known by what mechanism the vacancies are annealed out, 
nor has the cause of hardening been established. However, a number of 
theories have been suggested, which propose the formation of either cavities 
(Coulomb and Friedel 1957), collapsed discs of vacancies (Kimura, Maddin 
and Kuhlmann-Wilsdorf, unpublished) or jogs at dislocation lines leading 
to climb (Cottrell 1957). 

Recently two of us (R. E. 8S. and K. H. W.) have studied quenching 
effects by means of a small angle scattering technique. Specimens of 
polycrystalline aluminium were quenched from 630° into a calcium chloride 
solution at —45°o; the specimens were transferred immediately to the low 
angle scattering spectrometer (Smallman and Westmacott 1958) and the 
intensity of scattering measured in the region of 20= 4° (26=the angle- 
between the incident monochromatic ,copper Ka beam and the scattered 
beam). The intensity of scattering was found to increase with time till, 
after a few minutes, a constant value was reached. The angular breadth ~ 
of the scattering curve is about 4° and the intensities are very large. These 
intensities are too large to be explained in terms of true low angle scattering 
from a reasonable density of dislocations, but could be explained either in 
terms of small cavities, about 100° A diameter, or by the mechanism of 
double-Bragg scattering (Atkinson 1958) from dislocation loops, 100-500 4 
diameter. However, experiments by Takamura (1956) and by Koehler 
and his co-workers (see Cottrell 1957) indicate that density changes anneal 
at the same rate as the resistivity and this suggests strongly that the. 
vacancies are disappearing during the ageing process and that they are not. 
forming small cavities. 

Electron transmission microscope experiments were therefore carried 
out in an attempt to observe the dislocation loops directly, if they existed. 


§ 2, EXPERIMENTAL 


Electron transmission microscopy was carried out both at the Cavendish 
Laboratory, using a Siemens Elmiskop I operating at 80 kv, and at. 
A.E.R.E., using a Philips (100 kv model) microscope. Dislocation loops 
(see §3) were observed with both instruments, but the higher beam 
intensities, and the wider facilities for selected area diffraction experiments 
and for tilting the specimen enabled more detailed information to be 
obtained with the former. 

Specimens of polycrystalline aluminium, 99-995°%, pure, 0-075 mm thick, 
were quenched from a vertical tube furnace at ~ 600°C into iced brine. 
The foils were then thinned by electropolishing in a bath of 20% perchloric 
acid, 80% absolute alcohol (Murphy 1958) and examined in the microscope. 
For comparison purposes a furnace cooled sample was alsoexamined. The 
microscope specimens were tilted by means of a stereoholder to bring any 
dislocations into contrast. It is essential to follow this procedure to 
obtain a representative picture of the distribution of defects in the specimen. 
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§ 3. REesuLTS 


3.1. Contrast Effects 


Figures 1-3, Pls. 53-55 and figs. 6-9, Pls. 56-59 are micrographs of repre- 
sentative areas of the quenched specimens. All the photographs show 
contrast effects in the form of loops or of double lines. Experiments on the 
effect of tilting the specimen prove conclusively that the contrast arises 
from the Bragg diffraction mechanism (see for example Hirsch ef al. 1956). 
Two possible explanations have been considered. The contrast could be 
caused by suitably shaped cavities, the loops effectively corresponding to 
thickness extinction contours (Heidenreich 1949). However, if the 
cavities were discs, in general strong uniform contrast would be observed 
inside the extinction contour ring; if the cavities are spherical or if they 
have some other rounded shape the contours would shrink or expand on 
tilting the specimen. None of these features are observed; in fact the 
variation of the contrast as the specimen is tilted is typical of that observed 
from dislocations (Hirsch et al. 1956). Further, the loops have been 
observed to move by glide in a manner expected from dislocations (see 
§ 3.5). It is concluded therefore that the contrast is due to dislocation 
loops. When these loops are lying in a plane which is nearly parallel to 
the plane of the foil the contrast is expected to be uniform around the loop. 
However, when the loop lies on an inclined plane, the contrast will not be 
uniform. The kinematical theory of diffraction contrast from dislocations 
(Hirsch and Whelan, unpublished) shows that dislocations become narrower 
and less visible the steeper their inclination to the plane of the foil. The 
double lines frequently observed therefore presumably correspond to 
dislocation loops on inclined slip planes, for which the steeply inclined 
parts of the loop are invisible. 


3.2. Density and Sizes of Loops 


Figure 1 shows a representative area of one of the specimens (to be 
referred to as specimen X). The distribution of loops is not uniform, and 
the diameters of different loops are not all equal; the average number 
of loops is about 1015 em~3, and the average diameter ~ 200A. Assuming 
the loop to be a collapsed disc of vacancies, the corresponding initial 
concentration of vacancies is found to be ~10-*. This value agrees in 
order of magnitude with results reported from resistivity measurements 
(Bradshaw and Pearson 1957) and confirms the value of energy of formation 
of vacancy deduced from them. The results also show that the density of 
dislocations in loops is about 101° cm~?. 

Figure 2 shows a high magnification micrograph of another area of the 
same specimen. The wide distribution of sizes is evident; the smallest 
loops that can be definitely recognized as dislocation loops have a diameter 
of ~ 1004. 3 

In another specimen, (Y), quenched under apparently similar conditions, 
the loops are considerably larger, about 1000 A in diameter (figs. 3 (a, b, ¢)). 
The average number of loops is about 3 x 10! em~®, and the corresponding 
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initial vacancy concentration again about 1 in 10%. The density of dis- 
locations in loops is about 10° em~?. 

Apart from the loops a much smaller number of other dislocation lines, 
most of them in complicated configurations, are observed (see § 3.4). 

For comparison purposes a similar aluminium specimen was furnace 
cooled from 600°c over a period of three hours. Figure 4, Pl. 55 is a micro- 
graph from a typical area. The sharp bands are extinction contours ; 
by tilting the specimen these contours were swept across the field of view so 
that any dislocations present could be made visible. Most of the areas 
observed contain no dislocation loops at all; in one area a few ‘loops’ 
were observed, but these could not be identified with certainty as small 
dislocation loops. The density of small loops, if there are any at all, is 
certainly less than 10!2 cm~*. However some isolated dislocations running 
through the foil from one surface to the other, have been observed ; some 
can be seen on fig. 4 at A’, where they are shown up by kinks in the extinction 
contours. Although no exact figure can be given, there are about 108 
dislocations per cm? in this specimen. 


3.3. Nature of the Loops 


Many of the loops appear to have regular crystallographic forms. They 
often have the shape of parallelograms or of regular hexagons (A, figs. 
eo. 376800), 

In figs, 3 (a, b) a particularly well developed hexagon can be seen at A. 
The orientation of the plane of the foil, determined by selected area electron 
diffraction, is several degrees away from a (111) plane, but the hexagon 
lies very probably on this (111) plane, and its edges are parallel (at least to 
within a few degrees, the uncertainty in the orientation determination) 
to the [110] directions in this plane. The contrast from the loops B, and 
B, suggests that these are lying in or near two of the other (111) planes. 
Figure 5 (a) shows the hexagon A and the reference tetrahedron (Thompson 
1953) giving the orientation of the other (111) planes. The planes corres- 
ponding to loops B, and B, are labelled c and b respectively. The geometry 
of these loops is strong evidence that at least some of them correspond to 
collapsed dises on the (111) planes, although at present it is not possible to 
say whether loops are formed also on other planes. 

It should be emphasized that these loops are not Frank sessile dislocations 
(Frank 1949) as no stacking fault contrast is observed (Whelan et al. 1957). 
However, as Kuhlmann-Wilsdorf (1958) has pointed out, in metals of high 
stacking fault energy it is energetically favourable to nucleate a Shockley 
partial in the disc of stacking fault associated with the Frank sessile pro- 
duced by collapse of a dise of vacancies on the (111) plane. Following 
Frank’s calculation (Frank 1950) on the stress required to nucleate a disloca- 
tion loop at room temperature, and remembering that a stacking fault 
exerts an effective stress of y/b, where y is the stacking fault energy, and b 
the Burgers vector of the Shockley partial, the stacking fault will be removed 
if y > Gb/25, where G is the shear modulus. For aluminium this condition 
is satisfied if y > 180 ergs cm?; actually for aluminium, y is expected to be 
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at least as great as this value, so that it is clear that the faults will not be 
stable. The dislocation ring will therefore be a whole dislocation with 
Burgers vector $ [110] not lying in the (111) plane of the disc, the reaction 
being $ [111]+ 4 [112]=4 [110]. The dislocation loop is therefore a pris- 
matic dislocation (Read 1953) lying on the surface of a cylinder whose 


(a) The hexagonal loop A of figs. 3 (a, 6) and the reference tetrahedron giving 
the orientation of the other (111) planes. The heavy line of intersection 
between planes a and b is the Burgers vector of loops A and C (fig. 3). 


cross section is determined by the dislocation loop and whose axis is parallel 
to [110]. If the loop is in the form of a hexagon with sides parallel to 
[110] directions in the (111) plane (as in the case of A, figs. 3 (a, b)) it can 
glide on a cylinder bounded by pairs of (111), (111) and (001) faces (fig. 5 (0)). 
The dislocations can glide easily on the (111) type planes, but probably 
not so easily on the (001) planes, as slip has not been observed on (001) 
planes at room temperature. If the dislocations in the (001) plane are 
anchored, the other dislocations could act as Frank—Read sourcesy. 

Apart from the hexagonal, parallelogram, and many apparently irregular 
shapes, occasionally double loops (fig. 6 (a)) and possibly spirals W (figs. 
6 (b, c)) are observed. These are presumably produced by climb and are 
similar in nature to (but rather more irregular than) those observed by 
Amelinckx in NaCl (Amelinckx 1957). 


+The nature and properties of these dislocation loops are precisely those 
predicted and discussed by KuhImann-Wilsdorf (1958), who refers to them as 
R-dislocations. 
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3.4. Distribution of Loops 
The distribution of loops is not uniform throughout the specimen. In 
particular near grain boundaries there is a denuded zone extending to 
about one micron from either side of the boundary (fig. 7). This is to be 


Fig. 5 (0) 


(6) Geometry of movement of dislocation loop A. The hexagon A corresponds 
to the original position of A, lying in a (111) plane d. The segments 
of this loop lie in (111) planes a, b and the (100) plane e, respectively. 
These segments move along the hexagonal cylinder (whose axis and 
edges are parallel to the Burgers vector) so that A changes to the 
configuration LMKJ. The two segments LM and JK break through 
the surface forming the curved ends of the slip traces seen on fig. 3 (c) 
and copied on fig. 5 (c). The dislocations JL and KM move in planes b 
to produce the slip traces JR, LT and KN, MQ respectively. In this 


drawing planes u and | are respectively the upper and lower surfaces 
of the foil. 
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expected as the grain boundary is essentially a sink for vacancies so that 
in the region of the grain boundary the degree of supersaturation is always 
‘too low to cause the formation of loops. 


Fig. 5 (c) 
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(c) Slip traces produced by loops A and C. A gives rise to JR, LT and Be : 
MQ; C gives rise to FS, HU and GO, IP. Loops D and E are marked 


as reference points (ef. fig. 3 (c)). 


Figure 8 shows a denuded zone in the region of a complicated dislocation 
line. It is likely that here a dislocation line was present in the specimen 
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before quenching. It is suggested that the denuded zone around this 
dislocation line is evidence that a dislocation is a sink for vacancies, and 
that the complex configuration of the dislocation is the result of climb of 
various sections. However, some of the complexities may be due to 
coalescence of the dislocation line with loops such as at X on fig. 9. There 
is little doubt that the dislocation twists from one plane to another in an 
irregular way. These dislocations are quite different from the relatively 
straight dislocations introduced by deformation (e.g. Whelan et al. 1957). 


Fig. 10 


Junction by glide of two prismatic dislocation loops lying in parallel planes. 
(taken from Amelinckx 1957). 


Figure 9 shows a case where very large loops are present; the density 
of small loops is again relatively small in this area, particularly inside 
the large loops. Different segments of these loops are probably not all 
on the same planes. It is thought that these large loops arise as a result 
of the junction, by slip, of several loops not lying in the same plane as 
illustrated in fig. 10 taken from Amelinckx (1957). This reaction is 
expected to take place whenever the diameter of the loops is greater than 
the distance between them. The loops Y on fig. 9 have shapes consistent 
with that shown schematically on fig. 10. 

Complicated zig-zag dislocations can also be seen on fig. 9 at Z. 


3.5. Dislocation Movement 


Some of the loops have been observed to move in the electron microscope 
Figures 3 (a, b,c) are a sequence of pictures showing this moverionte 
comparing figs. 3 (b) and 3 (c) it appears that loops A and C have moved 
away. rom the slip traces which have appeared on fig. 3 (c) it is possible 
to deduce the nature of the slip which has taken place. The details are as 
follows: fig. 3 (a) shows a slip trace, one edge of which passes through loo ; 
D. This trace has disappeared almost completely on fig. 3 (6). Bares 
figs 3 (b) and 3 (c) loops A and C have disappeared and certain new slip traces 
have appeared. Because of difficulties of reproduction these slip traces 
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have been drawn schematically on fig. 5 (c). In order to avoid confusion 
only some of the points have been marked on fig. 3 (c). The striking feature 
of the slip traces is the presence of the curved segments FG, HI and LM. 
A very faint curved segment can be seen on the original at JK. The 
interpretation is as follows. Under the action of a shear stress loop A has 
slipped on its cylinder (see fig. 5 (b)) until opposite sides of the loop meet 
the surface (fig. 5 (b)). This is a case of prismatic slip (Read 1953) similar 
to that observed recently by Jones and Mitchell (1958) in silver chloride. 
For loop A the curved slip traces LM and JK thus correspond to the trace 
produced when the segments lying in planes a, e reach the surface. The 
directions JL, KM thus correspond to the directions of the dislocations 
in plane b when the dislocations in the other planes reach the surface. 
The direction of the line joining the corner X of the hexagon on figs. 3 (a, b) 
to M should correspond to the projection of the axis of the cylinder, i.e. 
of the Burgers vector of the prismatic dislocation, on to the plane of the 
foil. Although no accurate measurements are possible, this direction is 
consistent with the projection of a [110] direction not lying on the (111) 
plane of the dislocation loop on to the plane of the foil. Figure 5 (a) shows 
the hexagon and the tetrahedron showing the other (111) planes and [110] 
directions. The direction of the Burgers vector of the prismatic loop is 
indicated. This Burgers vector is just of the type expected for such a 
loop (see § 3.3). The dislocation JL then glides on plane b to produce the 
slip trace JR, LT, the dislocation KM is of opposite sign and moves in the 
opposite direction to produce the slip trace KN, MQ (see figs. 5 (b) and 
5(c)). For loop C the curved traces FG, HI can be seen more clearly. 
In this case the Burgers vector is found to be parallel to that corresponding 
to loop A; it must also be of the same sign as the dislocations move in the 
same sense as for loop A. However the shape of the loop C suggests that 
it lies on plane c rather than on the plane d of loop A. Figure 5 (a) shows 
that the Burgers vector of the loops could in fact be associated with pris- 
matic loops either on planes c or d, but not with those on planesa,b. For 
loop c the dislocations FH and GI finally glide away on planes b giving 
rise to the traces FS, HU and GO, IP respectively. The edges KN, GO 
and IP, MQ from traces from the two dislocation loops overlap along 
GO and IP and cannot be resolved individually. However the presence 
of GO and IP can be inferred from the enhanced contrast at G and I. 

It follows from above that although one segment of the hexagon lies 
on a (100) plane (plane e in fig. 5 (b)), this dislocation must be able to. 
glide on this slip plane, since it has evidently slipped out of the crystal. 
This is therefore direct evidence for slip on a (100) plane, but it must be 
remembered that the stresses are high in the foils (probably Zz G/1000), 
(see Whelan ef al. 1957) and the temperature may perhaps be as much as 
100°c. Slip on (100) planes in aluminium has been reported previously 
at elevated temperatures (see Maddin and Chen 1954). It can also be 
seen on figs. 3 (b) and 3 (c) that the dislocation V has moved slightly during; 
the interval between the times of the two exposures. 
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§ 4, DiscussION 


The results leave little doubt that the dislocation loops are produced 
by the collapse of discs of vacancies. The geometry of the dislocation 
loops, their Burgers vector, the nature of the slip produced by the loops, 
and the concentration of vacancies corresponding to the collapsed discs ~ 
all support this hypothesis. In addition the denuded zones at grain 
boundaries and at dislocation lines suggest that these act as sinks of 
vacancies. The complicated configurations of the dislocations (other than 
loops) are probably at least in part due to climb. 

Since dislocations (other than loops) and grain boundaries occur 
comparatively rarely in these specimens, most of the vacancies quenched 
in at room temperature form discs which collapse into loops. The annealing 
of the resistivity after quenching is therefore due to the disappearance 
of the vacancies by this mechanism as suggested for Au by Kimura, Maddin 
and Kuhimann-Wilsdorf (unpublished). 

With regard to the quench-hardening effect there is little doubt that 
the loops will cause considerable hardening, either by interaction with 
other dislocations passing through the crystals or possibly because they 
might act as very short length Frank—Read sources (Kuhlmann-Wilsdorf 
1958). The experiments also show that the other dislocations present 
have very complicated configurations probably owing to climb, and also 
by coalescence with loops, and their glissile sections probably correspond 
to some of the straight section seen for example on fig. 8, which extend 
only over a few 10004 at the most. If these sections act as sources the 
flow stress would be several times greater than that for a slowly cooled 
specimen. Such an increase in stress is of the same order of magnitude 
as that observed (Maddin and Cottrell 1955). It should be emphasized 
that as far as one can tell from the photographs the complex dislocations 
are made up of segments of lengths varying from 1004 to several 1000 A, 
lying in different slip planes. It is quite possible that the straight segments 
themselves may contain a number of single jogs; such a ‘joggy’ disloca- 
tion may be hard to move (Cottrell 1957); it is thought however that on 
application of a stress the jogs would move conservatively sideways, 
thereby increasing the free lengths of dislocations to a distance of the order 
of the distance between nodes. At present there is not sufficient evidence 
to decide which of the various hardening mechanisms discussed above is 
the most important. Furthermore, the hardening effect after quenching 
a depend on the dislocation configurations present at the high tempera- 

ure. 

The question arises whether these loops correspond to the dislocation 
sources in crystals, as suggested by Kuhlmann-Wilsdorf (1958). This 
is a very attractive hypothesis. It may be thought that during recrystal- 
ization at high temperatures all dislocations could be removed. Certainly 
more than 99% of all dislocations in the cold-worked aluminium are 
annealed out. If all dislocations are annealed out then new sources would 
be produced by the collapse of vacancy discs after cooling. The slower 
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the rate of cooling, the smaller will be the degree of supersaturation and 
hence the fewer and larger will be the stable loops. In the second specimen 
discussed in this paper the dislocation loops observed are already of the 
order of 10004 diameter. If for slower rates of cooling the diameters of 
the discs are larger and extend to about one micron, and if there are 
sufficient loops, these might grow and intersect. Then, since 2 of the dis- 
locations are glissile, in general the 4-dislocation nodes at the intersections 
will degenerate into two 3-dislocation nodes, thus forming the Frank— 
Read net-work. 

Finally it should be pointed out that the smallest loops which can be 
recognized as dislocation loops have diameters of the order of 1004. It 
is possible that smaller loops or possibly smaller cavities exist, but the 
concentration of vacancies corresponding to the observed loopsis sufficiently 
high to account for at least an appreciable fraction of the quenched-in 
defects. 

Experiments are now in progress on other metals, and on the influence 
of different quenching rates and temperatures, of annealing treatments 
and of deformation on the type and distribution of the loops formed. 
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A Generalization of MacDonald’s Inversion Theorem 


By F. N. H. Rosinsont+ 
Clarendon Laboratory, Oxford 


[Received May 29, 1958] 


MacDonaup (1949) has given an inversion theorem which allows one to 
express the power spectrum of w(f) of any fluctuating quantity X(¢) in 
terms of the ensemble average (X,”) where X, is the average of X over a 
finite time r. This theorem is especially useful in shot noise calculations 
(X,*) is usually directly calculable from simple physical and statistical 
arguments. The theorem states that if y(7) is defined by 


Wa=aCx 2) Sr ewes aly PR ECL ) 
then the power spectrum is: 


w(f)= anf | ° SE sin 2afrdr, tah eae (3) 


In this paper the theorem is extended to cover the case when we have 
two fluctuating quantities X(é) and Y(t) and w(f) is the power spectrum of 
X(t) Y(t). 

We take an interval of time 0—7' where we have in mind the eventual 
passage to the limit 7’ oo; then during this period we can expand X(é) 
and Y(t), the results of one particular measurement, as Fourier series of 
fundamental period 1/7. Adopting a unit of time 7'/27 this gives: 


X(t)= S a CORLL Ona umts i gute: «0 g(a) 
Y()= 5 eCOs (at.)s him ee test cc (30) 
n=0 


The power spectrum of X(t) Y(t) is then defined by 
w(f)=4 lim Tx,,y,, 008 ($y, — Fn) te. ae in (4) 
T>o 
where f=n/7'. Thus w(f)df is the time average of X(t) Y(t) due to com- 


ponents of frequencies within df. We note also that the interval df 
between terms in (3a) is 1/7 so that if we replace a sum by a continuous 


integral we have df=1/T, ie. 
[OP EH=EL enn C08 (bn On ss 8) 


+ English Electric Fellow. 


910 Correspondence 


We now take a period ¢ to t+7 in 7 and form the averages of X(t) and Y(t) 
over this period : 


Sinem ok = 5% [sin {n(¢+7) +$n}—sin (nt +$,,)], 


=-y% [sin {n(t +7) + 6,,}—sin (nt + 6,,)]. 
Rearranging the sine terms 


== <1 cos {n(t-+ 47) Yr ea EE, oo oe a 


ae ripe ren ea ie BR 
The 1 


aye > ann nul cos {m(t+ 47) +¢} cos {n(t + $7) + 6,,} sin 3n7 sin 37. 
pe mn 
(7) 
We now average this expression over all intervals t to t+7, i.e. over all 
t within 7. Neglecting the overlap at the end of 7’, i.e. assuming 7<T 
we have 


XY tee a y: st cos (¢,, — 6,,) sin? nr. atte ee 
ap 
We now define (7) by 
wr) =r X Yt £) 0 08S oa ee et 
Then 
= => =a cos(¢,—-@,)sinnz. «7, 0 nO) 


Replacing the expression on the right by an integral over the power 
spectrum w(f) defined by (4) we have 


oF 2 | Sop sin Daf df. Se) Se ee 
It follows from (11) that 05/d7 is an odd function of 7, in agreement with 

(8) and (9) which give 7 as an even function of 7 so that assuming the usual 

convergence requirements to be satisfied the inversion of (11) is 


w(f)= def | ot sin 2arfT dt. <. 2) le 

When X(t)= Y(t) this result is identical with MacDonald’s original 
theorem. An application of this theorem to the problem of correlations 
between current and velocity fluctuations in electron flow will be given in 
a forthcoming paper. 

When the theorem is applied to some limiting physical situations it 
may appear that the integral arising in (12) is not well defined from an 
elementary point of view. It therefore seems worthwhile to discuss this 
point here. | 
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If X(t) = Y(t) represent noise in a temperature limited current I » we have 
the well known result 


(A. )=el,/r 
and so 
(13): 
Equation (13) then gives 
w(f)=4nfel, ih BIN emyprar. ic A be . (14), 


0 
In evaluating the integral we remember that in practice we evaluate 


ia w(f)df over some finite range $f. 
af f+ oF rf 
| wihdt= | w(aydg— | wg) dg... . (15) 


0 0 


Now 
rf ap 20 
| w(g) dg =2el, | 2ag i sin 2g dr dq. 
0 0 0 


Carrying out the integral over g first we obtain 


ig 2 fj 
| ee oar {ss zat _ feos =t a 
0 Jo QaT Ct ‘ 
Integrating the second term by parts gives: 


f 7 ee) 
| w(g)dg= — 21] 2" | HE e a en OE): 
0 TT 0 


Thus in this case we obtain: 


[wf df =2elo(f + of) —2elaf = 2eLe8f see tect (LL) 


the well known result for temperature limited shot noise. 
We may summarize this result as follows. If 


b(r) = Kr Mia tr wee Ao eee (AS ap 
olen eee ae, OF ee (LSD) 
REFERENCE 


MacDona cp, D. K. C., 1949, Phil. Mag., 40, 561. 


912 Correspondence 


The Rotation of the First Russian Earth-Satellite 


By J. H. Taomsony 
Cavendish Laboratory, Cambridge 


[Received June 12, 1958] 


Many workers have reported the presence of complicated fading in the 
signals from the first Russian earth-satellite. In a previous communica- 
tion (Staff of the Mullard Radio Astronomy Observatory, 1957) it was 
pointed out that this fading showed both that the satellite was rotating 
at about 7 rev/min, and that the plane of polarization of the radiation 
was rotated by the Faraday effect in the ionosphere. 

Records of fading during transits over Cambridge at two different 
heights have now been examined in more detail. During evening transits 
when the satellite was at a height of about 200 km and below most of the 
ionosphere, the Faraday rotation expected was about a quarter of a turn, 
so that the observed fading must be attributed to rotation of the satellite. 
Both at 20 Mc/s and at 40 Mc/s the records show fading at the same 
uniform frequency, as in fig. 1 (a). During morning transits when the 
satellite was at a height of about 500 km, close to the maximum of the 
F, layer, the expected frequency of Faraday fading at 20 Mc/s was about 
ten times that of the rotation. The fading caused by rotation then 
appeared as an envelope to the Faraday fading, producing records such 
as that-shown in fig. 1 (5). 

A plot of the frequency of fading caused by rotation, as a function of 
time, for the period of the 6th to the 24th of October is shown in fig. 1 (c). 
The points of this figure have been derived from records of the types 
discussed above. It may be shown that for a simple aerial system such 
as would be expected on the satellite, a frequency of fading equal to that 
of the rotation would be rare. The great majority of the fading attributed 
to rotation was at approximately 13 fades per minute; fading at 6-5 f/min 
was occasionally observed, but never at 3-5 or at 26 f/min. It is therefore 
concluded that the satellite was rotating at 6-5 rev/min. The frequency 
of rotational fading showed a decrease in the first few days from 14-0 f/min 
to 13-0 f/min. It then remained constant without detectable change at 
13-0 f/min for about 15 days until transmissions ceased. 

It has been pointed out that the rotation of a satellite with a continuous 
metal skin would be slowed by eddy currents generated by the earth’s 
magnetic field. The change from 14-0 to 13-0f/min. in the first few 
days is believed to be real, but was evidently not a continuous change, 
so it appears that no such slowing was detected. 
et ee 5 a a le Ee 

+ Now at Jodrell Bank Experimental Station, University of Manchester. 
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Much of the observed fading during transits when the satellite was 
in the F, layer was caused by the Faraday effect in the ionosphere. 
Fading records obtained during such transits may be analysed to give 
information about the total ionization up to the satellite (Aitchison 
et al. 1958). We shall here discuss how information about the orientation 
of the rotation axis of the satellite may be obtained from them. 


Fig. 1 


20 Me/s 


Neti Sa LL ORT a RE 


40 Me/s 


October 1957 
(c) 


(a) Recordings at 20 Mc/s and 40 Me/s during an evening transit, when the 
the satellite was below the ionosphere. 

(b) 40 Me/s recording, showing rotational fading superposed upon Faraday 
fading. 

(c) The period of rotational fading versus time. 


The transits at Cambridge at about 03.30 U.T. from the 16th to the 22nd 
of October were of this type, the satellite being at a height of about 
500 km. An almost complete set of fading records at both 20 and 40 Me/s, 
received with single horizontal dipoles, were obtained during this period. 
An examination of the 40 Me/s records (fig. 2) shows that a discontinuous 
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change in the fading rate takes place at about the time of closest approach 
in the records for the 18th and succeeding days. On the 17th this 
effect is just apparent at the end of the record, and there is no sign of 
it on the 16th. Before this transition the ratio of fading rates at 20 and 
40 Mc/s is approximately 4 : 1, as is expected if the fading is caused by 
the Faraday effect, but after the transition it falls to about 2°5 ; 1. 


Fig. 2 


Se 


October 1957 0 30 60 
: B seconds 
16 LEN See 
J Yes dex FENIAN es = 


i ss eee TMA ene 


5 EDR ONGOING NIG BO ER 


Fading of 40 Mc/s radiation during transits at about 03.30 U.T. 


This change is believed to be an effect of the aspect of the satellite 
relative to the observer. Consider for simplicity the aerial on the satellite 
to be a dipole, the line of which does not coincide with the rotation axis; 
then as the satellite rotates this generates a cone as shown in fig. 3. 
When the line of sight is in a position such as A, outside the cone the 
projection of the dipole on the plane normal to the line of sight is seen to 
oscillate, but not to rotate. When the line of sight is in a position such 
as B, within the cone, the projected dipole is seen to rotate. In case B 
the rotation of the line of polarization by the Faraday effect will have the 
rotation of the satellite added to or subtracted from it, and the observed 
fading rate at both frequencies will be increased or decreased by 13 f/min 
Hence it should be possible to recover the true Faraday fading rate from 
the latter part of records on which the transition occurs by subtracting 13 
from the fading rate observed. When this is done both to the 20 aa 
40 Mc/s fading rates, the ratio of the corrected fading rates becomes 
approximately 4:1, as it was before the transition, i 
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Fig. 3 


A 


Cone generated by rotation of a dipole. From directions such as A the dipole 
appears to oscillate, whereas from B it is seen to rotate. 


Fig. 4 


Tracks of the satellite in the region of Cambridge for the transits at about 
03.30 U.T. from October 16th to 22nd. 
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Figure 4 is a map of the area round Cambridge, on which the paths 
of the satellite are shown for the 03.30 U.T. transits on the marked dates. 
The periods of observation are shown by the heavy lines, and the periods 
after the transition are indicated by the zig-zags. That the change of 
aspect occurred at almost the same position relative to Cambridge on 
each of five successive days, each with 15 orbits between, suggests that 
the rotation axis of the satellite was in a constant direction. If this 
were so, as the path changed from day to day, Cambridge may be regarded 
as making cuts across the cone discussed above. If fig. 4 is inverted it 
may be regarded as showing the paths of Cambridge relative to a fixed 
satellite at C': the dotted line then marks the intersection of the cone with 
the surface of the earth. On this assumption both the cone angle and the 
direction of the rotation axis have been calculated. The semi-angle 
of the cone is about 45°, and the direction of the rotation axis is close to 
the direction of the orbit at perigee, being at 23° to it. 

Published photographs of the satellite before launching show four 
aerial rods trailing. Rotation of the satellite about its axis of symmetry 
would generate a cone of semi-angle 45°, which would have the same 
polarization properties as the cone discussed above. Hence the observa- 
tions would be consistent with the satellite rotating about its axis of 
symmetry and roughly lining up with aerials trailing each time it passed 
through denser air close to perigee, and maintaining this direction of 
rotation axis round the rest of its orbit. 
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Measurements of Thermoelectricity Below 1°k.—II 


By D. K. C. MacDonatp, W. B. Pearson and I. M. Temeieron 


Division of Pure Physics, National Research Council, Ottawa, Canada 


[Received June 17, 1958] 


In a previous letter to this Journal (MacDonald et al. 1958 a) we presented 
some measurements of the absolute thermoelectric force of gold and silver 
below 1°k. We have now completed an initial survey of the alkali metals 
and wish to publish some preliminary results which appear to be 
interesting. 

A detailed study of the absolute thermoelectric power of the alkali 
metals between 2° and 20°K (MacDonald et al. 1958b) showed clearly 
that, even as low as 2°K, the thermoelectric power did not conform readily 
to theoretical expectations. This disagreement was a major stimulus for 
extending the work below 1°x. First let us note that we find that 
impurities can have a significant effect on the thermoelectric power even at 
these very low temperatures. The results reported here relate to specimens 
of the highest purity which were available to us for this series of measure- 
ments. (The residual resistance-ratios referred to room temperature 
mereras tjollows: Li: 2-5x102; Na: 4:7x 10-4: Ky 25x 10-32 Cs: 
6-5 x 10-3.) 

Our measurements have shown that the thermoelectric power of lithium 
continues positive down to the lowest temperature reached (perhaps 
0-05°K) and is approximately linear with temperature. We find that the 
thermoelectric power of the sample of caesium also remains positive below 
2°x and varies linearly with temperature below about 1°K. Sodium and 
potassium however present an interesting contrast, and more detailed 
results on potassium are shown in figs. 1(a) and 1 (0). 

Wilson (1953) gives 


(1) 


Belo 


as the limiting theoretical expression for absolute thermoelectric power, S, 
at low temperatures. Taking ¢)=3-l16ev for sodium, and 2-06ev for 
potassium (cf. e.g. Mott and Jones 1936) eqn. (1) intersects the observed 
data around 2°x. However we find that the magnitude of the observed 
thermoelectric power decreases more rapidly than the predicted linear 
variation with 7’ down to the lowest temperatures reached. The observed 
dependence on 7’ appears to be more or less quadratic in contrast to eqn. (1), 
and indeed rather closely so in the case of potassium (see fig. 1 (6)). This 
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¢ E > 
behaviour is very puzzling; it is to be remembered that see a 
effects should be relatively small in sodium and potassium at these very 
* ° 
low temperatures, certainly below 0-5°kK. 


Fig. 1 (a) 
T (®K) 


2x10-8 


3x10-8 
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ABSOLUTE THERMOELECTRIC POWER, S (V/°C) 


7x10-8 


Absolute thermoelectric power of potassium. (O: Present experiments. 
A: Data from MacDonald etal. 1958b. ———-—-— : Calculated from eqn. (1). 


These measurements should be regarded as a preliminary survey. 
Recognizing the particular importance of the temperature-dependence, 
for comparison with theory, we shall make more detailed measurements 


using alternative techniques, and we intend to publish the results more 
fully later. 
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Fig. 1 (6) 
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Absolute thermoelectric power of potassium. (Logarithmic scales ; full line 
drawn to show 7” dependence.) 
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REVIEWS OF BOOKS 
Geometry of Hinstein’s Unified Field Theory. By VActAv Hiavaty. (Groningen 

Holland: P. Noordhoff Ltd.) [Pp. 341.] £3 12s. 6d. 

Tue struggle to find a unified theory of gravitation and electromagnetism 
presents one of the most curious spectacles in the development of modern 
science—the majority of physicists looking on with amused contempt at the 
efforts of some of the most distinguished mathematicians and physicists of the 
twentieth century—men such as Cartan, Weyl and Veblen; Einstein, 
Schrédinger and Pauli. Einstein, in particular, ignored the general scepticism 
and devoted a major part of his life from 1923 to the search for the ‘ total field ° 
and the relations that it satisfies. His last theory, and the one by which he 
set most store, is based on a non-symmetric tensor g;;. It is this theory which 
is the subject of Professor Hlavaty’s book. 

The author’s aim is to ‘provide a detailed geometrical background for 
physical applications of the theory’. He studies in particular (i) the algebraic 
properties of a non-symmetric tensor (ii) the equation which replaces the 
Christoffel relation of Riemannian geometry, and (iii) some general mathematical 
properties of the field equations. Other mathematical developments are 
contained in two appendices on spinor algebra and line geometry. 

Physicists will naturally be most interested in the chapter entitled “ Physical 
Interpretation of Field Equations’. In this chapter the author attempts to 
decide which of the many symmetric and skew tensors thrown up by the theory 
should be identified with the gravitational and electromagnetic fields. His 
conclusions should be treated with reserve, since there is as yet no general 
agreement as to the correct identification, or even whether such an identifica- 
tion is possible. 

Because of the limits which the author sets himself, his book is not a balanced 
account of Einstein’s theory. For such an account one must turn to Madame 
Tonnelat’s “‘ La Théorie du Champ Unifié d’Einstein et Quelques—une de ses 
Développements ”, which lays greater emphasis on physical considerations. 
Final judgment on the theory itself must be deferred until its implications 
have been more fully explored. In this difficult task, Hlavaty’s detailed 
analysis may prove to be of value. D.W.S. 


Matrix Theory for Physicists. By J. Heapinc. (London: Longmans Green 
and Co.) [Pp. 242.] 35s. 

WHEN does mathematics cease to be higher ? Problems involving simultaneous 
equations are familiar to every physicist, and he knows how to solve them. 
He may even invent a condensed notation to represent them. The prose he 
is talking is matrix algebra. The general notions of this subject—trans- 
formation, symmetry, inverse matrix, determinant, principal axis, character- 
istic root, eigenvector—have as great a role in physics nowadays as the notions 
of analysis—function, derivative, integral, Fourier component, divergence, 
etc. To understand them requires no greater basic mathematical equipment 
nor more strenuous intellectual effort. This book sets out, as simply as possible, 
the elements of the formalism, without laboriously precise proofs, and then 
works through various general fields of physics—geometry, dynamics, electro- 
magnetism, special relativity and quantum theory—with their aid. If a 
student already knows his physics he should have no difficulty in reading this 
book and acquiring a feel for a great unifying mathematical technique or, at 
worst, a facility with a fashionable jargon. It seems a readable account 
enlivened with plenty of examples. Only the chapter on quantum theory is a 
bit difficult ; Heisenberg’s matrix mechanics is a trifle sophisticated, and 
students usually find an approach through the wave equation more natural. 
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R. PHILLIPS and N. C. WELSH Phil. Mag. Ser. 8, Vol. 3, Pl. 44. 


Equiaxed subgrains together with smaller elongated cells produced by etching 
in hydrofluoric acid. x11 000. 


Subgrain boundaries showing the individual dislocations at B and dislocations 
together with electron-optical fringes at CC. x 130 000. 


R. PHILLIPS and N. C. WELSH Phil. Mag. Ser. 8, Vol. 3, Pl. 45. 


Equiaxed hydrofluoric acid etching structure and larger equiaxed subgrains 
< 20 000. 


Etching structure produced by Alcoa Bright Dip R5 and subgrain boundaries. 
< 30 000. 


R. PHILLIPS and N. C. WELSH Phil. Mag. Ser. 8, Vol. 3, Pl. 46. 


Lack of sharp detail in the etching structure cell boundaries compared with 
the fringes and dislocations in sub-boundaries. x 160 000. 


Fig. 6 


Cell boundaries remaining at D while neighbouring foil has been completely 
etched away. X11 000. 
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(a) Low magnification micrograph of polycrystalline Mg, deformed 0:8% in 
tension, showing large transparent areas. (Magn. x1100, 100 kv.) 
The curved lines such as at A are extinction contours ; the very fine 
short lines near B are dislocations ; the grain boundaries are clearly 
visible. (b) Moving dislocation loops in same specimen ; note the 
cusps in some of the loops. (Magn. x29 200, 100 kv.) 
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Fig. 3 


(a) Stacking faults in polycrystalline Co annealed at 900°c and furnace cooled 
to room temperature. Some partial dislocations can be seen where 
stacking faults terminate. (For detailed accounts of stacking faults 
see Whelan et al. 1957, Whelan and Hirsch 1957 a, b.) (Magn. x 37 500, 
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HEATHER M. TOMLINSON Phil. Mag. Ser. 8, Vol. 3, Pl. 48. 


(a) Irregular arrangement of dislocations in polycrystalline Cu, deformed 
6-7%, in tension at room temperature. (Magn. x37500, 80 kv.) 
(>) Trace left after passage of dislocation showing cross-slip in a specimen 
of polycrystalline Cu, deformed 6-7°% in tension at room temperature. 
(For explanation of dislocation trace see Whelan et al. 1957.) (Magn. 
x 37 500, 80 kv.) 


Fig. 5 


(a) Dislocations arranged in poorly defined sub-boundaries in polycrystalline 
ED : a Sis 

Ni, deformed 25% by rolling at room temperature. (Magn. x37 500, 
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Trregular dislocation arrangement in a Dislocations arranged on slip planes 
single crystal of Al deformed 2°% in polycrystalline «-brass, deformed 1‘ 
tension at room temperature. A wide by rolling at room  temperatur 
(light) dislocation trace can also be Small regions of stacking faults ce 
seen. (Magn. x 40 000, 80 ky.) be seen in various areas. (Mag 


x 40 000, 100 kv.) 


(a) Irregular arrangement of dislocations in polycrystalline Fe, deformed 
7 0/ es Gea / 
nee by rolling at room temperature. (Magn. x40000, 100 kv.) 
(b) Curved trace left after passage of dislocation, showing apparently 
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(0) 
Microstructure of specimens compressed statically 4%, (x 1400). (a) Deformed 
at +15°c. (b) Deformed at —183°c. 
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(0) 
pies ute Oooe eras (anvil bar records). Calibration lines correspond 
to stresses of 0, 14, 19, 28, 56, 80, 112, and 140 x 10? Ib/in2 (a) Te ; 
; : 3 


temperature : + 15°c; Impact velocity : 430 in/sec. (6) Test tempera 
ture: —41°c; Impact velocity : 470 in/sec. 7 ae Later 
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Fig. 6 
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(a) 


(0) 
Stress-time oscillograms (weighbar records) obtained from tests on high- 
tensile steeldummy. Impact velocity : 238 in/sec. Test temperatures : 
(a) +15°o, (b) —41°c. 
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ye ees in polyerystalline aluminium quenched from ~600°c into 
ce ee (specimen X), _ Note that some of the loops are in the form of 
parallelograms, A, or distorted hexagons, B. Mag. x 55 000 
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i i i i ing wide distribution 
High magnification micrograph from specimen X showing wi 
i of Te eet. Note also the parallelogram shaped loop, A. 
Mag. x 111 000. 
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Fig. 3 (a) Fig. 3 (b) 
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3. Dislocation loops in specimen Y. The orientation of the specimen 
is such that the plane of the foil is several degrees away from a (111) 
plane. The hexagon at A lies apparently on this (111) plane and its 
sides are parallel to [110] directions. The appearance of loops B, and 
C suggests that they he on another (111) plane (plane e¢ on fig. 5 (a)) not 
in the plane of the foil, and loops B, probably lie on plane b (fig. 5 (a)). 
Figure 3 (a) shows a slip trace of plane b (fig. 5 (a)), which has practically 
disappeared on fig. 3 (b), taken subsequently. The third photograph 
in this sequence (fig. 3 (c)) shows that loops A and C have vanished, 
giving rise to the slip traces which can be clearly seen. See text for 
interpretation of ship movements. Note that loop V has moved in the 
interval between the exposures of fig. 3 (b) and 3 (c). Mag. x 52 000. 


4. Electron transmission micrograph of furnace cooled specimen of 
aluminium. The sharp bands are extinction contours ; the kinks (A’) in 
the extinction contours are produced bv dislocation lines. Mao x 298 ONO. 
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Fig. 6 
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(a) Double loop in a quenched specimen. Mag. x 10000. (6) Complex loop 
and a possible case of a spiral dislocation (W). Mag. x 57000.  (c) 
Distorted hexagons, B, and a possible example of a spiral dislocation 
at W. Mag. x 60 000. 
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Denuded zone near a grain boundary. Note also the parallelogram shaped 
loops at A and the more or less distorted hexagons at B. Mag. x 44 000. 


P. B. HIRSCH et al. Phil. Mag. Ser. 8, Vol. 3, Pl. 58, 
Fig. 8 


Complex dislocation line and surrounding zone denuded of dislocation loops. 
It is thought that the complex shape of the dislocation is the result of 
climb of various sections and of the coalescence of the dislocation with 
loops. Note also a parallelogram at A and distorted hexagons at B. 


Mag. x 45 000. 
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Fig. 9 
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Large loops probably formed by junction of several loops. Note the virtual 
absence of small loops inside the large loops. The loops marked Y 
have shapes similar to that shown in fig. 10. Note also the complicated 
zig-zag dislocations at Z. The appearance of these dislocations at points 


marked X suggests that they have reacted with dislocation loops. 
Mag. x 40 000. 


